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Improvements  to  the  Automated  Real-Time  Tidal  Elevation  System,  TR  DRP-93-2 


ISSUE:  As  presently  implemented,  the  Auto¬ 
mated  Real-Time  Tidal  Elevation  System 
(ARTTES)  uses  a  completely  deterministic 
procedure  for  predicting  the  tide.  When  the 
difference  between  the  predicted  tide  and  the 
observed  level  exceeds  some  predetermined 
value,  user  access  to  the  system  is  denied.  Pe¬ 
riods  of  user-access  denial  usually  occur  when 
sea  conditions  are  unfavorable  or  unsuitable 
for  survey  operations.  Dredging  operations, 
however,  may  continue  in  weather  conditions 
unsuitable  for  survey  operations.  Because 
ARTTES  is  sometimes  used  by  dredge  opera¬ 
tors,  it  is  desirable  to  minimize  periods  of 
user-access  denial. 

RESEARCH:  A  hybrid  approach  was  used 
to  improve  the  ARTTES  tide  predictions: 
“hybrid”  in  that  the  astronomically  forced 
component  of  the  tide  was  analyzed  and  pre¬ 
dicted  in  the  traditional  manner  via  harmonic 
analysis  while  the  short-term  meteorological 
component  (the  residual)  was  analyzed  and 
predicted  via  statistical  techniques.  Kriging 
and  autoregression  statistical  techniques  were 
investigated. 


SUMMARY:  The  two  statistical  techniques 
were  applied  to  improve  short-term  tide  pre¬ 
diction.  Evaluations  using  both  synthetic  and 
prototype  data  indicated  that  although  Kriging 
has  a  slightly  larger  root  mean  square  error 
than  the  autoregression  technique,  it  gives  an 
unbiased  estimate  and  the  autoregression  tech¬ 
nique  shows  bias. 

Because  dredges  are  costly  to  operate,  data-ac- 
cess  denial  should  only  occur  when  the  poten¬ 
tial  consequences  warrant.  The  short-term 
tide  predictions  provide  an  objective  basis  for 
a  decision  by  the  dredging  supervisor  to  use 
or  ignore  the  ARTTES  during  periods  of  de¬ 
graded  data. 

AVAILABILITY  OF  REPORT:  The  report 
is  available  through  the  Interlibrary  Loan  Ser¬ 
vice  from  the  U.S.  Army  Engineer  Waterways 
Experiment  Station  (WES)  Library,  telephone 
number  (601 )  634-2355.  National  Technical 
Information  Service  report  numbers  may  be 
requested  from  WES  Librarians.  To  purchase 
a  copy  of  the  report,  call  NTIS  at  (703)  487- 
4780. 
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1  Introduction 


Procedures  for  increasing  the  availability  of  real-time  tide  data  broadcast 
by  Automated  Real-Time  Tidal  Elevation  Systems  (ARTTES)  (Lillycrop  et  al. 
1991)  are  described  herein.  The  measured  water-level  data,  which  are 
broadcast  to  user  vessels,  are  compared  to  the  predicted  tide  at  the  water-level 
sensor.  As  presently  implemented,  ARTTES  uses  a  completely  deterministic 
procedure  for  predicting  the  tide.  When  the  difference  b^een  the  predicted 
tide  and  observed  level  exceeds  some  predetermined  value,  user  access  to  the 
system  is  denied.  The  difference  between  measured  and  predicted  values  (at 
the  measurement  site)  is  an  indicator  of  the  validity  of  the  assumption  that  the 
astronomical  tide  component  is  dominant.  Dominance  is  defined  herein  to 
mean  that  the  astronomical  component  of  the  water-level  variation  is 
approximately  one  order  of  magnitude  greater  than  the  sum  of  the  remaining 
components.  Zero  or  small  (a  few  centimeters)  differences  mean  that  the 
assumption  is  robustly  met  within  the  entire  specified  area  of  operation  and 
that  users  can  be  confident  of  the  broadcast  system  values.  When  differences 
between  the  predicted  and  observed  values  increase,  as  is  typically  the  case 
during  periods  of  disturbed  weather,  the  assumption  of  dominance  of  the 
astronomical  tide  component  grows  increasingly  weak.  At  some  difference 
between  the  predicted  tide  and  observed  water  level,  die  assumption  of  tidal 
dominance  is  considered  to  be  invalid  for  the  desired  degree  of  water-level 
accuracy  over  the  area  of  operation. 

Periods  of  user-access  denial  usually  occur  when  sea  conditions  are 
unfavorable  or  unsuitable  for  survey  operations.  Dredging  operations,  how¬ 
ever,  may  continue  in  weather  conditions  unsuitable  for  survey  operations. 
Because  ARTTES  is  sometimes  used  by  dredge  operators,  it  is  desirable  to 
minimize  periods  of  user-access  denial. 

The  present  procedure  for  predicting  the  tide  at  the  sensor  location  is  based 
upon  the  harmonic  method  of  tide  analysis  and  prediction  (Schureman  1985). 
Water-level  data  at  the  desired  location  are  acquired  for  several  months,  ana¬ 
lyzed  to  determine  the  amplitude  and  phase  of  the  dominant  tidal  constituents, 
then  used  to  predict  the  tide.  The  data  required  for  tide  prediction  are 
obtained  prior  to  initial  operation  of  an  ARTTES.  Although  water-level  data 
at  the  sensor  site  are  recorded  and  available  once  an  ARTTES  begins  opera¬ 
tion,  analysis  of  these  data  by  the  harmonic  method  and  incorporation  of  the 
results  into  the  prediction  routine  presentiy  being  used  typically  result  in  little 
improvement  in  the  predictions.  Ihe  tack  of  improvement  in  tide  prediction 
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steins  from  the  fundamental  assumption  undo'iying  the  harmonic  m^hod, 
mundy  that  the  tide  is  described  by  a  sum  of  sine  or  o}sine  functions  that 
rqtresent  die  astronomically  forced  component  (Godin  1972).  However,  a 
portion  of  die  water-level  record  is  not  described  by  this  sum.  That  portion  is 
termed  the  residual.  Widiin  the  context  of  this  rqiort,  most  of  the  residual  is 
the  result  of  short-tmn  meteorological  effects. 

The  approach  used  in  improving  the  ARTTES  tide  predictions  may  be 
described  as  "hybrid"  in  that  the  astronomically  forced  component  of  die  tide 
is  analyzed  and  predicted  in  the  traditional  manner  (via  harmonic  analysis), 
n^ile  foe  short-term  meteorological  componeit  (foe  residual)  is  analy^  and 
predicted  via  statisdcal  techniques.  The  statistical  techniques  investigated  are 
conventionally  called  "Kriging”  and  autoregression.  Because  tidal  harmonic 
analysis  is  wdl-documented  and  its  use  well-established  (e.g.,  Godin  1972; 
Schureman  1985),  no  fiirfoer  description  is  contained  herein.  The  chapters  in 
this  report  are  arranged  to  describe  (a)  the  problem  statement,  (b)  an  overview 
of  foe  Kriging  and  autoregressive  methods,  (c)  implementation  of  the  methods, 
(d)  results,  including  a  case  application,  and  (e)  a  summary  of  conclusions. 
The  Appendices  document  one  of  foe  case  applications.  Appendix  A  is  a 
synthesized  time  series  of  known  statistical  properties  that  simulates  a 
residual  time  series.  Appendix  B  is  a  listing  of  the  computer  program  used 
prqiare  both  Kriged  and  autoregessive  estimates  of  foe  residual.  Appoidix 
is  a  tabulation  of  foe  computed  results  using  foe  Kriging  estimate  method. 
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2  Background 


If  z(t)  represents  the  observed  water  levd  at  any  tinoe  t,  the  water  level  z, 
including  the  harmonic  representation  of  the  tide,  may  be  written  as 

^  cos(«*r  -  ♦  V(r)  (I) 

where 

Zo  =  mean  water  level  (or  other  datum) 

H  =  number  of  significant  harmonic  components 
A^,  at  =  amplitude  and  phase  lag,  respectively,  of  kth  component 
Ut  =  frequency  present  in  tide-generating  force 
V(f)  =  contribution  of  physical  forcing  other  than  astronomical 
The  residual  is  then 

V(r)  =  z(r)  -  ^0  ■  53 


The  observed  water  level  z(r)  is  normally  logged  by  the  ARTTES  at  6-min 
intervals  (10  values  per  hour).  The  residual  V(r),  the  time  series  upon  which 
the  statistical  procedure  operates,  is  obtained  by  subtracting  the  pr^icted  tide 
at  concurrent  times. 

Consider  a  discrete  time  series  of  residual  values 

{V,;  i  =  l,2,3,...,iV}  (3) 
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Having 


W  *  "i 

4*’.  -  <‘i) 


w 


whoe 

£  [  ]  =  is  die  expected  value  of  the  quantity  in  brackets 
fit  -  mean  of  series  V, 

Q  =  covariance  of  series  v;  at  log  k 
That  is, 

W,.V,-  (5) 

is  taken  as  being  covariance  stationary  with  mean  zero,  where  W,  is  the 
residual  with  its  mean  removed. 

The  purpose  of  this  study  is  to  detomine  the  values  of 

from  measured  values  of 
(V.-,.  V,.„  V.., . V.,). 

The  data  values  {V^  n-m-L  £  i  £  n>l}  for  a  user-selected  extmsion 
L  ^  0  ,  will  be  fitted  with  a  straight  line  to  estimate  aI(  as  a  linear  trend 


fif  -  a  bi 


(6) 


where 

fi  =  is  the  estimated  local  mean 
a,  b  =  parameters  calculated  at  each  time-step 
*  =  a  parameter  to  be  estimated 

The  covariance  function  will  be  estimated  as  mean  lag  products  of  deviations 
from  the  local  mean 

C.  -  ^  {'',-<■-*>1  {‘'w-a-Mi**))  < 
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3  Overview  of  Methods 


Two  methods  were  used  to  improve  die  ARTTES  tide  predictions:  Kriging 
and  autoregressive  estimation.  Kriging  is  based  on  the  minimization  of 
expected  square  error,  subject  to  the  constraint  that  the  estimation  is  unbiased. 
Autoregression  is  based  on  minimization  of  expected  squared  error  only. 

Both  methods  will  assume  that  the  estimate  is  covariance-stationary. 


Kriging  Procedure 


Here,  is  tbe  estimate  of  the  residual  at  time-st^  n-k-k  and  is  predicted 
as  a  linear  combination  of  the  data  for  time-stqis  preceding  time-stqi  n. 


<•1 


^1 


(9) 


So  as  to  minimize  Q,  the  expected  value  of  die  error  squared  between  the 
actual  and  estimated  residuals 

where 


Q  = 


Subject  to  the  unbiased  constraint 


(10) 


(11) 


which  requires  that  the  expected  value  of  the  estimator  be  equal  to  the 
expected  value  of  the  residual  at  time  step  n+k 

Since, 

0,  -  £[v;]  02) 
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the  "unbiased”  constraint  can  be  written 


Similarly,  ^  can  be  expanded  to 


Since  the  "unbiased”  conditions  are  going  to  be  imposed,  the  expr«iSsion  does 
not  change  if  the  "unbiased*  constraint  equation  (equal  to  zero)  is  inserted 
inside  the  square 


Q-E\  -  M.,.)  -  E 


If  the  algebra  of  squaring  is  completed 


-2E  -0(1'. 
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The  term  with  the  summation  sign  £  squared  can  be  expanded  as  follows: 


Nl 

OT 

•  < 

m  m 

j-t 

* 

J»1  J-\ 


(18) 


Hence,  the  expected  square  mor  tarn  ‘Q*  may  be  expressed  as 

a  -  e[(v.^-m,^)*] 

♦EE  (19) 

i-i  j-i 

$•1 

By  definition  of  the  covariance  function,  and  the  assumption  of  covariance 
stationarity  after  removing  the  trend 

C.  -  IM'',- 

It  follows  that 


(20) 


J«1  y-1  i»l 


(21) 


The  determination  of  die  coefficients  {a;}  for  Kriging  dien  involves  solving 
for  the  {ai}  ,  which  minimize  Q  subject  to  the  constraint  of  being  unbiased. 

By  the  method  of  Lagrangian  multipliers,  this  is  achieved  by  solving  the 
system  of  equations  given  by 
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* 

Ml 

■ 

Q-2\ 

. 

i-l 

d 

da. 


s  »  1,2,3,. ..w 

and  by  the  unbiased  constraint  Equation  (13) 


i«l 


(22) 


Performing  the  differentiation  gives 


2  a/:,.;  -  2  -  2V,.««-i  -  0 

>1 


a3) 


If  die  equations  are  listed  for  s  —  1,2,3,. ..,m,  ,  using  =  C^,  they  become 


p.^"l  ♦  p.— A  ••  p..l0«  *  K*k 


(24) 


The  first  m  equations  can  be  divided  by  Cq  to  express  the  relations  in  terms 
of  the  correlation  coefficients  p^. 


fl,  +Pia2+...+p„_ia„  -  (X/C^a) 

Pia,+a2  +  ...+p„.2a„+(X/q,)p«.i  -  p*^„.i 


(25) 
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This  cao  be  written  in  matrix  form  as 


• 

1  Pi  Pi  •  •  •  Pai-1  Pu-m 

•  • 

Pk*m 

Pi  1  Pi  ..  .  Pg,.2  Pm-m*l 

Pk*m-1 

P2  Pi  1  •  •  •  Pm-3  Pm~m*2 

^3 

m 

Pk*tm-2 

Pm-1  Pm-2  Pm~3  *  •  *  ^  P«-l 

Pk*l 

Pm-m  Pm-m-*1  Pii-m*2  •  •  •  Pm-1  ® 

x/q,_ 

Pm*k 

This  can  be  further  simplified  if  the  following  definitions  are  introduced 


1 

Pi 

Pi 

•••  P«-l 

Pi 

1 

Pi 

•  •  •  Pbi-2 

P2 

Pi 

1 

•  •  •  Pm-3 

P«-l 

Pm-2 

Pm-3 

...  1 

(27) 


Po 


Pk*m 

Pk*m~l 

Pk*m~2 


Pk*l 


(28) 


(29) 
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The  Kriging  equations  may  be  written  in  t^ms  of  these  arrays  as 


(30) 


(31) 


The  mean-square  error  of  estimate  is  the  value  of  Q  with  the  Oj,  which  are 
solutions  to  the  Kriging  equations  inserted.  If  (21)  is  multipli^  by  q,  and 
summed  over  s,  the  result  is  (dividing  by  2) 


<•1  >•!  l-l  1-1 


(32) 


By  the  "unbiased"  constraint,  the  summation  involved  may  be  written 
•  Consequently,  the  last  equation  becomes 


m  m 


m 


E  E  ~  E  ~t*l 

»-l  J‘l  <-I 


(33) 


This  relation  bolds  for  solutions  to  foe  Kriging  equations.  If  this  relation  is 
substituted  into  Q, 


Q  = 


i-J 


(34) 


The  root-mean-square  (rms)  error  of  estimate  for  foe  Kriging  method 
presented  is 


d  = 

B 


^0  ~  E  ~ 


1-1 


(35) 


where  foe  {Oj}  are  foe  solm.*....'  t-  he  Kriging  equations. 
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Autoregressive  Procedure 


The  autoregressive  technique  is  based  on  the  assumed  stnumiral  model 
equation  that,  for  all  j, 


♦  t. 


!•! 


(36) 


with  Wj  as  defined  previously  as  per  Equation  S  and 

where  the  sequence  is  tak«i  indq>endently,  with  mean  zero  and  constant 
variance.  The  {a,}  are  parameters  to  be  fit  via  the  methodology  to  be 
presented  below.  This  means  that  with  the  stationary  condition,  for  (  <  J 


Cj-I  -  E[V,Vl 


=  E 


W. 


E  *  <; 
#•1 


(37) 


Since  f  <  y  ,  the  nror  in  is  independent  of  W,  and  £[W,ej]  =  0. 
Consequently, 


Cy,  -  E  oflW.Wr--,) 

<•1 


Hie  model  equation  for  is 


W.  =  a. 


+  ♦ ...  +  ♦  c. 


(38) 


(39) 


Ify  =  n  whUe  n-m  ^  f  <  n  ,  the  covariance  becomes 

m 

~  X)  (40) 

/-I 

Now,  if  successively,  (  =  n-m.  n-m+l,  n-m+2, .....  n-l  ,  the  resulting 
system  of  equations  is  (using  =  Cj) 
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(41) 


*aS,.2  =  c.., 

+a.c..,  =  C..J 


®l^ip-l'*'®2^»-2'*‘**3^«-3'*'  ••  ■*‘**«Q  “  ^1 

The  system  of  equations  can  be  expressed  in  matrix  form,  after  division  with 
Q,  as 


■  1  p,  p,  ...  p..r 

’<*r 

p-* 

Pi  1  Pi  •••  P-j 

P-I 

Pi  Pi  1  •  •  P«-3 

= 

P-1 

P--I  P— 1  Pw-3  •••  1 

Pt . 

This  is  one  form  of  the  wdl*known  Yule-Walker  equations  for  the  coefficients 
{aj  to  predict  W,  from  {Wj;  n-m  n-l)  in  a  stationary  sequrace.  It  is 

convenient  to  express  die  prediction  of  W,  in  a  somewhat  r^undant  way. 


0  1  0  0  ...  0 

0  0  1  0  ...  0 

0  0  0  1  ...  0 


0  0  0  0...  1 

a,  flj  o,  a,  ... 


•  « 

w-1 

w,-. 

w.. 

w 

n 

This  can  be  written  as 


(43) 


(44) 


where  A  is  the  matrix,  and  the  bold  print  denotes  vectors.  The  actual  predic¬ 
tion  of  is  given  by  the  last  component  of  the  vector  .  The  matrix  A  is 
a  shift  operator  that  moves  the  sequence  up  by  one  step.  The  prediction  of 
is  the  last  component  of  W,+y ,  where 
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(45) 


W,,i  -  AW^  - 


Similarly, 


-  AW„,  .  A^W,., 


Alternatively,  this  sequence  of  equstfions  can  be  used  to  predict 
^+2»‘"*^+«)  {^-m*  ^-m+j'  Also,  the  mean 

square  error  of  estimate  for  each  of  these  predictions  can  be  obtained  from  a 
modification  of  this  same  sequence  of  equations.  Since  =  0  (by 
subtraction  of  trend),  if  is  defined  as  the  covariance  matrix  of  , 


Here,  the  superscript  "T  denotes  the  matrix  transpose.  It  follows  that 


for  ^  =  0,1,2,3 .  The  variance  of  die  prediction  of  is  die  (m,m) 

element  of  Since 

CJ,  Cj  C2  ...  C„_i 

^  Q  •"  Q»-2 

C2  Cl  Cq  ...  C„.3 


Cm-l  ^m-2  C„_2  "•  Cq 


the  variance  of  each  prediction  can  be  obtained  as  the  (m,m)  -  element  of  each 
of  the  sequence  of  covariance  matrices 
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(50) 


C.,, 


An  algorithm  that  takes  advantage  of  the  structure  of  A  where 

0  1  0  0  ...  0 

0  0  1  0  ...  0 

0  0  0  1  ...  0 


0  0  0  0...  1 
a,  a,  a,  ...  a. 


(51) 


where  0  is  an  (m-1)  >  component  vector  of  zeros,  I  is  an  (m  - 1)  by  (m  - 1) 
identity  matrix,  and 


(52) 


Let  a  general  step  in  the  iterative  computation  of  k  =  0,1,23,—}  be 
r^resented  by 

B  =  ACA^  (53) 
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Fuidimaore,  let  C  be  partitioned  as 
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C 


C„ 

c 


(54) 


where  Q]  is  a  scalar,  C22  is  an  (m-1)  by  (m-1)  matrix,  and  C  is  an  (m-1)  - 
component  column  vector. 

That,  similarly,  1^  B  be  partitioned  as 


B 


Bu  B 


(55) 


where  fi||  is  an  (m  - 1)  by  (m - 1)  matrix,  B  is  an  (m - 1)  -  component 
column  vector,  and  B22  is  a  scalar.  Then, 


B 

0  7 

C.,  C' 

0*“  Oj 

B'  B^ 

a^  AJ 

c  q, 

/  A 

It  follows  that 

Bji  -  Cjj 

B  at 

Bjj  M  Cuaf+2a|A^C+A’tj2A 


(56) 


(57) 


(Note:  The  scalar  B^  is  the  variance  of  the  prediction  at  that  particular 
iteration.) 

The  algorithm  consists  of; 

(1)  Compute  a^C  and  store  in  B. 

(2)  Multiply  B  by  2.0  and  store  in  work  vector  D . 

(3)  Compute  vector  CjjA  and  add  into  both  B  and  D . 

(4)  Calculate  the  scalar  A^D  and  store  in  B22  ■ 
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(6)  Move  C22  into  upper  left  (m-1)  x  (m-1)  submatrix  of  new  covariance 
matrix. 

(7)  Store  B  into  upper  right  (m-1)  by  1  submatrix  (i.e.,  vector)  of  new 
covariance,  and  its  transposition  into  the  low«-  left  1  by  (m-1)  row 
vector. 

(8)  Store  in  the  (m,m)  position  of  the  new  covariance  matrix. 

This  sequence  of  computations  will  allow  the  new  covariance  matrix  to  be 
stored  into  the  same  memory  as  the  preceding  covariance  matrix.  The  rms 
estimate  of  error  of  prediction  is  given  by 


'it* 


<•1  i*i 


(58) 
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4  Implementation 


The  Low-Pass  Filter 


An  initial  low-pass  filter  is  applied  to  the  time  series  to  eliminate  any  high- 
frequency  noise  diat  may  remain  in  the  signal  after  its  earlier  processing.  The 
time  series  for  times  from  JO-LPREV  to  JO+M-1  (i.e.,  M+LPREV  =  MPLP 
time-stq)s)  is  first  transformed  with  the  fast  Fourier  transform  (FFT) 
algorithm  to  die  frequency  domain.  An  FFT  subroutine  based  on  MPLP  time- 
stq>s  being  expressible  as  a  product  of  powers  of  2,  3,  and  5  is  used.  The 
fre^ency-domain  version  of  the  time  series  imerval  is  that  multiplied,  for 
eadi  frequency  st^,  by  a  filter  fonaion  that  is  almost  1.0  for  frequencies 
between  -FLPASS  and  -f  FLPASS,  and  almost  zero  for  frequencies  with 
absolute  value  greater  than  FLPA^+ WIDTH.  A  smoofo  transition  within 
the  interval  (FLPASS,  FLPASS -f  WIDTH],  and  its  mirror  image  in  n^ative 
frequencies  is  achieved  with  a  cumulative  normal  probability  shape.  The 
accuracy  with  whidi  1.0  on  low  frequencies  and  0.0  on  high  frequencies  is 
attmned  is  controlled  with  the  input  variable  FRAC  in  die  subroutine  call.  If 
FRAC  s  0.01,  for  exanqile,  die  error  will  be  less  dian  1  percent.  Since  the 
frequency  filter  ftmction  is  symmetric  about  frequency  equal  zero,  the  filter 
will  have  zero  phase  shift.  All  frequ«icies  (FLPASS  and  WIDTH)  are 
expressed  in  cycles  per  time  unit,  rather  than  time-step  number.  Thus,  these 
values  can  be  entered  independendy  of  the  number  of  stqis  in  the  FFT  call. 


The  Computer  Program 


The  computer  program  PREDICT.FOR  assumes  input  of  a  time  series, 

[Vj] ,  for  j  »  ,  NJ  as  read  from  file  TIMESER.DAT.  A  maximum 

value  for  m  is  specified  by  M,  although  the  program  may  choose  to  use  a 
shorter  m  value  than  M.  A  constant  JB  is  input  in  parameter  statements, 
which  controls  the  time-st^  for  the  first  prediction,  at  stq>  JB-i-M.  Estimates 
are  then  made  for  each  tiroe-stqi  from  JB+M  to  the  end  of  die  sequmce.  An 
rms  error  of  estimate  is  provided  at  each  stqi. 
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After  the  application  of  the  low-pass  filter  on  the  interval  [JO-LPREV, 
JO+M-1]  at  each  prediction,  a  straight  line  is  fitted  to  the  data  for  time-steps 

JO  -  LPREV  ^  JO  +  M- I  (59) 

whae  the  prediction  is  being  made  for  step  JO  +  M.  This  straight  line  is 
subtracted  out 

Wj^Vj-a-bJ  (60) 

and  covariance  function  estimates  are  made  over  the  same  interval  of  time- 
stq>s. 

The  computer  program  searches  for  the  lag  just  before  the  lag  at  which  the 
covariance  Unction  has  its  first  down-crossing  of  zero.  The  covariance 
function  at  all  higher  lags  is  set  to  zero,  and  m  is  set  to  equal  the  terminating 
lag.  This  seems  to  give  good  estimates  that  can  be  obtained  quickly  on  the 
computer.  However,  the  code  is  easy  to  adjust,  should  other  choices  appear 
better  in  subsequent  use. 

The  parameter  IFLAG  controls  which  prediction  method  is  used.  If 
IFLAG  is  1,  then  only  Kriging  estimates  are  prepared.  If  IFLAG  =  2,  only 
the  autoregressive  me^ods  are  used.  If  IFLAG  =  3,  both  procedures  are 
calculated.  Several  other  flags  are  included  in  the  program  for  the 
convenience  of  the  user.  The  flag  ID(X  causes  a  variety  of  check  output  to 
be  printed  on  DEVICE  2  if  IDOC=  I .  If  IDOC=0  the  ^eck  ouqmt  is 
suppressed.  Another  flag,  I6YPASS,  controls  whether  or  not  detailed  predic¬ 
tion  output  is  printed  into  DEVICE  2.  If  IBYPASS=0,  the  material  is  printed 
(i.e.,  the  output  section  is  not  bypassed).  If  IBYPASS= 1,  the  seaion  is 
bypassed.  A  brief  listing  of  the  low-pass  Altered  series  and  the  shutdown 
decisions  is  output  onto  DEVICE  4. 

Three  appendices  at  the  end  of  this  report  list,  respectively,  the  input  file 
GARSHOR.DAT  (in  compressed  form),  the  computer  program 
PREDICT6.for,  and  a  lengthy  listing  of  the  output  results  from  DEVICE  4. 

In  the  output,  the  autoregressive  estimates  have  slightly  less  rms  error,  but 
are  biased,  with  a  bias  of 

— «  1  -SUM  AR  COEF  «  1  -0.80226 

true  value  (61) 

-  0.19774 

The  Kriging  estimates  are  constrained  to  be  unbiased,  but  as  a  result,  they 
have  a  slightly  larger  rms  error.  Both  estimates  appear  good  and  involve 
about  the  same  computational  difficulty. 
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Filter  Characteristics 


Filter  characteristics  were  studied  by  an  approximate  procedure.  Since 
both  die  Kriging  and  the  autoregressive  approaches  depend  on  the  actual 
muitipliCT  coefficimts  and  the  covariance  function  of  die  input  data,  general 
results  cannot  be  given.  However,  ^proximate  behavior  can  be  deduced 
from  an  assumption  that  the  coefficients  produce  behavior  about  like  a  convo¬ 
lution  with  a  half-tent  function  given  in  the  log  domain  as 


w(*) 


0  s  A  i 
otherwise 


(62) 


where 


A,  the  fimction  width 


The  Fourier  transform  of  w(h)  and  its  associated  amplification  and  phase 
fimctions  rqiresent  an  approximate  description  of  the  filter  characteristics. 
The  Fourier  transform  of  w(h)  is 


After  integration  by  parts  and  some  alg^ra,  this  reduces  to 

m 


(63) 


l-« 

X/Rq 


(64) 


The  real  and  imaginary  parts  of  W0  are 
[sinCicyAo)? 


Ifeim) 


(65) 


Im{W0] 


1  sinCn^Ag)  cos(nfli^ 


The  amplification  and  phase  fiinaions  of  the  filter  are  then 


(66) 
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AMm  -  ARe{W{f)}y*Um{wm^ 

PHAZ(fi  -  arcUualIm{W(/)}/lte{W(f)}] 

All  of  diese  are  functions  of flio-  Hence,  a  table  of  frequmcy  characteristics 
of  the  half-tent  filter  can  be  listed  versus  values  of  flio  ,  which  is 
dimensionless.  The  filter  has  an  amplification  of  1.0  at  f  -  0 ,  and  a  half¬ 
peak  drop-off  zt  fho  —  0.774S.  Since  /v>  in  the  data  used  is  around  8,  the 
half-peak  frequency  would  be  0.774S/8  =  0.097  cycles  per  time-stq>,  or  a 
corre^nding  wave  Imgth  of  10  time-st^s.  Wavelengdis  shorter  than  this 
are  b^g  suppressed  by  the  filter. 


(67) 

(68) 


Stationarity  of  Prediction  Sequence 


The  sequence  of  predictions  {^,j=n,n-i-l,n-i-2,...},  as  based  on  {Vi,i=n- 
m,n-m-f  l,n-m+2,...,n-l},  is  not  a  statistically  stationary  sequence  in  either 
Kriging  or  autoregressive  estimation.  The  nns  error  varies  with  time-stq). 
There  are  other  differences  dq>ending  on  the  specific  properties  of  each 
method.  It  is  particularly  revealing  to  examine  the  asymptotic  behavior  of  die 
sequences  for  each  method. 


Asymptotic  Behavior  of  Kriging 


The  basic  structure  of  the  Kriging  estimation  is  given  by  Equation  31.  If 
the  time-stqi  at  the  estimation  is  sufficiently  far  ahead  of  the  data  time-steps 
used  in  making  the  estimate  that  the  value  Aere  is  uncorrelated  with  the  data, 
then  p«  =  0,  and  the  solution  to  Equation  31  can  be  written  as 

b  = 


a 

R-^nJb 

0 

\/C,_ 

fi^R-^/b 

-i.m 

Hence,  the  Kriging  coefficients  and  Lagrangian  coefficient  are 
a  =  n^^R-^tUb 


(70) 
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and  the  prediction  and  mean-square-error  is 


1.0 


(71) 


ConsequoiUy,  the  Kriging  error  does  not  go  to  zero,  but  has  some  value  that 
dqtends  on  R.  The  asymptotic  estimate  of  V  converges  to  a  weighted  average 
of  die  data  values  times  the  mean  at  die  prediction  time-step. 


Asymptotic  Behavior  for  the 
Autoregressive  Estimates 


From  Equation  46,  asymptotic  b^vior  of  depends  on  the  absolute 
value  of  the  largest  eigenvalue  of  the  matrix  A.  This  can  be  less  dian  1.0, 
equal  to  1.0,  or  greata-  than  1.0  for  reasonable  choices  of  correlation 
coefficients.  (For  example,  in  an  m=2  autor^essive  sequmce  with  correla¬ 
tion  coefficient  at  lag  1  equal  to  0.5  and  at  lag  2  equal  to  0.1,  the  absolute 
value  of  the  largest  eigenvalue  is  0.4472.  If  die  lag  1  corrdation  coefficient  is 
dianged  to  0.9  while  the  lag  2  correlation  coefficient  is  left  at  0.1,  the 
absolute  value  of  the  largest  eigenvalue  is  1.2247.  If  the  absolute  value  of  the 
eigenvalue  is  less  than  1.0,  then  die  asymptotic  value  of  the  estimate  of 
will  tend  towards  zero  as  k  increases  indefinitely.  Similarly,  if  die  absolute 
value  of  the  largest  eigenvalue  is  greater  than  1.0,  then  the  estimate  will  be¬ 
come  indefinitely  large  as  k  increases. 
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5  Results  and  Conclusions 


Case  Applications 


Application  of  the  Kriging  or  autoregressive  procedures  requires  the  user 
to  specify  two  parameters  related  to  the  physical  situation.  The  first  is  the 
thrc^old  value  of  the  difference  between  observed  and  predicted  water  levels 
(the  residual)  beyond  which  ARTTES  data  may  be  degraded  (variable  THOLD 
in  Appendix  B).  The  second  is  the  permissible  duration  of  the  threshold 
exce^ance  (variable  NPF  in  App^ix  B).  The  threshold  value  is  chosen 
using  engineering  judgment  and  knowledge  of  residual  values  from  previous 
watN-level  monitoring  at  die  site.  The  goal  is  to  choose  a  value  such  that 
smaller  residuals  represent  meteorological  effects  that  are  large  enough  in 
scale  to  have  a  uniform  effect  over  the  specified  area  of  ARTTES  operation. 
The  assumption  of  spatial  uniformity  becomes  somewhat  questionable  for 
residuals  exceeding  the  threshold. 

In  die  earlier  non-ad^tive  procedure,  user  access  was  denied  whenever  the 
residual  miceeded  a  preset  limit,  even  if  die  threshold  was  just  slighUy 
exceeded.  The  Kriging  or  autoregressive  procedures  provide  an  estimate  of 
the  duration  of  a  specific  threshold  exceedance.  The  magnitude  of  the 
estimate  is  related  to  the  size  and  severity  of  the  event,  causing  die  deviation 
from  predicted  water  levels,  thereby  providing  an  objective  basis  for  a 
decision  to  use  or  ignore  the  system  during  periods  of  degraded  data. 

The  essence  of  this  suggested  new  approach  is  that  the  assumption  of 
spatial  uniformity  is  likely  to  be  met,  even  though  the  threshold  is  exceeded, 
if  data  from  the  recent  past  are  used  to  distinguish  short-duration  events  from 
long-duration  events,  which  alter  the  fundamental  physical  setting.  Alteration 
of  the  physical  setting  for  extended  periods  affects  the  site-specific  tidal  ampli¬ 
tude  and  phasing  relationships  upon  which  the  ARTTES  calculations  are 
based,  thereby  producing  potentially  erroneous  values.  The  duration  limit  is 
based  upon  knowledge  of  the  history  of  site  water-level  data  and  engineering 
judgment. 

By  using  combinations  of  levels  of  threshold  exceedance  and  duration, 
practical  aspects  of  the  dredging  process  may  be  taken  into  account.  Because 
dredges  are  costly  to  operate,  data  access  denial  should  occur  only  when  the 
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potential  consequraces  warrant.  In  other  words,  data  degraded  to  a  known 
degree  are  usually  preferable  to  no  data  at  all.  For  example,  the  dredging 
supovisor  may  dean  it  preferable  to  permit  a  relatively  high  level  of  thresh¬ 
old  exceedance  for  a  relatively  short  time;  say,  an  estimated  duration  of  1  hr, 
rather  dian  deny  access. 

Two  applications  were  made  of  the  Kriging  and  autoregressive  procedures. 
The  first  was  to  a  syndiesized  time  series  (Appendix  A)  of  known  statistical 
properties.  The  values  of  the  variables  THOLD  and  NPF  for  this  tqtplication 
wae  0.25  m  and  10  time-steps,  respectively.  Results  of  this  sq>plication  are 
given  in  Appendix  C.  In  general,  periods  of  access  denial  for  this  simulation 
are  not  found  until  foe  latter  part  of  foe  time  series,  first  at  time-stq)  1200  and 
again  at  time-step  1600.  Beginning  at  time-step  1637,  there  is  an  extended 
period  of  access  denial  that  continues  until  time-stq)  1683,  or  47  time-steps. 
Examination  of  foe  input  time  series  for  this  period  shows  residual  values 
exceeding  foe  threshold,  often  by  factors  of  two  or  more.  Moreover, 
comparison  of  foe  input  and  oufout  time  series  shows  that  short  periods  of 
direshold  exceedance  in  foe  input  time  series  are  smoothed  out  by  foe  low-pass 
filter  and  therefore  do  not  invoke  application  of  foe  Kriging  procedure. 

The  Kriging  and  autoregressive  procedures  also  were  ^plied  to  a  residual 
time  saies  acquired  at  an  operating  ARTTES  system.  This  time  series  is 
typical  of  periods  when  we^er  conditions  are  unsuitable  for  survey  vessel 
operation  but  when  dredges  were  likely  to  continue  operating.  The  threshold 
value  and  accq>table  duration  for  this  sq)plication  were  set  at  0.2  m  and  3  hr, 
respectively.  Results  of  foe  applied  Kriging  procedure  are  shown  in  Figure  1. 
The  autoregressive  procedure  produced  almost  identical  results.  Both  periods 
during  whidi  user  access  is  denied  under  foe  present  system  and  foe  Kriging 
procedure  are  illustrated.  Note  foe  shorter  period  of  user  denial  with  foe 
Kriging  procedure  when  compared  with  foe  existing  system.  As  with  foe  syn¬ 
thesized  time  series,  a  low^ass  fUter  was  ^plied  to  foe  input  time  series  to 
remove  high-frequency  components  before  applying  foe  Kriging  or 
autoregressive  procedures. 


Summary  and  Conclusions 


Objective  estimation  procedures  have  been  developed  to  augment  the 
present  method,  which  denies  users  access  to  ARTT^  data  during  periods  of 
poor  weather  conditions.  The  procedures,  termed  "kriging"  and 
"autoregression,"  provide  an  estimate  of  foe  duration  that  a  given  residual 
threshold  will  be  exceeded.  This  estimate  can  be  used  by  project  managers  as 
a  decision  basis  for  permitting  dredge  operators  continu^  access  to  the  system 
during  periods  of  possible  degraded  data. 

A  computer  code  is  provided  for  determination  of  Kriging  and  autoregres¬ 
sive  estimates  of  future  residual  values  (VJ  based  upon  residual  values  at  foe 
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Figure  1 .  Comparison  of  ARTTES  access  periods  for  ad  hoc  and  adaptive  procedures  for 
0.2-m  threshold 


preceding  tin)e-stq)s.  The  code  includes  a  low-pass  filter  with  user-specified 
cut-off  frequmcy  (FLPASS). 

Two  estimation  procedures  wwe  evaluated  during  this  study;  autoregres¬ 
sion  and  Kriging.  It  was  determined  that  although  the  Kriging  procedure  has 
a  slightly  larger  rms  error  than  the  autoregression  procedure,  the  Kriging 
estimate  is  unbiased  whereas  the  autoregression  estimate  shows  bias. 


Chapter  5  Results  and  Conclusions 
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Appendix  A  Input  Tima  Series,  GARSHOR.DAT 


Table  A1 


(T ens  of  step  numbers  in  left  column,  units  for  steps  across 
columns,  time  series  values  in  body  of  table) 


0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

0 

0.07 

0.12 

0.14 

0.15 

0.14 

0.09 

0.09 

0.08 

0.10 

10 

0.15 

0.24 

0.31 

0.33 

0.31 

0.23 

0.16 

0.09 

0.10 

0.10 

20 

0.08 

0.07 

0.04 

0.05 

0.05 

0.07 

0.09 

0.08 

0.08 

0.09 

30 

0.02 

0.01 

0.01 

0.03 

0.08 

0.18 

0.29 

0.31 

0.27 

0.18 

40 

0.11 

0.08 

0.09 

0.06 

0.02 

0.04 

0.23 

0.40 

0.52 

0.53 

SO 

-0.43 

0.37 

0.28 

0.14 

0.12 

0.11 

0.07 

0.00 

0.01 

0.01 

60 

-0.03 

0.00 

0.06 

0.08 

0.08 

0.07 

0.07 

0.05 

0.02 

0.07 

70 

-0.06 

0.01 

0.05 

0.10 

0.14 

0.16 

0.17 

0.14 

0.12 

0.12 

80 

0.16 

0.22 

0.24 

0.23 

0.22 

0.24 

0.26 

0.26 

0.25 

0.22 

so 

0.19 

0.17 

0.16 

0.17 

0.16 

0.17 

0.15 

0.12 

0.10 

0.08 

100 

0.12 

0.16 

0.21 

0.28 

0.28 

0.23 

0.15 

0.13 

0.06 

0.04 

110 

0.03 

0.00 

0.01 

0.06 

0.08 

0.10 

0.10 

0.04 

0.01 

0.05 

120 

-0.10 

0.08 

0.10 

0.09 

0.08 

0.04 

0.03 

0.02 

0.00 

0.04 

130 

0.08 

0.11 

0.11 

0.06 

0.01 

0.04 

0.06 

0.01 

0.06 

0.10 

140 

0.07 

0.08 

0.08 

0.04 

0.02 

0.00 

0.07 

0.10 

0.11 

0.13 

ISO 

0.09 

0.08 

0.15 

0.20 

0.25 

0.28 

0.23 

0.12 

0.03 

0.02 

160 

-0.04 

0.02 

0.03 

0.05 

0.04 

0.04 

0.05 

0.03 

0.01 

0.04 

170 

-0.03 

0.01 

0.02 

0.02 

0.01 

0.01 

0.07 

0.16 

0.22 

0.26 

180 

0.23 

0.15 

0.07 

0.02 

0.05 

0.04 

0.02 

0.09 

0.10 

0.16 

190 

0.17 

0.18 

0.17 

0.18 

0.17 

0.20 

0.19 

0.17 

0.11 

0.12 

200 

0.12 

0.18 

0.33 

0.48 

0.59 

0.58 

0.45 

0.25 

0.12 

0.10 

210 

0.10 

0.10 

0.08 

0.11 

0.13 

0.18 

0.23 

0.22 

0.19 

0.16 

220 

0.09 

0.04 

0.03 

0.08 

0.11 

0.17 

0.22 

0.31 

0.35 

0.34 

230 

0.29 

0.19 

0.07 

0.01 

0.01 

0.03 

0.06 

0.10 

0.14 

0.20 

240 

0.23 

0.28 

0.27 

0.25 

0.21 

0.17 

0.12 

0.07 

0.07 

0.10 

250 

0.18 

0.24 

0.26 

0.26 

0.21 

0.09 

0.05 

0.10 

0.10 

0.03 

260 

-0.02 

0.04 

0.05 

0.03 

0.03 

0.03 

0.01 

0.01 

0.01 

0.07 

270 

-0.11 

0.13 

0.15 

0.15 

0.14 

0.12 

0.11 

0.08 

0.07 

0.10 

280 

-0.14 

0.19 

0.18 

0.15 

0.28 

0.24 

0.16 

0.04 

0.08 

0.10 

290 

0.12 

0.16 

0.19 

0.18 

0.12 

0.09 

0.13 

0.13 

0.12 

0.12 

300 

0.13 

0.19 

0.24 

0.26 

0.24 

0.21 

0.20 

0.16 

0.12 

0.09 

310 

0.06 

0.05 

0.05 

0.07 

0.05 

0.07 

0.08 

0.08 

0.03 

0.00 

320 

0.00 

0.04 

0.04 

0.03 

0.02 

0.02 

0.02 

0.07 

0.09 

0.08 

330 

0.05 

0.06 

0.11 

0.15 

0.17 

0.16 

0.17 

0.18 

0.18 

0.17 

340 

-0.11 

0.09 

0.11 

0.11 

0.11 

0.07 

0.08 

0.11 

0.15 

0.19 

350 

0.21 

0.23 

0.22 

0.17 

0.12 

0.09 

0.10 

0.13 

0.16 

0.13 

360 

0.11 

0.07 

0.06 

0,05 

0.05 

0.04 

0.04 

0.01 

0.09 

0.17 

370 

0.20 

0.21 

0.19 

0.14 

0.03 

0.04 

0.08 

0.05 

0.02 

0.07 

380 

0.07 

0.05 

0.02 

0.02 

0.05 

0.07 

0.09 

-0.10 

0.10 

0.09 

390 

0.08 

0.04 

0.01 

0.06 

0.11 

0.12 

0.11 

0.09 

0.03 

0.03 

(Continumtl 
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Table  A1  (Continued) 


0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

400 

-0.10 

-0.13 

-0.10 

-0.06 

-0.06 

-0.09 

•0.13 

-0.17 

-0.20 

-0.19 

410 

-0.17 

-0.17 

•0.21 

•0.28 

•0.28 

-0.23 

-0.16 

-0.07 

0.03 

0.11 

420 

0.10 

0.04 

•0.04 

-0.11 

-0.19 

-0.19 

-0.17 

-0.11 

-0.09 

-0.10 

430 

-0.14 

-0.16 

-0.17 

-0.19 

-0.17 

-0.13 

-0.11 

-0.14 

-0.18 

-0.17 

440 

-0.12 

•0.05 

-0.01 

0.02 

0.02 

0.01 

-0.04 

-0.10 

•0.16 

-0.19 

460 

•0.18 

-0.17 

-0.17 

-0.15 

-0.14 

-0.17 

-0.20 

-0.18 

-0.19 

-0.20 

460 

•0.18 

-0.18 

-0.18 

•0.10 

0.00 

0.10 

0.16 

0.22 

0.26 

0.24 

470 

0.20 

0.15 

0.06 

0.04 

0.02 

0.02 

0.09 

0.19 

0.24 

0.28 

480 

0.25 

0.19 

0.11 

0.07 

0.01 

0.06 

0.02 

0.04 

0.14 

0.25 

480 

0.38 

0.47 

0.48 

0.46 

0.37 

0.22 

0.13 

0.04 

0.04 

0.09 

500 

0.09 

0.15 

0.28 

0.42 

0.50 

0.44 

0.37 

0.26 

0.12 

0.03 

510 

•0.03 

0.00 

0.18 

0.34 

0.43 

0.52 

0.59 

0.58 

Q.SO 

0.38 

520 

0.22 

0.16 

0.16 

0.19 

0.23 

0.24 

0.28 

0.31 

0.31 

0.26 

630 

0.22 

0.22 

0.19 

0.17 

0.17 

0.21 

0.26 

0.33 

0.36 

0.37 

540 

0.34 

0.31 

0.23 

0.15 

0.11 

0.15 

0.16 

0.17 

0.16 

0.15 

SSO 

0.14 

0.16 

0.14 

0.11 

0.05 

0.02 

0.03 

0.03 

0.03 

0.09 

560 

0.14 

0.18 

0.20 

0.20 

0.16 

0.11 

0.04 

-0.01 

-0.04 

0.00 

670 

0.02 

0.04 

0.06 

0.04 

0.01 

0.01 

-0.01 

•0.03 

-0.07 

-0.08 

580 

•0.07 

-0.05 

-0.02 

0.04 

0.12 

0.17 

0.15 

0.15 

0.12 

0.08 

590 

0.01 

•0.06 

-0.08 

-0.05 

-0.01 

0.03 

0.05 

0.04 

0.02 

0.01 

600 

0.00 

-0.03 

•0.08 

-0.11 

-0.11 

-0.12 

•0.14 

-0.08 

0.00 

0.09 

610 

0.10 

0.08 

0.08 

0.06 

-0.01 

-0.06 

-0.11 

-0.10 

-0.03 

0.04 

620 

0.07 

0.08 

0.06 

0.04 

0.03 

•0.02 

-0.06 

-0.04 

0.00 

0.04 

630 

0.07 

0.11 

0.21 

0.31 

0.35 

0.37 

0.37 

0.34 

0.27 

0.20 

640 

0.12 

0.10 

0.14 

0.18 

0.17 

0.17 

0.13 

0.10 

0.08 

0.05 

650 

0.03 

0.01 

0.00 

0.00 

-0.01 

0.01 

0.07 

0.16 

0.19 

0.22 

660 

0.21 

0.15 

0.09 

0.01 

-0.04 

-0.05 

-0.04 

-0.01 

0.01 

0.03 

670 

0.04 

0.04 

0.01 

-0.02 

•0.03 

-0.05 

-0.08 

•0.10 

-0.12 

-0.10 

680 

•0.04 

0.04 

0.11 

0.14 

0.10 

0.03 

•0.06 

•0.13 

-0.17 

-0.18 

690 

•0.18 

•0.16 

-0.12 

-0.10 

-0.13 

-0.19 

•0.22 

-0.25 

•0.26 

-0.29 

700 

•0.30 

•0.28 

•0.27 

•0.24 

•0.16 

-0.06 

0.05 

0.09 

0.11 

0.09 

710 

0.06 

•0.01 

•0.08 

-0.10 

•0.08 

•0.02 

0.03 

0.04 

0.08 

0.07 

720 

0.04 

0.00 

•0.02 

•0.05 

-0.07 

-0.12 

•0.16 

-0.16 

-0.13 

•0.04 

730 

0.02 

0.06 

0.08 

0.06 

0.03 

•0.04 

-0.10 

-0.11 

-0.10 

-0.09 

740 

•0.09 

-0.06 

-0.02 

0.02 

0.01 

•0.03 

•0.08 

-0.12 

-0.15 

-0.18 

750 

•0.17 

-0.13 

-0.10 

-0.08 

-0.07 

•0.05 

•0.01 

0.01 

0.00 

-0.09 

760 

•0.18 

-0.22 

•0.22 

-0.22 

-0.23 

-0.21 

-0.14 

-0.11 

-0.11 

-0.12 

770 

•0.12 

-0.11 

-0.12 

-0.14 

-0.16 

-0.15 

•0.12 

•0.09 

•0.05 

•0.04 

780 

•0.03 

-0.02 

-0.05 

-0.11 

-0.17 

-0.19 

-0.20 

-0.18 

-0.18 

-0.19 

790 

•0.21 

-0.20 

-0.21 

•0.20 

-0.18 

-0.15 

-0.14 

-0.17 

-0.18 

-0.16 

800 

•0.10 

•0.05 

0.00 

0.00 

-0.04 

•0.03 

-0.01 

-0.01 

0.02 

0.03 

810 

0.06 

0.08 

0.11 

0.09 

0.09 

0.09 

0.12 

0.11 

0.08 

0.10 

820 

0.14 

0.17 

0.21 

0.25 

0.29 

0.32 

0.34 

0.31 

0.29 

0.22 

830 

0.15 

o.n 

0.11 

0.14 

0.14 

0.16 

0.16 

0.15 

0.14 

0.06 

840 

0.00 

-0.02 

0.01 

0.04 

0.09 

0.15 

0.19 

0.22 

0.27 

0.31 

850 

0.35 

0.30 

0.22 

0.11 

0.01 

-0.05 

-0.04 

0.01 

0.09 

0.14 

860 

0.17 

0.15 

0.12 

0.07 

-0.01 

-0.06 

•0.06 

-0.05 

-0.04 

0.00 

870 

0.07 

0.13 

0.18 

0.25 

0.29 

0.28 

0.18 

0.04 

-0.09 

-0.17 

880 

•0.19 

-0.13 

-0.05 

0.00 

0.02 

0.06 

0.07 

0.03 

-0.05 

-0.13 

890 

•0.15 

-0.14 

-0.15 

-0.14 

-0.10 

-0.02 

0.07 

0.15 

0.22 

0.27 

900 

0.23 

0.15 

0.04 

-0.08 

-0.15 

-0.12 

-0.03 

0.04 

0.04 

0.07 

910 

0.09 

0.09 

0.09 

0.07 

0.01 

-0.04 

-0.10 

-0.12 

-0.15 

-0.12 

920 

-0.06 

0.04 

0.15 

0.26 

0.32 

0.27 

0.13 

•0.02 

-0.10 

-0.11 

930 

-0.08 

-0.05 

-0.02 

-0.02 

0.03 

0.10 

0.12 

0.11 

0.10 

0.04 

(Continued) 
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0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

•40 

0.01 

-0.03 

-0.07 

■0.03 

0.01 

0.00 

■0.01 

0.05 

0.11 

0.16 

•SO 

0.08 

-0.05 

-0.14 

-0.14 

■0.13 

■0.09 

■0.06 

-0.02 

-0.02 

0.01 

•60 

0.01 

0.01 

0  04 

0.07 

0.03 

-0.03 

-0.06 

■0.06 

-0.04 

-0.03 

•70 

-0.05 

■0.03 

0.02 

0.04 

0.00 

•0.06 

■0.11 

-0.11 

-0.11 

•0.16 

•SO 

■0.17 

-0.13 

-0.14 

-0.16 

•0.18 

-0.16 

■0.12 

-0.09 

■0.05 

■0.04 

•SO 

■0.07 

-0.11 

-0.13 

-0.13 

■0.13 

■0.10 

■0.05 

•0.03 

■0.01 

■0.04 

1000 

-0.08 

-0.08 

-0.06 

-0.09 

-0.10 

-0.06 

-0.05 

■0.09 

■0.12 

•0.15 

1010 

-0.13 

■0.06 

0.01 

0.06 

0.07 

0.06 

0.05 

0.07 

0.05 

0.03 

1020 

0.02 

0.01 

■0.03 

-0.08 

■0.08 

■0.05 

■0.01 

■0.03 

■0.08 

■0.09 

1030 

■0.11 

■0.12 

-0.16 

■0.19 

-0.21 

-0.16 

■0.13 

-0.12 

-0.09 

-0.04 

1040 

0.00 

-0.01 

-0.03 

■0.07 

■0.13 

-0.16 

■0.21 

■0.25 

■0.25 

-0.17 

10S0 

■0.12 

■0.11 

-0.12 

-0.16 

•0.20 

■0.24 

■0.27 

-0.31 

-0.33 

■0.32 

1060 

■0.28 

-0.26 

■0.20 

-0.10 

-0.01 

0.03 

0.05 

0.01 

■0.04 

■0.11 

1070 

-0.19 

-0.25 

•0.24 

■0.17 

-0.10 

■0.06 

■0.06 

■0.09 

■0.14 

-0.20 

1080 

-0.23 

-0.26 

■0.27 

■0.27 

-0.28 

■0.27 

■0.22 

-0.15 

■0.08 

0.00 

1090 

0.05 

0.03 

0.00 

•0.04 

-0.11 

■0.19 

■0.25 

■0.22 

-0.13 

-0.07 

1100 

0.02 

0.11 

-0.02 

-0.09 

-0.12 

-0.12 

-0.16 

-0.18 

■0.16 

■0.10 

1110 

■0.03 

0.11 

0.28 

0.41 

0.48 

0.44 

0.37 

0.24 

0.15 

0.04 

1120 

-0.02 

-0.01 

0.07 

0.19 

0.32 

0.38 

0.37 

0.34 

0.27 

0.16 

1130 

0.08 

0.04 

0.05 

0.08 

0.14 

0.23 

0.37 

0.48 

0.54 

0.53 

1140 

0.47 

0.40 

0.31 

0.18 

0.12 

0.14 

0.21 

0.28 

0.34 

0.39 

1150 

0.41 

0.38 

0.28 

0.21 

0.14 

0.12 

0.14 

0.14 

0.15 

0.18 

1160 

0.25 

0.33 

0.40 

0.39 

0.36 

0.31 

0.24 

0.17 

0.11 

0.13 

1170 

0.17 

0.21 

0.24 

0.26 

0.26 

0.25 

0.19 

0.13 

0.09 

0.11 

1180 

0.14 

0.16 

0.20 

0.21 

0.24 

0.30 

0.34 

0.35 

0.31 

0.29 

1190 

0.27 

0.20 

0.17 

0.16 

0.19 

0.22 

0.24 

0.27 

0.30 

0.28 

1200 

0.20 

0.13 

0.08 

0.10 

0.14 

0.15 

0.18 

0.19 

0.20 

0.23 

1210 

0.25 

0.23 

0.17 

0.09 

0.00 

-0.05 

■0.07 

-0.05 

■0.02 

0.03 

1220 

0.06 

0.08 

0.12 

0.12 

0.10 

0.01 

■0.04 

0.02 

0.02 

0.04 

1230 

0.09 

0.15 

0.15 

0.14 

0.17 

0.17 

0.11 

0.06 

0.02 

0.00 

1240 

■0.01 

-0.03 

-0.01 

0.00 

0.00 

0.02 

-0.02 

■0.06 

■0.10 

■0.18 

1250 

-0.19 

-0.15 

-0.05 

•0.04 

•0.05 

■0.02 

0.00 

0.04 

0.06 

0.07 

1260 

0.02 

■0.01 

•0.04 

■0.11 

•0.09 

■0.12 

-0.09 

•0.10 

■0.18 

■0.22 

1270 

■0.23 

-0.26 

-0.35 

•0.43 

•0.42 

•0.29 

-0.08 

0.04 

0.12 

0.20 

1280 

0.22 

0.19 

0.13 

0.03 

-0.10 

■0.18 

■0.20 

■0.22 

■0.12 

■0.04 

1290 

0.07 

0.11 

0.02 

•0.09 

■0.19 

■0.22 

-0.30 

■0.42 

■0.59 

■0.67 

1300 

-0.56 

-0.37 

-0.13 

0.08 

0.15 

0.12 

0.07 

■0.01 

■0.09 

■0.17 

1310 

-0.20 

-0.23 

-0.19 

■0.08 

0.04 

0.17 

0.20 

0.19 

0.11 

0.00 

1320 

■0.15 

■0.29 

-0.36 

■0.34 

-0.24 

■0.11 

0.04 

0.23 

0.39 

0.40 

1330 

0.35 

0.26 

0.15 

0.04 

■0.05 

■0.09 

•0.05 

0.06 

0.19 

0.31 

1340 

0.41 

0.42 

0.33 

0.18 

0.02 

■0.14 

■0.22 

■0.24 

■0.21 

•0.11 

1350 

0.04 

0.20 

0.34 

0.38 

0.38 

0.31 

0.22 

0.04 

■0.10 

-0.22 

1360 

-0.24 

■0.18 

-0.04 

0.09 

0.19 

0.20 

0.11 

0.03 

■0.10 

•0.23 

1370 

-0.36 

■0.43 

-0.42 

-0.31 

■0.17 

0.01 

0.19 

0.32 

0.40 

0.37 

1380 

0.24 

0.05 

-0.09 

-0.20 

-0.22 

■0.14 

0.00 

0.14 

0.29 

0.39 

1390 

0.38 

0.29 

0.15 

•0.02 

-0.21 

■0.33 

■0.36 

-0.35 

■0.25 

-0.15 

1400 

0.03 

0.22 

0.32 

0.31 

0.21 

0.05 

-0.13 

-0.31 

■0.40 

-0.40 

1410 

-0.30 

-0.16 

0.05 

0.22 

0.28 

0.27 

0.17 

0.03 

-0.18 

-0.26 

1420 

-0.32 

-0.34 

-0.34 

■0.23 

-0.01 

0.23 

0.39 

0.48 

0.48 

0.33 

1430 

0.13 

-0.10 

-0.25 

-0.26 

-0.15 

-0.01 

0.14 

0.28 

0.45 

0.48 

1440 

0.45 

0.35 

0.17 

0.03 

-0.08 

•0.18 

-0.25 

-0.20 

-0.05 

0.16 

1450 

0.39 

0.58 

0.64 

0.59 

0.43 

0.21 

0.03 

■0.15 

-0.16 

-0.08 

1460 

0.09 

0.23 

0.40 

0.59 

0.68 

0.66 

0.52 

0.41 

0.25 

0.12 

1470 

0.02 

-0.03 

0.04 

0.14 

0.25 

0.41 

0.52 

0.58 

0.57 

0.36 

1480 

0.13 

-0.09 

-0.24 

-0.28 

-0.22 

-0.10 

0.05 

0.23 

0.40 

0.49 

1490 

0.47 

0.37 

0.25 

0.05 

-0.10 

-0.22 

-0.25 

■0.17 

-0.02 

0.13 

iContinumI) 
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Table  A1  (Continued) 


0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

1600 

0.31 

0.44 

0.45 

0.35 

0.17 

-0.03 

■0.20 

-0.31 

-0.33 

-0.26 

1610 

■0.17 

-0.03 

0.14 

0.30 

0.40 

0.46 

0.42 

0.31 

0.13 

-0.04 

1520 

-0.16 

-0.20 

-0.16 

-0.08 

0.07 

0.26 

0.38 

0.40 

0.34 

0.21 

1530 

0.03 

-0.17 

-0.31 

-0.37 

-0.37 

-0.25 

-0.06 

0.12 

0.31 

0.44 

1540 

0.50 

0.46 

0.32 

0.15 

0.01 

•0.08 

-0.14 

-0.19 

-0.16 

•0.02 

1550 

0.16 

0.30 

0.39 

0.40 

0.32 

0.13 

-0.11 

-0.33 

-0.45 

-0.47 

1560 

-0.37 

-0.19 

0.01 

0.20 

0.36 

0.43 

0.35 

0.23 

0.13 

0.01 

1570 

-0.10 

-0.17 

-0.17 

•0.09 

0.06 

0.24 

0.40 

0.51 

0.53 

0.42 

1530 

0.25 

0.07 

-0.04 

•0.14 

•0.19 

-0.09 

0.09 

0.26 

0.47 

0.60 

1590 

0.62 

0.56 

0.49 

0.43 

0.32 

0.26 

0.13 

0.08 

0.12 

0.20 

1600 

0.29 

0.44 

0.54 

0.53 

0.49 

0.37 

0.21 

0.08 

-0.04 

-0.05 

1610 

-0.04 

0.07 

0.26 

0.39 

0.51 

0.61 

0.63 

0.60 

0.48 

0.34 

1620 

0.19 

0.10 

0.13 

0.17 

0.26 

0.42 

0.56 

0.64 

0.68 

0.63 

1630 

0.49 

0.36 

0.23 

0.14 

0.11 

0.17 

0.26 

0.40 

0.54 

0.64 

1640 

0.70 

0.70 

0.63 

0.52 

0.40 

0.28 

0.17 

0.10 

0.13 

0.28 

1650 

0.50 

0.69 

0.77 

0.72 

0.61 

0.45 

0.30 

0.17 

0.09 

0.11 

1660 

0.18 

0.28 

0.42 

0.56 

0.67 

0.70 

0.67 

0.55 

0.39 

0.26 

1670 

0.13 

0.03 

0.00 

0.13 

0.34 

0.54 

0.70 

0.71 

0.63 

0.47 

1680 

0.27 

0.07 

-0.09 

-0.13 

•0.08 

0.02 

0.15 

0.32 

0.50 

0.60 

1690 

0.57 

0.46 

0.29 

0.14 

0.01 

-0.13 

-0.17 

•0.08 

0.10 

0.30 

1700 

0.47 

0.55 

0.53 

0.44 

0.28 

0.09 

-0.07 

-0.10 

-0.04 

-0.02 

1710 

0.11 

0.29 

0.49 

0.59 

0.61 

0.58 

0.48 

0.36 

0.22 

0.12 

1720 

0.05 

0.05 

0.11 

0.20 

0.36 

0.48 

0.57 

0.55 

0.44 

0.23 

1730 

0.03 

-0.12 

-0.15 

-0.13 

-0.01 

0.11 

0.28 

0.40 

0.49 

0.46 

1740 

0.43 

0.32 

0.14 

-0.01 

-0.12 

-0.14 

-0.07 

0.01 

0.14 

0.29 

1750 

0.40 

0.43 

0.39 

0.25 

0.06 

-0.10 

•0.22 

-0.28 

-0.25 

-0.17 

1760 

-0.03 

0.10 

0.23 

0.26 

0.24 

0.17 

0.05 

-0.08 

-0.20 

-0.25 

1770 

-0.27 

-0.23 

•0.16 

-0.07 

0.05 

0.16 

0.18 

0.15 

0.04 

•0.10 

1780 

-0.22 

-0.33 

-0.35 

•0.31 

-0.17 

•0.02 

0.10 

0.21 

0.25 

0.23 

1790 

0.16 

0.03 

-0.10 

•0.18 

-0.21 

-0.19 

•0.12 

•0.03 

0.07 

0.20 

1800 

0.34 

0.37 

0.34 

0.23 

0.09 

-0.02 

•0.06 

-0.07 

0.00 

0.10 

1810 

0.21 

0.30 

0.37 

0.41 

0.37 

0.27 

0.14 

0.04 

-0.04 

-0.10 

1820 

-0.11 

-0.10 

-0.03 

0.05 

0.16 

0.24 

0.27 

0.24 

0.14 

0.02 

1830 

-0.08 

-0.13 

-0.15 

-0.12 

•0.04 

0.04 

0.13 

0.23 

0.24 

0.23 

1840 

0.16 

0.05 

-0.05 

-0,12 

•0.18 

-0.20 

•0.18 

-0.07 

0.05 

0.12 

1850 

0.17 

0.22 

0.20 

0.10 

-0.02 

-0.10 

-0.14 

-0.18 

-0.14 

-0.01 

1860 

0.09 

0.21 

0.33 

0.37 

0.37 

0.28 

0.17 

0.12 

0.07 

0.06 

1870 

0.07 

0.08 

0.15 

0.25 

0.34 

0.40 

0.45 

0.43 

0.36 

0.26 

1880 

0.21 

0.16 

0.11 

0.11 

0.13 

0.21 

0.28 

0.33 

0.32 

0.28 

1890 

0.24 

0.20 

0.10 

0.04 

-0.03 

-0.07 

-0.06 

-0.02 

0.08 

0.19 

1900 

0.31 

0.36 

0.31 

0.18 

0.04 

-0.05 

-0.12 

-0.17 

-0.20 

-0.11 

1910 

-0.04 

0.03 

0.10 

0.16 

0.13 

0.06 

-0.06 

-0.17 

•0.25 

-0.31 

1920 

-0.32 

-0.29 

-0.14 

0.03 

0.23 

0.37 

0.35 

0.27 

0,21 

0.15 

1930 

0.11 

0.04 

0.05 

0.15 

0.24 

0.28 

0.32 

0.37 

0.40 

0.38 

1940 

0.35 

0.30 

0.24 

0.20 

0.21 

0.22 

0.24 

0.30 

0.40 

0.50 

1950 

0.58 

0.58 

0.55 

0.56 

0.49 

0.44 

0.35 

0.34 

0.33 

0.36 

1960 

0.42 

0.47 

0.53 

0.51 

0.45 

0.45 

0.31 

0.18 

0.09 

0.03 

1970 

0.01 

0.05 

0.20 

0.34 

0.44 

0.47 

0.39 

0.32 

0.19 

0.08 

1980 

-0.02 

-0.03 

•0.02 

0.01 

0.06 

0.17 

0.27 

0.32 

0.27 

0.16 

1990 

0.02 

-0.07 

•0.18 

-0.28 

-0.32 

-0.22 

-0.03 

0.14 

0.25 

0.31 

2000 

0.33 

0.26 

0.16 

0.06 

•0.03 

-0.08 

-0.10 

-0.06 

0.04 

0.17 

2010 

0.30 

0.37 

0.34 

0.23 

0.11 

-0.01 

-0.14 

-0.23 

-0.27 

-0.23 

2020 

-0.13 

0.03 

0.18 

0.28 

0.35 

0.34 

0.24 

0.11 

-0.03 

-0.12 
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Table  A1  (Concluded) 


2030 

2040 

2060 

2060 

2070 

2080 

2090 

2100 

2110 

2120 

2130 


0 

1 

2 

3 

4 

6 

6 

7 

8 

9 

-0.14 

-0.11 

-0.03 

0.12 

0.27 

0.39 

0.41 

0.31 

0.16 

0.00 

-0.16 

-0.28 

-0.36 

•0.35 

-0.25 

•0.10 

0.08 

0.21 

0.31 

0.32 

0.26 

0.13 

•0.04 

•0.16 

•0.24 

-0.25 

-0.20 

0.05 

0.14 

0.29 

0.39 

0.34 

0.23 

0.06 

•0.11 

-0.25 

-0.38 

0.43 

0.39 

0.25 

-0.09 

0.06 

0.19 

0.30 

0.28 

0.18 

-0.02 

0.19 

0.31 

0.34 

-0.31 

-0.20 

-0.02 

0.17 

0.31 

0.36 

0.30 

0.19 

0.05 

0.13 

-0.28 

•0.41 

■0.46 

-0.37 

•0.24 

-0.09 

0.10 

0.26 

0.29 

0.19 

0.06 

-0.13 

•0.28 

-0.38 

-0.41 

•0.36 

0.17 

0.04 

0.21 

0.34 

0.38 

0.36 

0.26 

0.09 

•0.13 

-0.29 

0.39 

0.42 

0.38 

0.27 

-0.03 

0.20 

0.31 

0.34 

0.26 

0.09 

0.12 

0.29 

0.42 

0.49 

-0.38 

•0.16 

0.11 

0.23 

0.32 

0.34 

0.32 
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ooooooooooooooooooooooooooooooooooooooooooooooooooooooo 


PROGRAM  PREDICre 


PREDICT  (WITH  ERROR  ESTIMATE)  JO  *  M.JO  M  1 , JO  -t-  M  -t-  2 . JO  M  NPF-1 

TIME-STEP-VALUES  BASED  ON  THE  VALUES  AT  TIME  STEPS 

(JO,JO-«-1,JO-f2.JO-»^3 . J04^M-1) 

AND  FLAG  ANY  SERIES  VALUES  MORE  THAN  K*SIGMA 
FROM  THEIR  PREDICTORS  (NORMALY.  K-3) 

THE  TIME  SERIES  IS  INITIALLY  STORED  IN  FILE:  TIMSER.DAT 
AS  {VU);J-1,NJ},  WHERE  JO >-JB>-1  AND  (JO-t^M-«^NPF-1)< -NJ 

MAKE  COMPLETE  SEQUENCE  OF  ESTIMATES  FOR  JB<  -JO<  -NJ-M-NPF4 1 
(NOTE:  THE  NUMBER  OF  TIME  STEP  SEQUENCES  PREDICTED  IS 
NPTSoNJ-JB-M-f  1).  THE  PREDICTIONS  ARE  STORED  IN 
(VPK(JJ,LU^U-1,NPTS:a-1,NPF}  (FOR  KRIGING  ESTIMATES) 

AND  {VPA(JJ,LU^JJ-1.NPTS  LL-1.NPF}  (FOR  AUTOREGRESSIVE 
ESTIMATES).  THE  CORRESPONDING  ESTIMATES  OF  R.M.S.  ERROR 
ARE  STORED  IN  {ERK(JJ,LL)^)Jx1.NPTS;a«1,NPF)  (FOR  KRIGING) 

AND  {ERA(JJ,LU;JJ>1.NPTS:LL-1,NPF)  (FOR  AUTOREGRESSIVE 
RESULTS).  THE  ASSOCIATED  TIME  STEP  INDICES  ARE  STORED  IN 
{JTS(JJ.a):JJ  - 1  ,NPTS:LL  =  1  ,NPF) 

(VPK(JJ.LL);JJ-1,NPTS:a«1.NPF}  <•  KRIGING  PREDICTIONS, 

USING  DATA  FROM  TIME  STEPS  JJ-M  TO  JJ-1,  AND  MAKING 
PREDICTIONS  FOR  TIME  STEPS  JJ  TO  JJf  NPF-1 
(ERK(JJ.a):JJ-1,NPTS:a-1.NPF}  -  KRIGING  R.M.S.  ERROR 
AS  BASED  ON  ANO  FOR  TIME  STEPS  USTEO  ABOVE 

{VPA(JJ,U):JJ-1,NPTS:LL-1.NPF}  -  AUTOREG.  PREDICTIONS, 

USING  DATA  FROM  TIME  STEPS  JJ-M  TO  JJ-1,  AND  MAKING 
PREDICTIONS  FOR  TIME  STEPS  JJ  TO  JJ  +  NPF-1 
{ERA(JJ,LL):JJ-1,NPTS:U°1,NPF}  -  AUTOREG.  R.M.S.  ERROR 
AS  BASED  ON  ANO  FOR  TIME  STEPS  LISTED  ABOVE 

{JTS(JJ,LL);JJ«1,NPTS,-LL«1,NPF}  «  TIME  STEPS  ASSOCIATED 
WITH  VPK  AND  VPA  ARRAYS  ABOVE 

IF  (IFLAG.Ea.1)  PREPARE  ONLY  KRIGEO  ESTIMATES 
IF  (IFLAG.EQ.2)  COMPUTE  ONLY  AUTOREGRESSIVE  ESTIMATES 
IF  (IFLAG.EQ.3)  PREPARE  BOTH  KRIGED  AND  AUTOREGRESSIVE 
ESTIMATES 

NJ  «  LENGTH  OF  TIME  SERIES  DATA.  (IN  TIME  STEPS) 

M  «  LENGTH  OF  PREVIOUS  HISTORY  USED  IN  PREDICTION  OF  NEXT 
VALUE,  (IN  TIME  STEPS) 

MP1  -  M-t^l 

M05  «  AN  INTEGER  LARGER  THAN  (MP1/5) 

JB-t-M  >  TIME  STEP  NUMBER  FOR  THE  FIRST  PREDICTION 
NPTS  «  NUMBER  OF  PREDICTED  TIME  STEPS 
JO  »  STARTING  TIME  STEP  IN  MAKING  A  CURRENT  PREDICTION, 

(JO  GREATER  THAN  OR  EQUAL  TO  JB  AND  JO  LESS  THAN  OR 
EQUAL  TO  NJ-M.) 

NPF  -  NUMBER  OF  TIME  STEPS  INTO  THE  FUTURE,  STARTING  AT  JO  +  M, 

FOR  WHICH  PREDICTIONS  ARE  SIMULTANEOUSLY  MADE 
-  DURATION 
THOLD  -  THRESHHOLD 

C  (BYPASS  -=  0,  DONT  BYPASS  SECTION 
C  IBYPASS  -  1,  BYPASS  SECTION 

C  MPF  -  M-t-NPF 

C  FRACT  «  CUMULATIVE  FRACTION  OF  TRACE  WHICH  IS  REGARDED  AS 
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C  NEGUGIBlf  RELATIVE  TO  THE  SUM  OP  SMAUER 

C  EIGENVALUES.  THAT  IS,  THE  EIGENVALUES  ARE  RANKED 

C  IN  ASCENDING  ORDER  OF  ABSOLUTE  VALUE,  AND  THE 

C  CUMULATIVE  SUM  OF  ABSOLUTE  VALUES,  FROM  SMALLEST  TO 

C  LARGEST  IS  COMPUTED.  THE  SET  OF  EIGENVALUES  WITH 

C  CUMULATIVE  ABSOLUTE  VALUE  LESS  THAN  FR  AC 'TRACE  IS 

C  IGNORED  IN  COMPUTING  THE  GENERAUZED  INVERSES. 

C  (JO  M.JO  M  NPF-1 )  .  INTERVAL  OF  TIME  STEPS  BEING  PREDICTED  AT 
C  STARTING  STEPM,  JO 

C  LPREV  -  NUMBER  OF  TIME  STEPS  PRIOR  TO  JO  TO  BE  USED  IN 
C  CORRELATION  ESTIMATES 

C  M>LP  -  Mt-LPREV  (MUST  BE  A  PRODUCT  OF  POWERS  OF  2,  3,  AND  5|. 

C  DT  -  TIME  INCREMENT 

C  FLPASS  «  FREQUENCY  (HERTZ)  AT  WHICH  LOWER  FREQUENCIES  ARE 
C  PASSED  WITH  DISTORTION  LESS  IXAN  FFRAC. 

C  WIDTH  >  WIDTH  OF  TRANSITION  ZONE  FOR  FILTER,  IN  HERTZ,  FROM 
C  FLPASS  TO  A  HIGHER  FREQUENCY  WHERE  THE  SIGNAL  IS 

C  EUMINATED  WITH  REMAINDER  LESS  THAT  FFRAC. 

C  FFRAC  «  FRACTION  PERMITTED  FOR  RLTER  DISTORTION. 

C  MVMAX  -  MP1  *  IMP)  t- 1  )/2 
C  NMAX2  -  MP1  *MP1 
C 

C  THE  CORRELATION  FUNCTION  IS  ESTIMATED  WITH  THE  VALUES  IN  THE 
C  TIME  STEP  INTERVAL  (JO-L,  JO-t-M-1).  WHERE  (JO-L)>  AND 

C  L>-1 

C 

C  JEXAM  -  J-VALUE  FOR  WHICH  INTERMEDIATE  OUTPUT  IS  PRINTED.  IF 
C  (IOOC.NE.O) 

C  IF  IDOC  -  0,  IT  DOESN'T  MATTER  WHAT  JEXAM  IS  SET  TO. 

C 

C  WRITTEN  BY  LEON  BORGMAN,  UNIVERSITY  OF  WYOMING 
C  . . . 


C  . 

C  FFT  COMMON  BLOCK . . . 

PARAMETER  (MAXN  «  3000,MAXN2  «  2048.MAXN3  -  21 87,MAXNS  -  625) 
PARAMETER  (MAXN02  « 1024.MAXN03  -729,MAXN05  >  1 25) 

C  . . . . . 

DIMENSION  W(MAXN),WR(MAXN2),WI(MAXN2).WRD(MAXN3),WID(MAXN3) 
DIMENSION  WRT(MAXNS),WIT(MAXN5) 

DIMENSION  IRB(MAXN02I,C0(MAXN02I.SI(MAXN02) 

DIMENSION  IRB0(MAXN03),C0D(MAXN03).SID(MAXN03) 

DIMENSION  IRBT(MAXN05),C0T(MAXN05).Srr(MAXN05) 

DIMENSION  C01  (MAXN05),SI1  (MAXN05),C02(MAXN05),SI2(MAXN05) 
DIMENSION  C03(MAXN05),SI3(MAXN0S) 

COMMON  /FFT/W,WR,WI,WRD,WID.WRT.Wrr,IRB,CO,SI,IRBD.COD,SID 
COMMON  /FFT/IR8T,COT,SIT,C01  ,SI1  ,C02,SI2,C03,SI3 
C  . . . 


C  . PARAMETERS  TO  BE  SET  BY  USER  '  *  ' 

PARAMETER  (NJ  «  21 36,M  *  1 2.JB  «  53,NPTS  -  2073,LPREV  >  52) 
PARAMETER  (MP1  « 13,MPLP-64,IFLAG  «3) 

PARAMETER  (MVMAX  «496,NMAX2  -961  .NPF -  lO.MPF  - 22,FRAC  -0.01 ) 
PARAMETER  (M05-7,TH0L0«0.25,IBYPASS- 1) 

PARAMETER  (DT - 1 .0.FLPASS  -0.04,WIDTH «0.02,FFRAC  -0.05) 

C  . DIMENSION  LIST  *  '  * 

DIMENSION  V(NJ),VV(NJ),C0RR(MPF),C0RMAT(MP1  .MP1 ) 

DIMENSION  BKMPI  ,MP1  ),C0(MP1  ),CK(MP1  ),CR(M) 

DIMENSION  VPK(NPTS.NPF),ERK(NPTS,NPF),MLAG(NPTS) 
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DIMENSION  VPA(NPTS,NPF).ERA(NPTS.NPF),JTS(NPTS,NPF) 
DIMENSION  A<MVMAX),EIVEC(NMAX2) 

DIMENSION  INOEXdOO) 

OINENSION  lOBIMOSI.iaTIMOSI.IOAIMPI  I.QAIMPI  I 
DIMENSION  TIMSER(MPU>).WRKVEC(MPLP> 


C  . OUTPUT  FILE  ••• 

OPEN  (2.FILE  -  'PREOICTS.OUT'.STATUS  -  'UNKNOWN') 
OPEN  (4,RLE  -  RESULTS.OUT', STATUS  « 'UNKNOWN') 

C  . INPUT  DATA  FILE 

OPEN  I3,FILE  -  'aarSHOR.OAT', STATUS  « 'OLD') 

C  . PROGRAM  CONSTANTS  •  *  * 

IDOC-0 

JEXAM-S 


WRITE  (4,160)  THOLO 
ISO  FORMAT  (2X, THRESHOLD-  ',F8.2.//) 


C  . READ  INPUT  DATA  *  *  * 

READ  (3,111)  (W(l),l-1,JB^M-2) 

1 1 1  FORMAT  (F6.2) 

C  PRINT*,'INPUT  DATA  READ' 

•  •  *  -  . . START  COMPUTATIONS  *  •  * 

. LOOP  OVER  PREDICTION  TIME  STEPS  •  •  * 

DO  1  JJ«1,NPTS 

IST0P1  -0 
ISTOP2-0 

JO-JB-fJJ-1 
JJJ-JO+M-I 
READ  (3,111)  VVUJJ) 

00  80LTS-1,MPLP 
TIMSER(LTS)  -  VV(JO-LPREV  +  LTS-1 ) 

WRKVEC(LTS) -0.0 
80  CONTINUE 

CALLLOPASS(MPLP,FLPASS,WIOTH,TIMSER,WRKVEC,DT,FFRAC) 


00  81  LTS-1,MPLP 
V(JO-LPREV  +  LTS-1 )  -TIMSERILTS) 
81  CONTINUE 

WRITE  (4,151)  JJJ,V(JJJ) 

151  FORMAT  (2X,'JJJ,V(JJJ)-  ',I5,F8.2) 
PRINT*, 'JJJ,V=  ',JJJ,V(JJJ) 


MAXLAG-M 

. IF  READ  VALUE  IS  ABOVE  THRESHHOLD,  ESTIMATE  *** 

.  NPF  -  DURATION  FUTURE  VALUES  BEYOND  STEP  *  *  * 

.  JO  +  M-1 
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IF  (VgO-^M-D.GE.THOLO)  THEN 
WRITE  (4.1S2)  V(JJJ) 

152  FORMAT  (2X. 'EXCEEDANCE  OF  THRESHOLD  FOUND.  VALUE  > 
A  F8.2) 

PRINT*, ‘EXCEEDANCE  OF  THRESHOLD  FOUND.  VALUE- 
•i  VUJJ) 

. CONKRiTE  MEANS  OF  SET  •** 

PRINT'.'COMPUTING  MEANS' 

SUMV-0.0 

SUMJ-0.0 

DO  3  J  -  JO-LPREV.JO  M-1 
SUMV-SUMV-«-V<J) 

SUMJ«SUMJ>J 
3  CONTINUE 

AVV  -  SUMV/(M + LPREV) 

AVJ  -SUMJ/IM  +  LPREV) 

. COMPUTE  UNEAR  TREND  *•* 

PRINT'/FITTING  UNEAR  TREND' 

SXX-0.0 

SYY-0.0 

SXY-0.0 

DO  5  J  -  JO-LPREV,JO  *  M-1 
SXX  -  SXX  g-AVJ)  •  IJ-AVJ) 

SYY  -  SYY  +  (VIJ)- AVV)  •  (Vg)-AVV) 

SXY  -  SXY  +  (J-AVJ)  •  (V|J)-A  VV) 

5  CONTINUE 
BH-SXY/SXX 
AH-AVV-BH*AVJ 

IF  (IOOC.NE.O.ANO  JJ.EQ.JEXAM)  THEN 
WRITE  (2.200)  JJ.AH.BH 

200  FORMAT  (/,2X, 'JJ.AH.BH-  ',I5.2F12.5,//) 

END  IF 


. COMPUTE  CORRELATION  FUNCTION  FOR  LAG  -  0 . M 

PRINT*. 'CORRELATION  FUNCTION' 

DO  2KK-1.MPF 
K-KK-1 

C  . LAG  K  COMPUTATIONS  *  *  * 

IF  (K.LE.MAXLAG)  THEN 
SUM -0.0 

DO  4  J  -  JO-LPRE  V  >  K,  JO -f  M- 1 
JF-J-K 

SUM  -SUM  +  (V(JF)-AH-BH*  JF)  ‘(VOI-AH-BH  *  J) 

4  CONTINUE 

CORR(KK)  -SUM/(M+  LPREV-K) 

IF  (KK.EQ.1)  THEN 
COVO-CORR(KK) 

CORR(KK)-1.0 
ELSE  IF  (KK.GT.1)  THEN 
CORRIKK)  -  CORR(KK)/COVO 
END  IF 

IF  ICORR(KK).LE.O.O)  THEN 
MAXLAG-K 
C0RR(KK)-0.0 
END  IF 
ELSE 

C0RR|KK)-0.0 

END  IF  _ 
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IF  (IOOC.NE.O.ANO  JJ.EQ  JEXAM)  THEN 
WRITE  12.205)  COVO 

205  FORMAT  (/.2X.‘C0V0- '.FI  2.5) 

WRITE  (2.201)  JJ.KK.CORR(KK) 

201  FORMAT  (2X.'JJ.KK.C0RR-  '.2IS.F12.5) 

END  IF 

2  CONTINUE 

C  . LOAD  UPPER  LEFT  MAXX  BY  MAXX  PART  OF  CORMAT  •  •  * 

C  PRINT*.'SET  UP  CORRELATION  MATRIX' 


IF  (MAXLAG.GT.2)  THEN 
MAXX-MAXLAG 
ELSE 

MAXX-2 
END  IF 

MLAGfJJI-MAXX 
D0  6L-1.MAXX 


00  7  a*1.MAXX 
LAG-IABSIL-a) 

K«LAG  +  1 

CORMAT(L.LL)  -  CORR(K) 

7  CONTINUE 

6  CONTINUE 

C  . IF  IFLAG  - 1  OR  3.  DO  POINT  KRIGING  *  •  * 

IF  (IFLAG .Ea.1.0R.IFLAG.EQ.3)  7XEN 
C  PRINT'.'PROCEEO  TO  MAKE  KRIGED  ESTIMATE' 

C  . .  IN  THE  RIGHT  AND  BOTTOM  MARGIN  OF  CORMAT  *  *  • 

C  .  AND  THE  RIGHT-HAND  VECTOR  ••• 

C  PRINT*, 'PREPARE  KRIGING  MATRIX' 


DO  8  J  -  JO + M-MAXX.JO  M-1 
L-J-JO-M-fMAXX-f  1 
MAXXP-MAXX-t-1 
CORMATIUMAXXP)  •  AH  4  BH*  J 
CORMATIMAXXP.U  -  CORMATIUMAXXP) 

8  CONTINUE 

CORMATIMAXXP.MAXXP)  -0.0 

IF  (IOOC.NE.O.ANO  JJ.EQ.JEXAM)  THEN 
WRITE  (2.202) 

202  FORMAT  (//,2X. 'KRIGING  MATRIX:') 

CAa  WRMATICORMAT.MAXXP.MAXXP.MPI  ,MP1  .INDEX) 
WRITE  (2.'(//)'J 
END  IF 


C  . COMPUTE  THE  GENERALIZED  INVERSE  OF  AUGMENTED  *  *  * 

C  .  CORMAT 

CALL  EGINV2(C0RMAT,BI.FRAC.MAXXP,MP1  ,MVMAX,NMAX2, 

&  EIVEC.A,QA,IQA.IQB,IQT,M05) 

C  . COMPUTE  KRIGING  PREDICTIONS  FOR  NPF  FUTURE  STEPS  *  *  * 

C  PRINT*, 'COMPUTING  KRIGING  COEFFICIENTS' 

C  . SET  FLAG  FOR  STOPPING  *  *  * 

IST0P1  -1 
D0  30  JPF-1,NPF 

C  . CREATE  CO  VECTOR  *  *  * 

DO  31  L-1,MAXX 
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CO<U  «CORR(MAXX-L  *1+JPF) 

31  CONTINUE 

A  VP  -  AH  BH  *  (JO  M -f  JPF- 11 
CO(MAXXP)-AVP 

IF  (IDOC.NE.O.ANDJJ.EQJEXAM)  THEN 


WRITE  (2,134) 

134  FORMAT  U.2X/C0  VECTOR:*) 

00  32L-1.MAXXP 
WRITE  (2,206)  1„C0(U 

206  FORMAT  (2X,'L,C0- ‘,IS,F12.S) 

32  CONTINUE 

END  IF 

C  . MULTIPLY  BI*CO  TO  GET  KRI6ING  COEF.  *  •  • 

IF  (IDOC.NE.O.ANOJJ.EQJEXAM)  THEN 
WRITE  (2,136) 

136  FORMAT  (/,2X,'CK  VECTOR:*) 


END  IF 

00  9L-1,MAXXP 
SUM -0.0 

00  10  LL-1,MAXXP 
SUM  -SUM  BI(L.LL)  “COIU) 

10  CONTINUE 
CK(L)-SUM 

IF  (IDOC.NE.O.ANO.JJ.EQ.JEXAM)  THEN 
WRITE  (2,204)  L,CK(L) 

204  FORMAT  (2X,'L,CK-  *,I5.F12.S) 

END  IF 

9  CONTINUE 

. COMPUTE  ESTIMATE,  VPK,  AND  ERROR,  ERK  ••• 

PRINT*,*CALCUUTING  KRIGEO  ESTIMATE* 

SUM -0.0 
SUMM-0.0 
00  11  L-1,MAXX 
SUM  -SUM  -I-  CK(L)  *CO(LI  'COVO 
J-JO-i-M-MAXX-t-L-1 
SUMM  -  SUMM  >  CK(L)  *  V(J) 

11  CONTINUE 

ERRR  -  COVO-CKIMAXXP)  'COVO  *  AVP-SUM 
ERK(JJ.JPF)  -SQRT(ABS(ERRR)) 

JTS(JJ,JPF)  -  JO  +  M  +  JPF-1 
VPK(JJ.JPF)  -  SUMM 


C  . IF  VPK  LESS  THAN  THOLD,  SET  IST0P1  -0  •  •  • 

IF  (SUMM.LT.THOLO)  IST0P1  -0 

30  CONTINUE 

C  . STOP  IF  ALL  VPK.GE.THOLD 

IF  (IST0P1.EQ.1)THEN 
C  STOP 

WRITE  (4,153)  JJJ 

153  FORMAT  (2X,’KRIGING  SHUTDOWN  STEP  AT  *,16) 

PRINT*, *KRIQING  SHUTDOWN  STEP  AT  *,JJJ 


END  IF 
END  IF 


App«ndix  B  Pradietion  Computer  Program 


B7 


ooo  oooooo 


c 


& 

33 

137 

13 

208 

12 

209 


51 

so 

212 


. IF  IFLAG  -2  OR  3.  DO  AUTOREGRESSIVE  ESTIMATE  *  *  * 

IF  (IFLAG.Ea.2.0R.IFLAG.EQ.3)  THEN 

PRINT  V'PROCEEO  WITH  AUTOREGRESSIVE  ESTIMATES' 

. COMPUTE  THE  GENERAUZEO  INVERSE  OF  •  *  * 

.  NON-AUGMENTEO  CORMAT 

PRINT*,‘INVERT  MATRIX  IN  YULE-WALKER  EQUATIONS' 
CAU  EGINV2(C0RMAT.BI,FRAC.MAXX,MP1  ,MVMAX.NMAX2, 
EIVEC,A,QA.IQA,IQB,I0T,M05) 

. SET  UP  RIGHT  SIDE  OF  YULE-WALKER  EQUATION  *  •  • 

00  33  L-1.MAXX 
COIL) «  CORRIMAXX-L  4  2) 

CONTINUE 

. MULTIPLY  Bi'CO  •** 

PRINT*, 'COMPUTING  AUTOREGRESSIVE  COEFRCIENTS' 

IF  (IOOC.NE.O.ANO  JJ-EO  JEXAM)  THEN 
WRITE  (2,137) 

FORMAT  (/,2X,'CR  VECTOR:') 

END  IF 

SUMCR-0.0 
00  12  L->1,MAXX 
SUM -0.0 
DO  13  a-1,MAXX 
SUM  -  SUM  4  BIIUU)  ‘COILU 
CONTINUE 
CRIU-SUM 

SUMCR-SUMCR4CR(L) 

IF  (IDOC.NE.O.ANO.JJ.EQJEXAM)  THEN 
WRITE  (2,208)  L,CR(L) 

FORMAT  (2X,'L,CR-  ',(5,Ft2.6) 

END  IF 

CONTINUE 

IF  (IOOC.NE.O.ANO.JJ.EQ.JEXAM)  THEN 
WRITE  (2,209)  SUMCR 
FORMAT  (//,2X,'SUM  AR  COEF-  ’,F12.5,//> 

END  IF 


. ESI  IMATE  VALUE  AND  ERROR  FOR  AR  EST.  •  •  * 

. COMPUTE  ESTIMATE,  VPA,  AND  ERROR,  ERA  **• 

PRINT*, 'CALCULATING  AUTOREGRESSIVE  ESTIMATES' 

. CHANGE  CORRELATION  MATRIX  TO  COVARIANCE  MATRIX  *** 

00  50L-1,MAXX 
DO  51  LL-1,MAXX 
CORMAT(L,a)  -C0RMAT(L,U)*C0V0 
CONTINUE 
CONTINUE 

IF  (IOOC.NE.O.ANO.JJ.EQ.JEXAM)  THEN 
WRITE  (2,212) 

FORMAT  (//,2X,'INrriAL  COV  MATRIX:') 

CAU  WRMAT(C0RMAT,MAXX,MAXX,MP1  ,MP1  ,INDEX) 

WRITE  (2, '(//)') 

END  IF 
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ISTOP2-1 
D0  40  JPF-1.NPF 
SUMV-0.0 
DO  41  L-2,MAXX 
CCC  -  CR(1)  *CORMAT(L.1) 

CORMAT(L,MAXXP)  >CCC 
CORMAT(MAXXP,U  -  2.0  *  CCC 
SUM  >0.0 
DO  42a-2.MAXX 
SUM  -  SUM  *  CORMAT(L.ai  ‘CRtU) 

42  CONTINUE 

CORMATIL,MAXXP)  >  CORMAT(L.MAXXP)  SUM 
CORMAT(MAXXP.U  >CORMAT(MAXXP.U  +  SUM 
41  CONTINUE 

SUM  >0.0 
00  43L>2.MAXX 

SUM  >  SUM  -f  CR(L)  *CORMAT(MAXXP,L) 

43  CONTINUE 
CORMATIMAXXP.MAXXP)  >SUM 
00  44L>2.MAXX 

00  4Sa-2,MAXX 
CORMATIL- 1  ,LL- 1 1  -  CORMATIL.LU 

45  CONTINUE 

44  CONTINUE 

00  46  L>2,MAXX 

CORMATIL- 1  ,MAXXP-1 )  >  CORMATIL.MAXXP) 
CORMATIMAXXP- 1 X- 1 1  -  CORMAT<L,MAXXP| 

46  CONTINUE 

CORMAT(MAXX,MAXX)>CORMAT<MAXXP,MAXXP) 

IF  (IOOC.NE.O.ANO.JJ.EQJEXAM)  THEN 
WRITE  (2,213) 

213  FORMAT  <//,2X,'NEXT  COV  MATRIX;') 

CALL  WRMAT(C0RMAT,MAXX,MAXX,MP1  ,MP1  .INDEX) 
WRITE  |2,'(/n') 

END  IF 

AVP»AH  +  BH*(J0  +  M+JPF-1) 

00  47  L>1,MAXX 
JV  >  JO  -f  M  JPF-MAXX  *  L-2 
AVV>AH  +  BH»JV 
IF  IJV.LT.JO  +  M)  THEN 
WW>V(JV)-AVV 
ELSE 

WW  «  VPAIJJ,  JPF-MAXX  +  L- 1  )-A  VV 
END  IF 

SUMV  >  SUMV  +  CfllL)  •  WW 

47  CONTINUE 

VPAIJJ, JPF)  «SUMV  +  AVP 
JTSIJJ,  JPF)  -  JO  +  M  +  JPF- 1 
ERR  >  CORMATIMAXX.MAXX) 

ERA|JJ,JPF)  -SQRT(ABSIERR)) 


C  . IF  VPA  LESS  THAN  THOLD.  SET  IST0P2  >0  •  *  • 

IF  (VPAIJJ, JPR.LT.THOLD)  IST0P2-0 
40  CONTINUE 

C  . STOP  IF  AU  VPA.GE.THOLD  ”  •  • 

IF  IIST0P2.EQ.1)  THEN 
C  STOP 

WRITE  14, 1 54)  JJJ  _ 
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1 64  FORMAT  (2X,‘AR  SHUTDOWN  STEP  AT  MSI 

PRINT*, 'AR  SHUTDOWN  STEP  AT  ‘JJJ 
END  IF 

END  IF 
END  IF 

1  CONTINUE 

IF  (IBYPASS.Ea.OI  THEN 

•  •  •  *  *  PRINT  PREDICTIONS  *  •  * 

PRINT*.‘PRINT  PREDICTIONS' 

DO  17JJ-1.NPTS 


WRITE  (2,130)  JJJ,MLAG{JJ) 

130  FORMAT  (//,1X,'LAST  TIME  STEP  USED  FOR  PREDICTION-  ',IS, 
&  SX,'MAXLAG-  ',15,/) 

C  . WRITE  TIME  SERIES  PRECEEDING  PREDICTIONS  •  *  * 

JENO-JB-i-M-fJJ-2 
JSTART- JENO-MLAG(JJI  -t- 1 
WRITE  (2,131) 

131  FORMAT  (/,2X,'0ATA  PRECEEDING  FIRST  PREDICTION:') 

DO  18  J-JSTART,JENO 

WRITE  (2,113)  J,V(J) 

1 1 3  FORMAT  (2X,I5,F1 0.3) 

18  CONTINUE 

WRITE  (2.'(1X)') 

C  . WRITE  PREDICTIONS  ••• 


I  WRITE  (2,120) 

120  FORMAT  (2X,'  J',8X,'V'.7X.'VPK'.7X,'ERK',7X,'VPA', 

&  7X,'ERA:') 

00  48  JPF-1,NPF 

WRITE  (2,112)  JTS(JJ,JPF),V(JTS(JJ,JPF)),VPK(JJ,JPF), 

<1  ERKUU.JPR,VPA(JJ,JPF),ERA(JJ.JPF) 

1 1 2  FORMAT  (2X,I5,SF10.3) 

48  CONTINUE 
17  CONTINUE 
END  IF 

STOP 

END 

C 

SUBROUTINE  EGINV2(B,BI,FRAC,N.NMAX,MVMAX,NMAX2,EIVEC,A, 
&QA,IQA,IQB.IQT,NMOS) 

C  . 

C  COMPUTE  GENERAUZEO  INVERSE  FOR  A  SQUARE,  SYMMETRIC 
C  MATRIX,  WITH  EIGENVECTOR  METHODS. 

C 

C  B  -  N  BY  N  INPUT  MATRIX 
C  Bl  -  INVERSE  OF  B 

C  FRAC  «  VALUE  OF  EIGENVALUE  CUTOFF.  NEGUGIBLE  EIGENVALUES 
C  SUMMING  TO  A  VALUE  LESS  THAN  TRACE  *  FRAC  ARE  IGNORED. 

C  (EIGENVALUES  ARE  ELIMINATED  IN  ASCENDING  ORDER  OF 

C  ABSOLUTE  VALUES.) 

C  NMAX  -  INTEGER  >  -  N,  FOR  DIMENSIONING 
C  MVMAX  -  NMAX*(NMAX-»^1)/2 
C  NMAX2  -  NMAX 'NMAX 


BIO 
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ElVEC  •  NMAX2-C0M>0NENT  WORK  VECTOR 
A  -  MVMAX-COMPONENT  WORK  VECTOR 

WRITTEN  BY  LEON  BORGMAN,  UNIVERSITY  OF  WYOMING 


DIMENSION  B(NMAX,NMAX).BI(NMAX,NMAX) 

DIMENSION  A(MVMAX).EIVEC(NMAX2) 

DIMENSION  IQB(NMOS).IQT(NMOS),IQA(NMAX),QA(NMAX) 

. SET  1C  -0  TO  SIGNAL  SUBR  EIG  TO  COMPUTE  BOTH  •  •  • 

.  EIGENVALUES  ANDS  EIGENVECTORS 

IC>0 

. CREATE  A  VECTOR  FOR  SUBR.  EIG  AND  COMPUTE  TRACE  *  *  * 

TRACE  «  0.0 
L-0 

DO  1  J-1.N 
D0  2U1,J 
L-L-l-1 
A(L)-B(I,J) 

IF  tI.EQ  J)  TRACE  -TRACE  -t-  B(UI 

2  CONTINUE 
1  CONTINUE 

CUT-TRACE‘FRAC 

. COMPUTE  EIGENVECTORS  AND  EIGENVALUES  **• 

CALLEIG(A,EIVEC,N.IC,MVMAX,NMAX2) 

. LOAD  ABSOLUTE  VALUE  OF  EIGENVALUES  FOR  SORT  ••• 

L-0 

00  8  J-1,N 
D0  9I-1,J 
L-L+1 

IF  (I.EQJ)  THEN 
aA(ll-ASS(A(L)) 

IQA(I)-I 
END  IF 

9  CONTINUE 
8  CONTINUE 

CAaFQS0RT(QA,IQA,N,NM05,IQB.IQT,NI) 

. ZEROBI*** 

00  3  l-1,N 
D04  J-1,N 
BI(l,JI-0.0 
4  CONTINUE 

3  CONTINUE 

. COMPUTE  VIL  INVERSE) ‘(V  TRANSPOSE),  STORE  IN  Bl  '•* 

ESUM-0.0 
00  5KK-1,N 

. COMPUTE  WHERE  THE  K-TH  EIGENVALUE  IS  IN  A  ••• 

K-IOAIKK) 

L»K»(K+1)/2 
ESUM  -  ESUM  *  QAIKK) 

IF  (ESUM.GT.CUT)  THEN 

. IF  K-TH  CUM.  SUM  EXCEEDS  FRAC'TRACE,  ADD  DYADIC  • 

.  INTOBI 

00  6I-1,N 

. COMPUTE  WHERE  THE  l-TH  COMPONENT  OF  THE  •  *  • 

.  K-TH  EIGENVECTOR  IS  IN  ElVEC 
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00  7  J-1,N 

. COMPUTE  WHERE  THE  J-TH  COMPONENT  OF  THE  * 

.  K-TH  EIGENVECTOR  IS  IN  ElVEC 

JK-IK-D’N-f  J 

C  . SUM  INTO  B1 ‘o* 

Bid,  J)  -  Bid  >  El  VECdK)  ‘EIVECUIO/AIU 
7  CONTINUE 
6  CONTINUE 
END  IF 
5  CONTINUE 

RETURN 

END 

C 

SUBROUTINE  EIG(A,EIVEC,N.IC,MVMAX.NMAX2) 

C  . 

C  CONFUTE  EIGENVALUES  AND  EIGENVECTORS  OF  A  REAL  SYMMETRIC 
C  MATRIX 

C  A  -  ORIGINAL  MATRIX  (OF  ORDER  N)  IS  DESTROYED  IN  COMPUTATION. 

C  THE  ELEMENTS  OF  A  ARE  ORDERED  INTO  A  VECTOR  OF  UPPER 
C  TRIANGULAR  ELEMENTS  AS  SHOWN  IN  THE  SEQUENCING  PATTERN 

C 

C  1  2  4  7  11  16  THE  EIGENVALUES  ARE  CONTAINED 

C  3  5  8  12  17  IN  THE  POSITIONS  OF  THE  DESTROYED 

C  6  9  13  18  ETC.  MATRIX  A  CORRESPONDING  TO  THE 

C  10  14  19  MAIN  DIAGONAL.  THAT  IS:  1.3.6,10. 

C  IS  20  15,21,  ETC. 

C  21 

C  THE  >TH  EIGENVALUE  IS  IN  POSITION  J«(J  4 1  )/2 

C 

C  ElVEC  •  MA  TRIX  OF  EIGENVECTORS  .STORED  IN  A  SINGLE  VECTOR. 

C  ONE  AFTER  THE  OTHER,  IN  SAME  SEQUENCE  AS  EIGENVALUES. 

C  THE  >TH  EIGENVECTOR  WOULD  OCCUPY  THE  SEQUENCE 

C  OF  COMPONENTS  OF  ElVEC,  STARTING  AT  COMPONENT 

C  U-1)*N  1  AND  ENDING  AT  COMPONENT  J*N 

C  N  -  ORDER  OF  MATRICES  A  AND  ElVEC 

C  NMAX  -  MAXIMUM  SIZE  OF  N.  FOR  DIMENSIONING  ONLY.  (NOTE:  NMAX 
C  IS  NOT  REQUIRED  AS  INPUT.  BUT  DETERMINES  MVMAX 

C  AND  NMAX2.I 

C  MV  -  THE  DIMENSION  OF  THE  VECTOR  A.  -N*IN *W2 
C  MVMAX  -  MAXIMUM  SIZE  OF  MV,  FOR  DIMENSIONING  ONLY 
C  «  NMAX* (NMAX 1)12 

C  NMAX2  -  THE  DIMENSION  OF  THE  VECTOR  ElVEC,  -  NMAX*  *2 
C  1C  -  0  COMPUTE  EIGENVALUES  AND  EIGENVECTORS 
C  1C  «  1  COMPUTE  EIGENVALUES  ONLY  (ElVEC  NEED  NOT  BE 
C  DIMENSIONED  BUT  MUST  STILL  APPEAR  IN  CALUNG 

C  SEQUENCE) 

C 

C  DIAGONAUZATION  METHOD  FROM  'MATHEMATICAL 

C  METHODS  FOR  DIGITAL  COMPUTERS',  EDITED  BY  A.  RALSTON  AND 

C  H.S.  WILF,  JOHN  WILEY  AND  SONS,  NEW  YORK.  1962,  CHAPTER  7 

C 

C  FOR  A  DOUBLE  PRECISION  VERSION,  REMOVE  THE  LEADING  'C' 

C  FRO'  THE  FOaOWING  STATEMENT: 

C  DOUBLu  ;  RECISION  A,EIVEC,AN0RM,ANRMX,THR,X,Y,SINX.StNX2.C0SX, 
C  &  C0SX2,SINCS.C0NST 

C  ALSO  CHANGE:  SORT  TO  OSQRT,  ABS  TO  DABS.  AND  CONST  IN 
C  FIRST  STATEMENT  BELOW  TO  1  .OD-1 2 

C  . * . 

DIMENSION  A(MVMAX),EIVEC(NMAX2) 
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CONST- 1.0E-6 
IF  (IC.NE.1)  THEN 
IQ— N 

00  20  J-1,N 
IQ-IQ-t-N 
OO20I-1.N 
U-IQ-fl 
EIVEC(U)-0.0 
IF  (I.EQ.J)  THEN 
EIVEC(U)-1.0 
END  IF 

20  CONTINUE 
END  IF 
ANOftM-0.0 
D0  3SI-1.N 
DO  35  J-I.N 
IF  (I.NE  J)  THEN 
IA-l-t-U*J-J)/2 

ANORM  «  ANORM  *  A(IA)  *  Ad  A) 

END  IF 

35  CONTINUE 
IF  (ANORM.GT.0.0)  THEN 
ANORM  - 1 .414*SORT(ANORMI 
ANRMX  -  ANORM*CONST/N 
INO-0 

THR- ANORM 
45  THR-THR/FUOATIN) 

50  L-1 

55  M-L-fl 
60  MQ-(M*M-M)/2 

LQ-(L*L-U/2 
LM-L-t-MQ 

IF  |ABS(A(LM)).6E.THR)  THEN 
IND-1 
U.-L4LQ 
MM-M4^MQ 
X-0.5MA(LU-A(MM)) 

Y  -  •  AILMI/  SQRTIAILM)  *  A(LM)  X  *  X) 

IF  (X.LT.0.01  THEN 
Y— Y 
END  IF 

SINX-Y/  SQRT«2.0*ll.0  +  (  SQRT(1.0-Y*Y)m 

SINX2-SINX*StNX 

COSX-  SQRT<1.0-SINX2) 

C0SX2-C0SX*C0SX 

SINCS  «SINX*COSX 

IL0-N‘(L-1) 

IMa-N*(M-1| 

00  12SI»1,N 
IQ-d*l-l)/2 
IF  (I.NE.L)  THEN 
IF  (I.NE.M)  THEN 
IF  d.LT.M)  THEN 
IM-I  +  MQ 

ELSE  IF  (I.GT.M)  THEN 
IM-M-flQ 
ENO  IF 

IF  d.LT.L)  THEN 
IL-I-»^LQ 
ELSE 
IL-L-t-IQ 
ENO  IF 

X  -A|IL)*COSX-AdM)*SINX 
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A(IU-X 
END  IF 
END  IF 

IF  (IC.NE.1)  THEN 
ILR-iLQ-fl 
IMR>IMQ-fl 

X  -  EtVEC(ILR)  *COSX-EIVEC(IMR)  *SINX 
ElVECOMR)  -  ElVECdLRI'SINX  EIVEC(IMR)  *COSX 
EtVEC(ILft)-X 
END  IF 

125  CONTINUE 

X>i2.0*A(LM)*SINCS 
Y  -  A(U)  ■  COSX2  *  A(MM)  ‘SINX2-X 
X  -  A(a)*SINX2  ■¥  A(MM)  *COSX2  *  X 
A(LM)  -  (A(a>- A(MM))  *  SINCS  *  A(LM)  *  (COSX2-SINX2) 
A(a)-Y 
A|MM)«X 
END  IF 

IF  (M.NE.N)  THEN 
M-M-t-1 
GO  TO  60 
ELSE 

IF  (L.NE.N-11  THEN 
L-Lt-l 
GO  TO  55 
ELSE 

IF  (IN0.EQ.1)  THEN 
INO-0 
GO  TO  SO 
ELSE 

IF  (THR.GT.ANRMX)  GO  TO  45 
ENO  IF 
END  IF 
ENO  IF 
END  IF 
IQ—N 

00  185  l«1.N 
IQ-IQ-t'N 
a-U(l*l-l)/2 
JQ-NMI-2) 

DO  185  J-I,N 
JQ>JQ4-N 
MM-J-^U*J-J)/2 
IF  (A(LUXT.A(MM))  THEN 
X-A<UI 
A(LU-A(MM) 

A(MM)-X 
IF  (IC.NE.1)  THEN 
00  180K-1,N 
ILR-IQ-t-K 
IMR-JQ-fK 
X«EIVEC(iLR) 

EIVEC(ILR)-EIVEC(IMR) 

EIVEC(IMR)-X 

180  CONTINUE 
ENO  IF 
END  IF 

185  CONTINUE 
RETURN 
ENO 
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C  PRINT  A  MXN  MATRIX  WITH  COLUMN  AND  ROW  NUMBERING 
C  A  -  MATRIX,  WITH  10  COLUMNS 

C 

C  INDEX-  N-OIMENSIONAL  WORK  VECTOR 
C  THE  MATRIX  A  IS  DIMENSIONED  IN  MAIN  PROGRAM 
C  AS  A  MMAX  BY  NMAX  MATRIX 

C  (M.LE.MMAX.  N.LE.NMAX) 

C 

C  WRITTEN  BY  LEON  BORGMAN,  UNIVERSITY  OF  WYOMING 

c  . . 

DIMENSION  AIMMAX.NMAXI.INOEXINMAX) 

1  FORMAT  <1X.n3,9l12) 

2  FORMAT  (1X.I5,10F12.S) 

9  FORMAT  (1X1 
NR-N>10 
ICB— 9 

DO  3  K-1,N 
INOEX(K)-K 

3  CONTINUE 
4NR-NR-10 

ICB-ICB-t^lO 
IF  (NR.LE.10I  THEN 
ICE -ICB -f  NR- 1 
ELSE 

ICE-ICB-f9 
END  IF 
WRITE  12,9) 

WRITE  (2,1)  (INOEX(K),K-ICB,ICE) 

DOS  IR-1,M 

WRITE  (2,21  IR,(A(IR,iC).IC-ICS,ICE) 

8  CONTINUE 
IF  (NR.GT.10)GOTO4 
RETURN 
END 
C 

SUBROUTINE  FQS0RT|A,IA,N,N05,IB,IT,NI) 

. . 

C  SUBROUTINE  SORTS  THE  COMPONENTS  OF  THE  REAL  VECTOR  A(J), 
C  1  .LE.J.LE.N,  IN  ORDER  OF  INCREASING  SIZE. 

C 

C  AIN)  IS  THE  VECTOR  TO  BE  SORTED,  lA  IS  A  VECTOR  OF  INTEGERS 
C  CO-SORTED  WITH  AIN). 

C  BIS)  IS  A  WORK  ARRAY. 

C  IB(N/5)  AND  IT(N/5)  ARE  INTEGER  WORK  VECTORS  OF  DIMENSION 
C  AT  LEAST  N/5 

C  Nl  IS  A  WORK  INTEGER 
C 

C  PROGRAM  WRITTEN  BY  L.  E.  BORGMAN.  LARAMIE.  WYOMING 
C 

C  THE  SUBROUTINE  USES  SUBSIDIARY  SUBROUTINES:  FMID,  FRANK, 

C  FSWAP,  AND  FSORT. 

c . . . 

DIMENSION  IB(N05),IT(N05),IAIN),A(N),B(5) 

IB(1)  =  1 
IT(1)-N 
NI-1 

10  IF(NI.EQ.O)  RETURN 

CALL  FS0RT|N,N05,A,IA,B,IB.IT.NI) 

GOTO  10 
END 
C 

SUBROUTINE  FSWAPIA.B.IA.IB) 

C  . . 
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T-A 

A-B 

B-T 

I- IA 
lA-IB 
IB-1 
RETURN 
END 

C 

SUBROUTINE  PRANK(N,NS,NT,A,IA) 

C . . 

C  SUBSIDIARY  SUBR  TO  SUBR  FQSORT 

C  . 

DINCNSION  IA(N),A(N) 

NTOT-NT-NS^I 
IF  (NTOT.Ea.1i  RETURN 
NTM-NT-1 
DO  10I-NS,NTM 
IP-1+1 

DO  10  J-IP,NT 

IF  (A(I).GT.A(J))  CALL  FSWAP(A(I).AU).IA(I),IAU)) 
10  CONTINUE 
RETURN 
END 
C 

SUBROUTINE  FMIOIN.NS.NT.A.B.IA) 

C  . 

C  SUBSIDIARY  SUBR  TO  FQSORT 

C  . 

DIMENSION  AIN),B(S),IAIN) 

II - NS 
IS-NT 

12- 11 +(NT-NS)/4 

13- 11  +(NT-NS)/2 

14- 11 +(3*(NT-NS))/4 
B(t)-A(l1) 

BI2)-A|I2) 

B(3)-A(I3) 

B(4)-A(I4) 

B(S.'-A(I5) 

DO  1  J-1,4 

JP-J+1 

DO  1  JJ=JP,5 

IF  (B(J).LE.B(JJ))  GO  TO  1 

Z-BUJ) 

B(JJ)-B(J) 

B|J)-2 
1  CONTINUE 
AA-A(II) 
ll-IA(ll) 

BB-B(3) 

IF(BB.EQ.A(I2))  CALL  FSWAP(AA,A(I2I,II,IA<I2)) 
IF(BB.EQ.A(I3I)  CALL  FSWAP<AA,A<I3),IMA(I3)) 
IF(BB.EQ.A(I4))  CALL  FSWAP(AA,A(I4),II,IA(I4)) 
IF(BB.EQ.A(IS))  CALL  FSWAP(AA,A(I5),II,IA(I5)) 
A(I1)-AA 
IA|I1)-II 
RETURN 
END 
C 

SUBROUTINE  FS0RT(N.N05,A,IA,B.IB,IT,NI) 

C  . 
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DIMENSION  IB(NOS),mNOS),IA(N).A(N).B(S) 
NS«IB(1) 

NT-rr(1) 

NTOT-NT-NS  +  1 
IF(NT0T.GE.10)  GO  TO  10 
CAU  FRANK(N.NS.NT.A.iA) 

NI-NI-1 

IF(NI.GE.1)  GO  TO  17 
IB(1)-0 

mi)-o 

RETURN 

17  DO  18  JK-I.NI 
JP-JK-i-1 
IB(JKI-IB(JP) 
iT(jK)-rryp) 

18  CONTINUE 
IBUPI-O 
rrwPi-o 

RETURN 

10  CONTINUE 

CAa  FMIO(N,NS,NT,A,B,IA) 

JS-NS 

JE-NT 

JM-NS 

IGO-1 

1 1  CONTINUE 
IF(IGO)  13,12,12 

12  CONTINUE 

IF(AgEI.GE.A(JM).ANO.JE.GT.JM)  JE  « JE-1 
IF(A<JE).GE.A(JM).ANOJE.GT.JM)  GOTO  12 
IFUE.LEJM)  GOTO  14 
CAaFSWAP(A(JE),A(JM),IA(JE),IA(JM» 
JM-JE 
14IG0— 1 
GOTO  15 

13  CONTINUE 

IF(A(JSI.LE.AMMI.ANO JS.LTJMI  JS  > JS 1 
IF(A(JS).LE.A(JM).ANDJS.LTJM)  GO  TO  13 
IF  (JS.GEJM)  GOTO  16 
CAUFSWAP(A(JS),A(JM),IA(JS),IA(JM)) 
JM*JS 
16IGO«1 

15IFUE.GT.JS)  GOTO  11 
NI-NI-1 
IBOT-JM-NS 
lUP-NT-JM 

IF(IBOT.EQ.O.OR.IUP.EQ.O)  GO  TO  19 

D0  20  JK-1,NI 

J-NI-JK-t-2 

JJ-J  +  1 

IB(JJ)«IB(J) 

IT(JJ)=mJ) 

20  CONTINUE 
IBID -NS 
IT(1)-JM-1 
IB(2)-JM-t^1 
m2)-NT 
NI-NI-t-2 
RETURN 
19NI-NI-f1 
IFIIUP.EQ.O)  GO  TO  21 
IBID-JMfl 
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RETURN 
21  IBID-NS 
IT(1)-JM-1 
RETURN 
END 


SUBROUTINE  LOP ASSILPSTEP.FLPASS.WIDTH.TIMSER.WRKVEC.OT.FFRACI 

LOW  PASS  FILTER  {TIMSERUhJ-I.LPSTEP}  TO  PASS  CONTENT  AT 
FREQUENCIES  LESS  THAN  FLPASS  AND  EUNNNATE  SIGNAL  AT 
FREQUENCIES  GREATER  THAN  FLPASS WIDTH 

LPSTEP  -  NUMBER  OF  TIME  STEPS.  ASSUMED  TO  BE  FACTORABLE 
INTO  A  PRODUCT  OF  POWERS  OF  2,  3,  AND  5. 

OT  -  TIME  INCREMENT 

FLPASS  -  FREQUENCY  <HERT23  AT  WHICH  LOWER  FREQUENCIES  ARE 
PASSED  WITH  DISTORTION  LESS  THAN  FFRAC. 

WIUTh  -  WIDTH  OF  TRANSITION  ZONE  FOR  RLTER,  IN  HERTZ,  FROM 
FLPASS  TO  A  HIGHER  FREQUENCY  WHERE  THE  SIGNAL  IS 
EUMINATEO  WITH  REMAINDER  LESS  THAT  FFRAC. 

SUBROUTINE  USES  A  GAUSSIAN  SHAPE 

WRITTEN  BY  LEON  BORGMAN,  UNIVERSITY  OF  WYOMING 


DIMENSION  TIMSERILPSTEPl.WRKVECILPSTEP) 

C  •  •  •  COMPUTE  PREUMINARY  CONSTANTS  •  •  * 
0F-1.0/|LPSTEP*0T| 

FMP  -  FLPASS  -f  WIOTH/2.0 

XL02-LPSTEP/2.0 

U-1.0-FFRAC 

CAU  RNORMIU, CUTOFF) 

SIG  -  WI0TH/<2.0»CUT0FF| 

C  •••  FFT  THE  TIME  SERIES  ••• 

CAU  FT23S(LPSTEP,-1.0,0,TIMSER.WRKVEC) 

C  •  •  •  RLTER  THE  FREQUENCY  SEQUENCE  •  *  * 

DO  1  J-1, LPSTEP 
JMmJ-1 
XJM>JM 

IF  (XJM.LT.XL02)  THEN 

C  •••  NOTE:  F-ABSOLUTE  VALUE  OF  FREQUENCY 

F-JM*OF 
ELSE 

F-(LPSTEP-JM)*OF 
END  IF 

IF  (F.LE.FLPASS)  THEN 
WGT«1.0 

ELSE  IF  <F.GT.FLPASS  + WIDTH)  THEN 
WGT-0.0 
ELSE 

Z«(F'FMP)/SIG 
CALL  DFNORM  (Z.WGT) 

END  IF 

TIMSER(J)  >  WGT*TIMSER(J) 

WRKVEC(J)  -  WGT*  WRKVECU) 

1  CONTINUE 

C  •  •  •  INVERSE  FFT  THE  FREQUENCY  SEQUENCE  *  *  * 
CALL  FT235(LPSTEP,1 .0.O.TIMSER.WRKVEC) 
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00  2J-1.LPSTEP 
TIMSERU)  -TtMSERUt/LPSTEP 
2  CONTINUE 

C  . TIMSERU)  NOW  CONTAINS  THE  FILTERED  RESULT  *  ■  * 

RETURN 

END 

C 

SUBROUTINE  FT23S(N,SGN.UST,XR,Xn 
DIMENSION  XRID.XId) 

C . . . 

C  <XRU),  J  - 1  .N)  -  REAL  PART  OF  SEQUENCE  TO  BE  FOURIER  TRANSFORMED. 

C  (XIU).J - 1  ,N) -IMA6.  PART  OF  SEQUENCE  TO  BE  FOURIER  TRANSFORMED. 

C  THE  DISCRETE  FOURIER  TRANSFORM  IS  OF  THE  FORM: 

C  AIM)  -SUM  FOR  J  - 1  TO  N  OF 

C  (XR-fl*XI)*EXP(l>2*PI*S6N*(M-1)*U-1)/N), 

C  FOR  1.LE.M.LE.N 

C  THUS,  SON  OERNES  THE  SIGN  OF  THE  EXPONENT.  AND  IS  EITHER  1 .0 
C  OR -1.0. 

C  A  TABLE  OF  TRIGONOMETRIC  ANO  BIT  REVERSAL  FUNCTIONS  ARE 
C  USED  IN  THE  FFT.  IF  UST-0.  THE  TABLE  IS  COMPUTED.  IF 

C  UST  - 1 .  THE  TABLE  IS  NOT  COMPUTED  AND  IS  PRESUMED  TO  BE 

C  PRESENT  FROM  EARLIER  CALLS  TO  THE  SUBROUTINE. 

C  CONSEQUENTLY,  UST  -0  ON  FIRST  CALL  TO  SUBR. 

C  UST  - 1  FOR  SUBSEQUENT  CALLS  AT  SAME  N  VALUE. 

C 

C  PARAMETER  VALUES: 

C  MAXN  -  LARGEST  VALUE  TO  BE  USED  FOR  N  (LENGTH  OF  SEQUENCE) 

C  MAXN2  -  LARGEST  POWER  OF  2  LESS  THAN  OR  EQUAL  TO  M..XN 

C  MAXN3  -  LARGEST  POWER  OF  3  LESS  THAN  OR  EQUAL  TO  MAXN 

C  MAXNS  -  LARGEST  POWER  OF  S  LESS  THAN  OR  EQUAL  TO  MAXN 

S  MAXN02-MAXN2/2 

C  MAXN03-MAXN3/3 

C  MAXN05-MAXNS/S 

C 

C  COMMON  BLOCK . . . . . . . 

PARAMETER  (MAXN  - 3000,MAXN2  -  2048,MAXN3  -  21 87,MAXN5  -  625) 
PARAMETER  (MAXN02  -  1024,MAXN03  -729,MAXN05- 1 25) 

C  . 

DIMENSION  W(MAXN).WR(MAXN2),WI(MAXN2),WRD(MAXN3),WID(MAXN3) 
DIMENSION  WRT(MAXN5),WIT(MAXN5) 

DIMENSION  IRB(MAXN02I.C0(MAXN02),SI(MAXN02) 

DIMENSION  IRB0(MAXN03),C00(MAXN03),SI0(MAXN03) 

DIMENSION  IRBT(MAXN0S),C0T(MAXN0S),SIT(MAXN05) 

DIMENSION  C01  (MAXNOSI.SII  (MAXN05).C02(MAXN05),SI2(MAXN05) 
DIMENSION  C03(MAXN0S).SI3(MAXN0SI 

COMMON  /FFT/W.WR.WI,WR0.WI0.WRT,WIT,IRB,C0,SI,IRB0.C0D,SI0 
COMMON  /FFT/IRBT,C0T.SIT,C01  ,SI1  ,C02,SI2.C03,SI3 

C  . 

C  PARAMETER  (MAXN«S000,MAXN2-4096.MAXN3-2187.MAXN5-3125) 

C  PARAMETER  (MAXN02  -  2048.MAXN03  -729.MAXN05-625) 

C  NOTE:  THE  ABOVE  COMMON  BLOCK  PARAMETERS  WOULD  ALLOW  THE  USE  OF 
C  ANY  N.LE.5000  OF  THE  FORM  N-(2»*K)*(3*»L)*(5**M)  FOR 
C  INTEGER  K,  L,  AND  M  WITHOUT  CHANGING  DIMENSIONS. 

C  THE  CHOICES  OF  N  WHICH  MATCH  THESE  REQUIREMENTS  ARE: 


c 

2 

3 

4 

5  6  8 

9  10 

c 

12 

15 

16 

18  20 

24  25 

27 

c 

30 

32 

36 

40  45 

48  50 

54 

c 

60 

64 

72 

75  80 

81  90 

96 

c 

100 

108 

120 

125  128 

135 

144 

150 

c 

160 

162 

180 

192  200 

216 

225 

240 

c 

243 

250 

256 

270  288 

300 

320 

324 
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c 

486 

600 

612 

640 

676  600  626  640 

c 

648 

676 

730 

728 

750  768  800  810 

c 

864 

800 

860 

872  1000  1024  1080  1125 

c 

11S2 

1200 

1216 

1250 

1280 

1286 

1350 

1440 

c 

1468 

1600 

1636 

1600 

1620 

1728 

1800 

1876 

c 

1820 

1844 

2000 

2026 

2048 

2160 

2187 

2250 

c 

2304 

2400 

2430 

2500 

2560 

2582 

2700 

2880 

c 

2816 

3000 

3072 

3125 

3200 

3240 

337S 

3456 

c 

3600 

3646 

3760 

3840 

3888 

4000 

4050 

4086 

c 

4320 

4374 

4600 

4608 

4800 

4860 

5000 

C  HOWEVER.  FOR  A  HXED  N.  LESS  COMPUTER  STORAGE  CAN  BE 
C  REQUIRED  IF  THE  FOUOWING  "EXACT*  DIMENSIONING  IS 

C  INTRODUCED: 

C  N2-2**K 

C  N3-3**L 

C  NS-S**M 

C  N02-N2/2 

C  N03-IN3/3 

C  NOS -NS/S 

C  DIMENSION  W<N),WR(N2).WI(N2),WRD|N3),WIO(N3),WRT(NS), 

C  W1T(N5),IRB(N02),C0(N02).SI(N02i.lRBD(N03),C0D(N03). 

C  SI0(NO3),IRBT(NOS),COT(NOS).SIT(NOS).CO1  (NOS). 

C  SI1(N0S),C02(N0S),SI2(N0S).C03(N0S).SI3(N0S) 

C  THE  MAXIMAL  DIMENSIONING  WAS  TAKEN  FROM  THE  LARGEST 
C  POWER  OF  2,  OF  3,  AND  OF  S  WHICH  IS  LESS  THAN  SOOO. 

C 

C  PROGRAMS  WRITTEN  BY  A.  YFANTIS,  UNIV.  NEVADA,  LAS  VEGAS, 

C  AND  LEON  BORGMAN,  UNIV.  WYOMING,  LARAMIE,  WYO. 

C 

. . . . 

NN-N 

K-0 

1  NO-NN-(NN/2)*2 
IF  (ND.NE.O  r'O  TO  2 
K-K-f  1 
NN-NN/2 

IF  (NN.EQ.1)  GOTO  2 
GO  TO  1 
2L-0 

3  NO-NN-(NN/3)*3 
IF  (NO.NE.O)  GO  TO  4 
L-L  +  l 
NN-NN/3 

IF  (NN.EQ.1)  GO  TO  4 
GO  TO  3 
4M-0 

5  NO-NN-INN/5)*5 

IF  (NO.NE.O)  GO  TO  6 
M-M-f  1 
NN-NN/5 

IF  (NN.EQ.il  GO  TO  6 
GOTOS 

6  CONTINUE 

IF  (NN.NE.1)  THEN 
WRITE  (2,100) 

100  FORMAT  ('WARNING:  THE  ASSUMPTIONS  OF  THE  SUBR.  ARE  VIOLATED'/ 
A'  N  IS  NOT  EQUAL  TO  I2**K)*(3**L)*(S**M)') 

RETURN 

ELSE  IF  (K.NE.O.AND.L.EQ.O.AND.M.EQ.O)  THEN 
N02-N/2 

CAaFFT2(K.N,N02,SGN,LIST,XR,XI,IRB,C0.SI) 

ELSE  IF  (K.EQ.O.ANO.L.NE.O.ANO.M.EQ.O)  THEN 
N03-N/3 
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ELSE  IF  (K.EQ.O.AND.L.EQ.O.AND.M.NE.Ot  THEN 
NOS-N/5 

CAa  FFTSIM.N.NOS.SGN.UST.XR.XI.IRBT.COT.SIT.COI  ,SI1 . 

A  C02,SI2.C03,SI3) 

ELSE  IF  <X.NE.O.AND.L.NE.O.ANO.M.EQ.O)  THEN 
N2-2»*K 
N3»3**L 
N02-N2/2 
N03-N3/3 

CAaFFT23(K,L,N,N2,N3,N02,N03,SGN,LIST,XR,XI.IRB.IRB0, 

A  CO,SI.COO,SIO,WR,WI.WRO,WIO,W) 

ELSE  IF  (K.NE.O.AND.UEQ.O.ANO.M.NE.O)  THEN 
N2-2**K 
N5-5**M 
N02-N2/2 
NOS -NS/5 

CAaFFT2S(K.M.N,N2,N5,N02,N0S,SGN.LIST.XR,XI,IRB,IRBT, 

A  C0.SI,COT,SIT,WR.WI,Vy/RT,WIT,W,CO1  .SI1  .C02,SI2,C03,SI3| 

ELSE  IF  (K.Ea.O.AND.L.NE.O.ANO.M.NE.O)  THEN 
N3-3*»L 
NS-5**M 
N03-N3/3 
NOS -NS/5 

CALLFFT3S(L,M,N.N3.NS,N03,N05.SGN,UST.XR.XI.IRBD,IRBT, 

A  COO.SIO,COT,SIT,WRO,WID,WRT.WIT.W,C01  .SI1  ,C02,SI2,C03,SI3) 

ELSE 

N2-2«*K 

N3-3»*L 

NS-S**M 

N02-N2/2 

N03-N3/3 

NOS -NS/S 

CAU  FFT23S(K,UM,N,N2.N3,NS,N02,N03.N0S,SGN,UST,XR,Xt. 

A  IRB,IRBO,IRBT,CO,SI,COO,SIO,COT.SIT.tVR,WI.WRO,WIO,WRT,WiT, 

A  C01,SI1,C02,SI2,C03,SI3.W) 

END  IF 
RETURN 
END 
C 

SUBROUTINE  FFT2(K,N,N02,SGN.UST,XR,XI,IRB,C0,SI) 

C  . 

C  N  -  LENGTH  OF  DATA  SERIES  BEING  TRANSFORMED. 

C  K-LOG2(N). 

C  N02-N/2 

C  SGN  -  1  OR  -1  ACCORDING  AS  THE  EXPONENT  IN  THE  FFT  TRANSFORMATION 

C  IS  POSITIVE  OR  NEGATIVE. 

C  IF  UST  -0,  THE  REVERSEO-BIT  UST  AND  THE  ASSOCIATED  SINES  AND 
C  COSINES  ARE  COMPUTED.  IF  UST  - 1 ,  THE  UST  IS  NOT  COMPUTED. 

C  INSTEAD,  THE  UST  COMPUTED  ON  THE  PREVIOUS  CALL  IS  USED. 

C  XR  -  REAL  PART  OF  THE  INPUT  DATA  SERIES.  INITIALLY  AND  EQUALS  REAL 
C  PART  OF  OUTPUT  ON  RETURN  FROM  THE  SUBR. 

C  XI -IMAG.  PART  OF  DATA  SERIES  ON  INPUT  AND  IS  IMAG.  PART  OF  OUTPUT 
C  ON  RETURN  FROM  THE  SUBR. 

C  IRB-N02-0IMENSI0NAL  VECTOR  OF  REVERSED  BITS  ON  OUTPUT. 

C  C0-N02-DIMENSI0NAL  VECTOR  OF  COSINES  ON  OUTPUT. 

C  SI  -N02-DIMENSIONAL  VECTOR  OF  SINES  ON  OUTPUT. 

C  . . . . . . . . . 

DIMENSION  IRB1 1 1.CO1 1  ),Sin  ),XR1 1  ),Xin ) 

IF  (LIST.EQ.1)G0T0  3 

C  THE  REVERSED-BIT  LIST  AND  TRIG.  FUNCTIONS  ARE  GENERATED  BY  THE 
C  PROCEDURE  GIVED  IN  TABLE  7-VI. 

IRB(1)-0 
DO  1  J-1.K 
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DO  1  1-1,10 
IRB<U-IRB(I)*2 

IF  aLT.IO  IRWUIOI-IRBII)-*- 1 

1  CONTINUE 
FN-N 

W-6.28318S3/FN 
00  21-1,10 
Fm-mB(ll/2 
A-FIR*W 
C0<l)-C0S(A) 

SI(I)-SIN(A) 

2  CONTINUE 

3  CONTINUE 

C  FOR  EACH  COLUMN  OF  THE  FFT  (NC- COLUMN  NUMBER),  CALCULATE 
C  NB-NUMBER  OF  BLOCKS  IN  THE  COLUMN,  LB-LENGTH  OF  BLOCKS,  LB2-HALF 
C  OF  BLOCK  LENGTH. 

00  4NC-1,K 
NB-2**(NC-1) 

LB-N/NB 

LB2-LB/2 

C  THE  BLOCKS  OF  THE  COLUMN  ARE  LOOPED  OVER.  IS -SEQUENCE  NUMBER 
C  AT  THE  START  OF  THE  BLOCK.  IFF  -  SEQUENCE  NUMBER  HALF  WAY  THROUGH 

C  BLOCK. 

00  4IB-1,NB 
C-CO(IB) 

S-SGN*SI(IB) 

IS-<IB-1)*LB-t-1 

IFF-(IB-1)*LB-fLB2 

C  THE  VALUES  IN  THE  NEXT  COLUMN  ARE  COMPUTED  USING  THE  TRIG. 

C  FUNCTIONS  FROM  THE  PRE-GENERATED  UST. 

00  4I-IS,IFF 
l2-l-t-LB2 

QR-XR(I2)*C-XI(I2)*S 

QI-XR(l2)*S-»-XI<t2l*C 

XR(I2)-XR(I>>QR 

XI(I2)-XI(I)-QI 

XR(I)-XR(I)4QR 

XKD-XKD-t-QI 

4  CONTINUE 

C  THE  FFT  COEFFICIENTS  ARE  UNSCRANBLED  INTO  NORMAL  ORDER. 

00  6I-1,N02 
DO  6L-1,2 
IR-IRBIDi-L 
ll-l•f(L-1)■N02 
IF  (IR.LE.il)  GO  TO  6 
2R-XR(IR) 

ZI-XKIR) 

XR(IR)-XR(II) 

XI(IR)-XI(II) 

XR(II)-ZR 
Xl(ll)-Zt 
6  CONTINUE 
RETURN 
END 
C 

SUBROUTINE  FFT3(K,N.N03,SGN,LIST.XR.XI.IRB,C0,SI) 

DIMENSION  IRBd  ),C0(1  ),SI(1  ),XR(1  ).XI(1) 

. . : . 

C  THE  INPUT  TO  THE  SUBROUTINE  CONSISTS  OF  THE  VARIABLES 
C  N,  K,  N03,  SGN,  UST,  XR,  AND  XI.  N  IS  THE  NUMBER  OF 
C  MEASUREMENTS  WHICH  IS  OF  THE  FORM  N  -  3  >  >K,  N03  IS  EQUAL 
C  TO  N/3.  SGN  IS  -1  IN  THE  CASE  OF  THE  FOURIER  TRANSFORM, 

C  AND  -f  1  IN  THE  CASE  OF  THE  INVERSE  FOURIER  TRANSFORM. 
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C  OF  THE  INPUT  COMPLEX  SEQUENCE  OF  NUMBERS.  AFTER  THE 
C  FFTS  IS  PERFORMED,  XR.  XI,  ARE  THE  REAL  AND  IMAGINARY 
PARTS  OF  THE  TRANSFORMED  DATA.  IF  LIST- 0,  THE  REVERSE- 
DIGIT  LIST  AND  THE  ASSOCIATED  SINES  AND  COSINES  ARE 
COMPUTED.  IF  UST-1  THE  UST  IS  NOT  COKN>UTED.  INSTEAD 
THE  LIST  COMPUTED  OF  THE  PREVIOUS  CALL  IS  USED.  IRB, 

CO,  ANO  SI.  ARE  WORK  VECTORS. 


T -6.2831 8S3/3.0 
T1  -2.*T 
FN-N 

IFHJST  .EQ.  1)  GO  TO  3 

IRB(1)-0 

00  1  J-1.K 

I0-3**<J-1I 

00  1  1-1,10 

IRB(I)-IRB(I)*3 

IFU  .EQ.  K)  GO  TO  1 

IRB(UI0I-IRB(I)-«^1 

IRB(U2*IOI-IR8(l|-t-2 

1  CONTINUE 
00  21-1,10 
FIR-IRB(II/FN 
A-nR*T 
CO(l)-COS(A) 

SI(I)-SIN(A) 

2  CONTINUE 

C  FOR  EACH  COLUMN  OF  THE  FFT  <NC -COLUMN  NUMBER!,  NB-NUMBER 
C  OF  BLOCKS  IN  THE  COLUMN,  LB-LENGTH  OF  BLOCKS.  LBS-  A  THIRD 
C  OF  BLOCK  LENGTH. 

3  ci-cosm 

SI  -SGN*SIN(T) 

C2-C0Sm) 

S2-SGN*SINmi 
00  4NC-1.K 
NB»3**(NC-1I 
LB-N/NB 
LB3-LB/3 

C  THE  BLOCKS  OF  THE  COLUMN  ARE  LOOPED  OVER.  IS  -  SEQUENCE  NUMBER 
C  AT  THE  START  OF  THE  BLOCK,  IFF  -  SEQUENCE  NUMBER  A  THIRD  WAY 
C  THROUGH  BLOCK 
00  4IB-1,NB 
C-CO(IB) 

S-SGN*SI(IBI 

IS-(IB-irLB-«-1 

IFF-(IB-1)*LB-»-LB3 

C  THE  VALUES  IN  THE  NEXT  COLUMN  ARE  COMPUTED  USING  THE  TRIG. 

C  FUNCTIONS  FROM  THE  PRE-GENERATED  UST. 

00  4I-IS,IFF 

12- l-t-LB3 

13- U2*LB3 

QR1  -XR<I2I*C-XI<I2)*S 

QI1  -XR(l2l*S-»-XI(l2)*C 

QR2  -  XRII3)  MC  ‘C-S  •  S)-2.  •  XI(I3)  •  S  •  C 

012  -  1C  *  C-S  ’S)  •  XI(I3)  +  2.  ‘S  *  C  •  XR(I3) 

XR(I2)  -  XR(I)  +  Cl  •  QR 1  -S 1  •  Ql  1  +  C2  •  QR2-S2*  QI2 
XI|I2)-XI|I)  +  S1*QR1  +C1*QI1  +S2*QR2  +  C2*QI2 
XR(I3)-XR(I)+C2*'QR1-S2*QI1  +C1 ’QRS-SI'QIS 
XI(I3)  -  XKII + S2  'QR  1  +  C2  'QI1  +  S 1  •QR2  Cl  *012 
XR|I)-XR(I)  +  QR1+QR2 
XI(I)-XI(IKQI1  ■•■012 

4  CONTINUE 

C  THE  FFT  COEFFICIENTS  ARE  INSCRAMBLED  INTO  NORMAL  ORDER 
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0061-1,3 

IR-mB(l)>L 

ll-U(L-1)*N03 

IF  IIR.LE.il)  GO  TO  6 

ZR-XR(IR) 

a-XKIR) 

XR(IRI-XR(II) 

XIOR)-XIUI) 

XR(m-ZR 

xiuii-a 

6  CONTINUE 
RETURN 
ENO 
C 

SUmOUTINE  FFT5(K.N,N0S,SGN,LIST,XR,XI,IRB.C0.SI,C01  ,SI1 , 

+  C02.SI2,C03,SI3) 

OINCNSION  IRB(1  ).C0(1 1.SKI  I.XRd  I.XWI  I.C0 1  ( 1 1. 

+  SI1  (1  ),SI2(1  ).C02(1),C03(1),SI3(1) 

. . 

C  THE  INPUT  TO  THE  SUBROUTINE  CONSISTS  OF  THE  VARIABLES 
C  N,  K,  NOS,  SGN,  LIST,  XR.  XI.  N  IS  THE  NUMBER  OF  MEAS- 
C  UREMENTS  WHICH  IS  OF  THE  FORM  N-5**K,  NOS  IS  EQUAL  TO 
C  N/S.  SGN  IS  -1  IN  THE  CASE  OF  THE  FOURIER  TRANSFORM. 

C  ANO  1  IN  THE  CASE  OF  THE  INVERSE  FOURIER  TRANSFORM. 

C  XR  ANO  XI  ARE  TRE  REAL  ANO  IMAGINARY  PARTS  RESPECTIVELY 
C  OF  THE  INPUT  COMPLEX  SEQUENCE  OF  NUMBERS.  AFTER  THE  FFTS 
C  IS  PERFORMEO  XR.  XI  ARE  THE  REAL  ANO  IMAGINARY  PARTS  OF 

C  THE  TRANSi-ORMEO  OATA.  IF  LIST -0.  THE  REVERSE-0I6IT  LIST 
C  ANO  THE  AS30CIATE0  SINES  AND  COSINES  ARE  COMPUTEO.  IF 
C  UST-1,  THE  UST  IS  NOT  COMPUTED.  INSTEAD  THE  UST  COM- 
C  PUTEO  ON  THE  PREVIOUS  CAU  IS  USED.  IRB,  CO,  SI,  C01 , 

C  SI1 ,  C02,  SI2,  COS,  ANO  SIS,  ARE  WORK  VECTORS. 

. . 

T- 6.283 18S3/S. 

T1  -2.‘T 
T2-3.»T 
T3-4.«T 
FN-N 

IF(UST.EQ.  1)G0T03 
IRB(1)-0 
00  1J-1,K 
ID-5»*|J-1) 

DO  1  1-1,10 
IRB(I)-IRB(I)*5 
IFU  .EQ.  K)  GO  TO  1 
IRB(l•^ID)-iRB(l)•f  1 
IRB(U2*IDI-IRB(l)-f2 
IRB(l-)-3*IO)-IRB(l|-t-3 
IR8(l-^4*IO)-IRB(l)-i-4 
1  CONTINUE 
DO  2  1-1,10 
FIR-IRB(I)/FN 
A-FIR»T 
CO(l)»COSIA) 

SI(I)-SIN(A| 

A1  -2.*A 
C01(I)-C0S(A1) 

SI1(I)-SIN(A1) 

AS-S.oA 

C02(I)-C0S(A2) 

SI2(I)-SIN(A2) 

A3-4.*A 

C03(I)-C0S(A3) 
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2  COHTINUE 

C  FOR  EACH  COLUMN  OF  THE  FFT  (NC  -  COLUMN  NUMBER),  CALCULATE 
C  NB-NUMBER  OF  BLOCKS  IN  THE  COLUMN.  LB-LENGTH  OF  BLOCKS, 

C  LBS-  A  FIFTH  OF  BLOCK  LENGTH 

3  SI -SGN'SlNm 

Cl  -cosm 

S2-S6N*SIN(T1I 


C2-COSm) 

53- S6N*SIN(T2) 

C3-C08(T2) 

54- SGN*SIN(T3I 
C4-C0S(T3) 

00  4NC-1,K 
NB-5**(NC-1) 

LB-N/NB 

LBS-LB/5 

C  THE  BLOCKS  OF  THE  COLUMN  ARE  LOOPED  OVER.  IS -SEQUENCE  NUMBER 
C  AT  THE  START  OF  THE  BLOCK,  IFF  -  SEQUENCE  NUMBER  A  FIFTH  WAY 
C  THROUGH  BLOCK. 

00  4IB-1,NB 

IS-(IB-1)*LB-i-1 

IFF-(IB-1)*LB-»-LB5 

C  THE  VALUES  OF  THE  COLUMN  ARE  COMPUTED  USING  THE  TRIG.  FUNC- 
C  HONS  FROM  THE  PRE  GENERATED  UST. 

D0  4I-IS,IFF 
SAA-SGN*SI(IB) 

SA1  -SGN'SIKIB) 

SA2-SGN*SI2(IB) 

SA3-SGN*SI3(IB) 

12- 1 -t- LBS 

C  THE  FFT  COEFFICIENTS  ARE  UNSCRAMBLED  INTO  NORMAL  ORDER 

13- l-«-2*LBS 

14- I4-3*LBS 

15- |44*LB5 

QR1  -XR(I2)*C0(IB>-XI(I2)*SAA 
Oil  -XR(l2)*SAA-^XI(t2)*CO(IB) 

QR2  -  XR(I3)*C01  (IB)-XI(I3)*SA1 
QI2-XR(I3|*SA1  -t-XKISl'COKIBI 
QR3  -  XR(I4)  *  C02IIB)-XI(I4)  *  SA2 
013  -  XR(I4)  *SA2  *  XI(I4|  *C02(IB) 

QR4  -  XRIISI  "COSdBI-XKIS)  *SA3 
QI4  -  XR(I5)  *SA3 -f  XKISI  *  C03(IBI 

XR(l2>-XR(l)-t-QR1*C1-QI1*S1  -f  QR2*C2-QI2*S24-QR3*C3-QI3*$3-t-QR4’C4-QI4* -t-SA 
XI(I2)  -XKI)  -I-  QR1  *S1  -f  on  *C1  +  QR2*S2  4- QI2*C2  *  QR3*S3  *  QI3*C3 1- QR4*S4 f  014*  +  C4 
XR(t3)  - XR(I)  QR1  •C2-QI1  *S2 QR2*C4-Qi2'S4 4- QR3*C1 -QI3*S1  4^  QR4*C3-QI4' 4- S3 
XI(I3)  -  XKI)  4-  QR1  *S2  4-  on  *C2  4-  aR2*S4  4  Qi2*C4  4-  QR3  *S1  4-  013  'Cl  4  OR4*S3  4  014*  4C3 
XR(I4)-XR(I)4  0R1  'C3-0I1  *S3  4  0R2*C1QI2*S1  40R3*C4-0I3'S440R4*C2-0I4*  4S2 
XKI4)-XI(t)'4-0R1*S3  4  0t1*C3  4  0R2*Sl4  0)2*C1  40R3*S440I3'C440R4*S240I4*  4C2 
XRIIS)  - XRII)  4  0R1  'CA-On  '344  0R2'C3-0I2'S3  4  OR3*C2-OI3*S2  4  OR4'C1-OI4'  4 SI 
XKIS)  -  XKI)  4  OR1  *S4  4  011  *C4  4  aR2  *S3  4  012  'C3  4  OR3  *S2  4  013  *02  4  ORA  'S 1  4  014*  4  Cl 
XRII)  -  XRII)  4  OR  1  4  0R2  4  QR3  4  ORA 
XKI)  -  XKI)  4  on  4  012  4  013  4  014 
4  CONTINUE 

C  THE  FFT  COEFRCIENTS  ARE  UNSCRAMBLED  INTO  NORMAL  ORDER 
D0  6I-1,N05 
DO  6  L-1.5 
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m-muiHL 
il-U(L-1)*N05 
IF  (m.LE.H)  GO  TO  6 
ZR-xnimi 
Zi->XHm) 

XR(m)-XR(H) 

xHim-xidi) 

XR(iU«ZR 
X1«I)-ZI 
6  CONTINUE 
RETURN 
END 
C 

SUBROUTINE  FFT23(K1  ,IC2,N,N1  ,N2.N02.N03.SGN,UST.XR,XI,IRB,IRB0,C0, 
*  SI,COO,SIO,WR.WI.WRD,WIO,W) 

DIMENSION  XR(1 1.XId  I.IRBd  I.IRBOd  I.COd  I.SId  I.CODd  I. 

SlOd  I.WRd  I.  Wl<  1 1.WRO1 1 ),  WlOd  ).Wd ) 

. . . . 

C  THE  SUBROUTINE  FFT23  JOINS  THE  FFT2  AND  THE  FFT3.  WHEN 
C  USED,  BOTH.  FFT2  AND  FFT3,  SUBROUTINES  MUST  BE  PRESENT. 

C  N  IS  THE  NUMBER  OF  MEASUREMENTS  WHICH  IS  OF  THE  FORM 
C  N-(2**K1)M3**K2).  N1  <-2**K1  AND  N2-3**K2,  N02-N1/2, 

C  NOS  -N2/3.  SGN  — 1  IF  THE  FOURIER  TRANSFORM  IS  TO  BE  PER- 
C  FORMED.  IFUST-O.THEREVERSED-DIGITIJST  ANDTHE  AS- 
C  SOCIATED  SINES  AND  COSINES  ARE  COMPUTED.  IF  UST«  1,  THE 
C  UST  IS  NOT  COMPUTED.  INSTEAD  THE  UST  COMPUTED  ON  THE 
C  PREVIOUS  CALL  IS  USED.  XR  IS  THE  REAL  PART  AND  XI  IS  THE 
C  IMAGINARY  PART  OF  THE  COMPLEX  SEQUENCE  TO  BE  TRANSFORMED. 

C  AFTER  THE  FFT23  IS  PERFORMED  XR  AND  XI  ARE  THE  REAL  AND 
C  IMAGINARY  PART  OF  THE  TRANSFORMED  DATA.  IRB.  CO,  AND  SI 
C  ARE  WORK  VECTORS. 

C . . . . . 

TPI-6.28318S307 

XN-N 

OO30l-1,N1 

Flil-1 

TPNF-TPI*FI/XN 
DO  10  J-1.N2 
L-J-1 

WR0(J)-XR(N1*L-«^I) 

WI0(J)-XI(N1*L-t-l) 

10  CONTINUE 

CAaFFT3(K2,N2.N03,SGN.UST,WR0.WID,IRBD,C00,SID) 

D0  20  J-1,N2 

L-J-1 

XJ-J-1 

THET-SGN*TPNF‘XJ 

C-COSaHET) 

S-SIN(THET) 

XRINI  •L-t-l)-WRO(JI*C-WID|J)*S 
XIIN1  •L-*-l)-WRD|J)'S+WID(J)*C 
20  CONTINUE 
30  CONTINUE 
DO  60  l-1,N2 
M-l-1 

D0  40J>1,N1 
WR(J)>XR|N1*M-«-J) 

WI<J)>XI(N1*M-f  J) 

40  CONTINUE 

CALL  FFT2(K1  ,N1  ,N02,SGN.LIST.WR.WI,IRB.C0,SI) 

D0  50  J-’I.NI 
XRCN1*M-t^J)«WR(J) 

XKNIoM-t-JI-WKJI 
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so  CONTINUE 
0O70l-1,N 
W(I)-XR(I) 

70  CONTINUE 
D0  80I-1,N2 
OOSO  J«1,N1 

80  CONTINUE 
OOSOi-I.N 
W(I)-XI(I) 

80  CONTINUE 
D0  95I-1.N2 
00  9S  J-I.NI 

XI(U-1 )  ■  N2 I)  -  W(N  1  *  (I- 1 1  +  Jl 
95  CONTINUE 
RETURN 
ENO 
C 

SUBROUTINE  FFT2S(K1  .K2,N,N1  .N2.N02.NOS,SGN,UST,XR.XI.IRB,IRBO, 

*  CO.SI,COO.SIO,WR,WI,WRD.WIO,W,C01  ,SI1  ,C02.SI2,C03.SI3) 
DIMENSION  XR(1),XI(1),tRB(1|,IRB0(ll 
DIMENSION  C0(1  i.SId  i.COOd  I.SIOdl.WmdI.WId) 

DIMENSION  WRO1 1 ), WIO1 1 ). Wd  ).C0 1  d  ).SI1  d ) 

DIMENSION  C02d  ).C03d  ),SI2d  ).SI3d  I 

. . 

C  THE  SUBROUTINE  FFr25  JOINS  THE  FFT2  AND  THE  FFTS.  HENCE 
C  WHEN  USED  THE  SUBROUTINES  FFT2  AND  FFTS  MUST  BE  PRESENT. 

C  N  IS  THE  NUMBER  OF  MEASUREMENTS,  WHICH  IS  OF  THE  FORM 
C  N-I2**K1)‘|S»*K2).  N1-2**K1,  N2-S<>*IC2.  N02-N1/2,  N05-N2/5, 
C  SGN  —  1,  IF  THE  FOURIER  TRANSFORM  IS  TO  BE  PERFORMED.  IF 

C  UST  >0  THE  REVERSED-DIGIT  UST  AND  THE  ASSOCIATED  SINES  AND 
C  COSINES  ARE  COMPUTED.  IF  UST- 1,  THEUST  IS  NOT  COMPUTED. 

C  INSTEAD  THE  UST  COMPUTED  ON  THE  PREVIOUS  CAU  IS  USED. 

C  XR.  XI  ARE  THE  REAL  AND  THE  IMAGINARY  PARTS  RESPECTIVELY 
C  OF  THE  INPUT  COK.c>JEX  SEQUENCE  OF  NUMBERS.  AFTER  THE  FFT25 
C  IS  PERFORMED  XR  AND  XI  ARE  THE  REAL  AND  IMAGINARY  PARTS 
C  OF  THE  TRANSFORMED  DATA.  IRB,  IRBD,  CO.  COO.  SID,  WR,  Wl, 

C  WRO,  WIO,  W.  C01,  Sd,  C02.  SI2.  C03.  AND  SIS  ARE  WORK 
C  VECTORS. 

. . . . 

TPI- 6.283 185307 
XN-N 

00  30I-1.N1 
FI -1-1 

TPNF=TPfFI/XN 
DO  10  J-1,N2 
L-J-1 

WR0U)-XR(N1*L-«-l) 

WID(J)-XI|N1>L-«-l) 

10  CONTINUE 

CALL  FFT5|K2,N2,N05,SGN,UST.WRD.WID,IRBD.C0D,SID,C01  .SI1 , 
-t-C02,SI2,C03,SI3) 

)0  20  J-1,N2 
L-J-1 
XJ-J-1 

THET-SGN*TPNF*XJ 

C-COSITHET) 

S-SIN(THETI 

XRIN1*L  +  I)-WRDIJ)‘C-WI0IJ)*S 
XIIN 1  •  L  +  II  -  WRD(J)  'S  -t-  WIDU)  •  C 
20  CONTINUE 
30  CONTINUE 
00  60I-1,N2 
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D0  40  J-1,N1 
WRU»-XR(N1*M-*^J) 

WIU)>XliN1*M>J) 

40  CONTINUE 

CAU  FFT2(K1  ,N1  .N02.SGN,UST.WR,Wt.lRB.C0.SI) 

00  80  J-1.N1 
XR(N1*M-f  J)-WRU> 

XKNI'M-t-JI-WIUI 
SO  CONTINUE 
60  CONTINUE 
00  70  I- 1,N 
W(I)-XR(U 
70  CONTINUE 
00  80I-1,N2 
00  80  J«1,N1 

XR(U-1 )  *N2 -t- 1)  -  WINI  *  (1-1 K  J) 

80  CONTINUE 
D0  90I-1.N 
W<I)-XI(I) 

90  CONTINUE 
00  9SI*1,N2 
DO  9S  J-1,N1 

XI((J-1)*N2f  l)*W<N1  *(I-1>+J) 

95  CONTINUE 
RETURN 
END 
C 

SUBROUTINE  FFT3S(K1  ,K2.N.N1  .N2.N03,N0S.S6N.LIST,XR,XI.IRB.IRB0, 
C0.SI,C00.SI0,WR,WI,WR0.WI0.W<C01  ,SI1  ,C02.SI2.C03,SI3) 
DIMENSION  XR(1  ),XI(1 1.IRBII  I.IRBOd  I.COII  I.SId  I.COOd  I 
DIMENSION  SlOd  I.WRd  i.Wld  ),WROd  I.WiOID.Wd  >,C01  d  > 

DIMENSION  SI1  (1  ),C02d  i.SI2d).C03d  ).SI3d  I 


C  THE  SUBROUTINE  FFTSS  JOINS  THE  EFTS  AND  THE  EFTS.  HENCE  WHEN 
C  IS  USED,  THE  SUBROUTINES  FFT3  AND  FFT5  MUST  BE  PRESENT.  N  IS 
C  THE  NUMBER  OF  MEASUREMENTS,  WHICH  IS  OF  THE  FORM 
C  N«(3**K1)MS**K2i  N1 -3**K1,  N2-5**K2,  N03-N1/3,  N05-N2/5, 

C  SON  — 1 ,  IF  THE  FOURIER  TRANSFORM  IS  TO  BE  PE:^F0RME0.  AND  SON  -  *  1 
C  IF  THE  INVERSE  FOURIER  TRANSFORM  IS  TO  BE  PERFORMED.  IF  UST  -0, 

C  THE  REVERSED-DIGIT  UST  AND  THE  ASSOCIATED  SINES  AND  COSINES. 

C  IFUST-1.THEUSTISN0TC0IWUTE0.  INSTEAD  THE  UST  COMPUTED 
C  ON  THE  PREVr  w^  CALL  IS  USED.  XR,  XI,  ARE  THE  REAL  AND  THE 
C  IMAGINARY  PAR'^i  FLSPECTIVELY  OF  THE  INPUT  COMPLEX  SEQUENCE  OF 
C  NUMBERS.  Ar  TER  THE  FFT35  IS  PERFORMED  XR  AND  XI  ARE  THE  REAL 
C  AND  IMAGINARY  PARTS  OF  THE  TRANSFORNKD  DATA.  IRB,  IRBD,  CO.  SI. 

C  COO,  SID,  WR,  Wl,  WRO,  WIO,  W,  COI,  Sd.  C02,  SI2,  C03.  AND  813 
C  ARE  YfORK  VECTORS. 

. . . . 

TPI- 6.2831 85307 
XN-N 

00  30I-1.N1 
FI-t-1 

TPNF-TPI»FI/XN 
DO  10  J-1,N2 
L-J-1 

WR0IJ)-XR(N1*L-fI) 

WI0(J)-XI(N1*L-»-l) 

10  CONTINUE 

CAU  FFT5(K2,N2,N05,SGN,LIST.WR0.Wt0,mB0.C00,SI0,C01  .SI1  ,C02, 
•fSt2,C03,SI3) 

DO  20  J«1,N2 

L«J-1 

XJ-J-1 


B28 


Appendix  6  Prediction  Computer  Program 


C-COS(THET) 

S-SmCTHED 

XWN1  •L  +  I>*WRDU)*C-WIDU)*S 
XI<N1  *1  -f  I)  -  WRDU)  *5  WID(J)  *C 
20  CONTINUE 
30  CONTINUE 
00  60I-1,N2 
M-l-1 

D0  40  J-1,N1 
WRU)«iXfl(N1*M-i^J) 

WIU)-XI(N1*M-f  J| 

40  CONTINUE 

CAU  FFT:<;tK1  ,N1  ,N03,SGN.UST.WR,W(,IRB.C0.SI) 

00  50J-1.N1 
XR<NI*M-fJ)-WR(J) 

XI(N1*M4^J)-WI(J) 

80  CONTINUE 
60  CONTINUE 
00  70I>1.N 
W(I)-XR(I) 

70  CONTINUE 
0O80l>1,N2 
00  80  J-1.N1 

XRKJ- 1 )  *  N2  •»- 1)  -  WIN  1  ‘  (I- 1 1 -t' Jl 
80  CONTINUE 
00  90I-1.N 
W(I)-XI(I) 

90  CONTINUE 
00  9SI«1,N2 
00  95  J«1,N1 

XKU- 1 )  •  N2  + 1)  ■  WIN  I  •  <1- 1 1  +  J) 

95  CONTINUE 
RETURN 
ENO 
C 

SUBROUTINE  FfT235|K1  ,K2.K3.N,N1  ,N2,N3,N02,N03,N05.SGN,UST,XR,XI, 

•f  IRB.IRBO,IRBT,CO,SI,COO,SIO,COT.SIT.WR,WI,WRO,WIO.WRT.WIT,C01  ,SI1 , 
•*-C02,SI2,C03,SI3,W) 

OIMENSION  XRI1 1.XIII  I.IRBII  I.IRBOII  I.IRBTII  I.COII  ),SI|1 ), 

+ coon  ),sioii  i.coTii  ),srrn  i.wRiD.wiiD.wRoiD.wion ), 

WRTI1 1.WITII  I.COI  II  I.SII  n  ).C02n).SI2|1  ).C03|1 ). 

51311 1.WII) 

C . . . . . 

C  THIS  SUBROUTINE  JOINS  THE  FFT2,  THE  FFT3,  AND  THE  FFT5.  WHEN 
C  IT  IS  USED,  THE  SUBROUTINES  FFT2,  FFT3.  AND  FFT5,  MUST  BE  PRESENT. 

C  N  IS  THE  NUMBER  OF  MEASUREMENTS.  N  IS  OF  THE  FORM  N >12* *K1  )* 

C  I3**K2)*|5**K31.  N1«2**K1,  N2-3**K2.  N3-5**K3.  N02-NV2. 

C  N03>N2/3,  N05-N3/5.  SGN  >-1,  IF  THE  FOURIER  TRANSFORM  IS  TO  BE 
C  PERFORMED,  AND  SGN  >  -»^  1  IF  THE  INVERSE  FOURIER  TRANSFORM  IS  TO  BE 
C  PERFORMED.  IF  LIST-O,  THE  REVERSED-OIGIT  UST  AND  THE  ASSOCIATED 
C  SINES  AND  COSINES  ARE  COMPUTED.  IF  LIST  - 1,  THE  LIST  IS  NOT 
C  COMPUTED.  INSTEAD,  THE  LIST  COMPUTED  ON  THE  PREVIOUS  CALL  IS 
C  USED.  XR  IS  THE  REAL  PART  AND  XI  IS  THE  IMAGINARY  PART  OF  THE 
C  COMPLEX  SEQUENCE  TO  BE  TRANSFORMED.  AFTER  THE  FFT235  IS  PERFORMED 
C  XR  AND  XI  ARE  THE  REAL  AND  IMAGINARY  PARTS  OF  THE  TRANSFORMED  DATA. 
C  IRB,  IRBD,  IRBT,  CO.  SI,  COD,  SID,  COT,  SIT,  WR,  Wl,  WRD,  WID,  WRT, 

C  WIT,  C01.  C02,  C03,  SI1,  SI2,  SI3,  AND  W  ARE  WORK  VECTORS. 

C . . . 

TPI- 6.283 185307 
XN>N 
XN2-N2 
XN3-N3 
DO  65  l«1,N1 
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TP1  -(XL*TW)/XN 
TP3-TPTXN3 
D0  30J«1,N2 
XJ-J-1 

TP2  -TP»*XJ/(XN2*XN3)  +  TP1 
DO  10K«1,N3 

WRTK)  -XR(N1  •N2»(K-1 » +  N1  ‘(J-D-l-l) 

Wrr(K)  -  XKN 1  •  N2  •  (K- 1)  +  N 1  *  (J- H  +  0 
10  CONTINUE 

CAU  FFT5(K3,N3.N0S,SGN,UST,WRT,WtT.IRBT.C0T,SIT,C01  .SI1  ,C02,SI2, 
•fC03,SI3) 

D0  20K«1,N3 
XK-K-1 

THET-SGN»TP2*XK 

C>C0S|THET) 

S-SIN(THET) 

XR(N  1  •  N2  •  (K-1 )  +  N1  •  (J- 1 1  +  II  -  WRTIKI  'C-WITIKI  *S 
XMN 1  •  N2  •  IK- 1)  +  N  T  ( J- 1 )  + 1)  -  WIT(K)  *0  +  WRTIK)  •  S 
20  CONTINUE 
30  CONTINUE 
00  60K-1,N3 
KM-K-t 
00  40  J«1,N2 

WROUl  -  XRIN 1  *  N2  *  KM -f  N 1  *  (J- 1 )  + 1) 

WI0(J1-XI(N1*N2»KM  +  N1*U-1»  +  I| 

40  CONTINUE 

CALL  PFT3(K2,N2.N03,SGN.LIST.WR0,WI0.IRB0,C0D,SID) 

00S0J-1.N2 

XJ«J-1 

THET«SGN»TP3*XJ 

C-COSITHET) 

S«SIN(THET) 

XRIN 1  •  N  2  •  KM  +  N 1  •  ( J- 1 )  + 1) «  WROI  J)  •  C- WID(  J)  •  S 
XKN  1  •  N2  •  KM  +  N 1  •  IJ- 1 1  +  II  ■  WRDIJ)  •  S  +  WIOIJ)  •  C 
SO  CONTINUE 
60  CONTINUE 
65  CONTINUE 
D0  90K«1,N3 
KM-K-1 
00  80  J>1,N2 
KJ-J-1 
D0  70I-1,N1 

WR(I)-XRINTN2*KM  +  NTKJ  +  I) 

WKII  =  XKN  1  •  N2  •  KM  +  N 1  •  KJ  +  II 
70  CONTINUE 

CAU  FFT2(K1  ,N1  .N02,SGN.UST,WR,WI,IRB,C0,SI| 

DO  80I-1,N1 

XR|NTN2«KM  +  NTKJ  +  II-WR(I| 

XKN  1  •  N2  •  KM  +  N 1  •  KJ +  11 »  Will} 

80  CONTINUE 
90  CONTINUE 

C  THE  FFT  COEFFICIENTS  ARE  UNSCRAMBLED  INTO  NORMAL  ORDER 
DO  100  '  =  1,N 
WIII-XRIII 
100  CONTINUE 
DO  110  l  =  1,N1 
DO  no  J  =  1,N2 
MJ  =  J-1 

DO  110  K-1,N3 
MK-K-1 

XRI(I-1|*N2*N3  +  N3*MJ  +  KI=WIN1*N2*MK  +  NTMJ  +  II 
110  CONTINUE 
DO  120  l>1,N 
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120  CONTINUE 
00  130I>1,N1 
00  130  J-1,N2 
MJ-J-1 

00  130K-1,N3 
MK-K-1 

Xl«t>  1 )  *N2  *N3 -f  N3  *  MJ  t- K)  -  WIN  1  *  N2  *MK N 1  “MJ l> 

130  CONTINUE 
RETURN 
ENO 
C 

SUBROUTINE  OFNORM  (Z,P) 

C  . 

C  SUBROUTINE  COMPUTES  PROBABILITY  ,P.  THAT  A  STANOARO  NORMAL 
C  RANOOM  VARIABLE  IS  LESS  THAN  OR  EOUAL  TO  Z 
C  THE  METHOO  USEO  IS  FROM:  HANDBOOK  OF  MATH.  FUNCTIONS, 

C  EOITEO  BY  ABRAMOWITZ  AND  STEGUN.  PAGE  932,  EQ.26.6.17, 

C  IN  THE  CHAPTER  BY  ZELEN  ABO  SEVERO 

. . . . * . . . 

zz-z 

IF  (Z.LT.0.0)  ZZ— ZZ 
T»1.0/(1.0>0.2316419*ZZ) 

E— l2Z*2ZI/2.0 
IF  IE.LT.-20.0)  P  =  1.0 
IF  IE.LT.-20.0)  GO  TO  1 
ZX-EXP(E)/2.S0662828 
P  - 1 .0-ZX  •((« 1 .330274429  'T- 1 .82 1 255978)  *T 
+  + 1 .78 1 477937)  *T-0.356563782)  'T  +  .31 9381 530)  *T 
1  IF  IZ.LT.0.0)  P-1.0-P 
RETURN 
ENO 
C 

SUBROUTINE  RNORM(U,Z) 


Appendix  B  Prediction  Computer  Progrem 


B31 


Appendix  C 
Output  Results 


Appendix  C  Output  Results 


THRESHOLD^  0.25 


jjj,vajj)= 

64 

0.09 

JJJ,V(IJJ)= 

65 

0.12 

66 

0.13 

jjj,vajj)= 

67 

0.14 

JJJ,V(IJJ)= 

68 

0.14 

JJJ,V(JJJ)= 

69 

0.13 

JJJ,V(JJJ)= 

70 

0.13 

JJJ,V(jjJ)= 

71 

0.13 

72 

0.14 

JJJ,V(JJJ)= 

73 

0.14 

JJJ,V(JJJ)= 

74 

0.15 

JJJ,V(JJJ)= 

75 

0.15 

jjj,vajj)= 

76 

0.14 

JJJ,V(JJJ)= 

77 

0.12 

JJJ.V(JJJ)= 

78 

0.10 

JJJ,V(JJJ)= 

79 

0.09 

JJJ,V(JJJ)= 

80 

0.09 

JJJ.V(JJJ)= 

81 

0.09 

JJJ,V(JJJ)= 

82 

0.10 

JJJ,V(JJJ)= 

83 

0.10 

JJJ,V(JJJ)= 

84 

0.11 

JJJ,V(JJJ)= 

85 

0.12 

JJJ,V(JJJ)= 

86 

0.13 

JJJ,V(JJJ)= 

87 

0.14 

88 

0.15 

JJJ,V(JJJ)= 

89 

0.16 

JJJ,V(JJJ)= 

90 

0.17 

JJJ,V(JJJ)= 

91 

0.17 

JJJ,V(JJJ)= 

92 

0.18 

JJJ,V(JJJ)= 

93 

0.19 

JJJ,V(JJJ)= 

94 

0.20 

JJJ,V(JJJ)= 

95 

0.21 

JJJ,V(JJJ)= 

% 

0.22 

jjj,vajj)= 

97 

0.23 

JJJ,V(JJJ)= 

98 

0.22 

JJJ,V(JJJ)= 

99 

0.20 

JJJ,V(JJJ)= 

100 

0.17 

JJJ,V(JJJ)= 

101 

0.13 

JJJ,V(JJJ)= 

102 

0.10 

JJJ,V(JJJ)= 

103 

0.07 

JJJ,V(JJJ)= 

104 

0.05 

JJJ,V(JJJ)= 

105 

0.03 

JJJ,V(JJJ)= 

106 

-0.01 

JJJ,V(JJJ)= 

107 

-0.04 

JJJ,V(JJJ)= 

108 

-0.08 

JJJ,V(JJJ)= 

109 

-O.ll 

JJJ,V(JJJ)= 

no 

-0.13 
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jjj.vajj)*  120 

-0.04 

JJJ,V(JJJ)=  121 

-0.05 

JJJ,VaJJ)=  122 

-0.05 

123 

-0.05 

JJJ,V(JJJ)=  124 

-0.04 

JJJ,V(JJJ)=  125 

-0.04 

JJJ,V(JJJ)=  126 

-0.04 

JJJ,V(JJJ)=  127 

•0.03 

JJJ,VaJJ)=  128 

-0.03 

JJJ,V(JJJ)=  129 

-0.02 

JJJ,V0J])=  130 

-0.01 

JJJ,V(JJJ)=  131 

0.01 

132 

0.03 

JJJ,V(JJJ)=  133 

0.06 

JJJ,VaJJ)=  134 

0.07 

JJJ,V(JJJ)=  135 

0.08 

JJJ,VaJJ)=  136 

0.09 

JJJ.V(JJJ)=  137 

0.09 

JJJ.V(JJJ)=  138 

0.10 

JJJ,VaJJ)=  139 

0.10 

JJJ,V(JJJ)=  140 

0.11 

JJJ,V(JJJ)=  141 

0.12 

JJJ,V(JJJ)=  142 

0.13 

JJJ,V(JJJ)=  143 

0.14 

JJJ,V(JJJ)=  144 

0.14 

JJJ,V(JJJ)=  145 

0.13 

JJJ,V(JJJ)=  146 

0.13 

JJJ,V(JJJ)=  147 

0.13 

JJJ,V(JJJ)=  148 

0.13 

JJJ,V(JJJ)=  149 

0.13 

JJJ,V(JJJ)=  150 

0.13 

JJJ,V(JJJ)=  151 

0.12 

JJJ,V(JJJ)=  152 

0.12 

JJJ,V(JJJ)=  153 

0.13 

JJJ,V(JJJ)=  154 

0.14 

JJJ,V(JJJ)=  155 

0.16 

JJJ,V(JJJ)=  156 

0.18 

JJJ,V(JJJ)=  157 

0.18 

JJJ,V(JJJ)=  158 

0.17 

JJJ,V(JJJ)=  159 

0.16 

JJJ,V(JJJ)=  160 

0.15 

JJJ,V(JJJ)=  161 

0.14 

162 

0.14 

JJJ,V(JJJ)=  163 

0.14 

JJJ,V(JJJ)=  164 

0.13 

JJJ,V(JJJ)=  165 

0.12 

JJJ,V(JJJ)=  166 

0.11 

JJJ,V(JJJ)=  167 

0.08 

JJJ,V(JJJ)=  168 

0.06 
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JJJ,V(IJJ)= 

178 

•0.01 

jjj,vajj)= 

179 

0.02 

JJI,V(IJJ)= 

180 

0.04 

jjj.vajj)= 

181 

0.05 

jjj,vajj)= 

182 

0.05 

JJJ,V(JJJ)* 

183 

0.04 

JJJ,V(JJJ)= 

184 

0.04 

jjj,vajj)= 

185 

0.03 

jjj,vajj)= 

186 

0.03 

JJJ,V(JJJ)= 

187 

0.04 

JJJ,V(JJJ)= 

188 

0.05 

jjj,vajj)= 

189 

0.06 

jjj,vajj)= 

190 

0.07 

jjj,vajj)= 

191 

0.08 

JJJ,V(JJJ)* 

192 

0.08 

jjj,vajj)= 

193 

0.09 

JJJ,V(JJJ)* 

194 

0.09 

JJJ.V(JJJ)* 

195 

0.09 

JJJ.V(JJJ)= 

196 

0.09 

JJJ.V(JJJ)= 

197 

0.09 

jjj.vqjj)* 

198 

0.09 

JJJ,V(JJJ)* 

199 

0.09 

JJJ,V(JJJ)* 

200 

0.10 

JJJ.V(JJJ)= 

201 

0.11 

JJJ.V(JJJ)* 

202 

0.12 

JJJ.V(JJJ)* 

203 

0.15 

JJJ,V(JJJ)* 

204 

0.18 

JJJ,V(JJJ)= 

205 

0.22 

JJJ,V(JJJ)* 

206 

0.24 

JJJ,V(JJJ)* 

207 

0.25 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE*  0.25 

JJJ,V(JJJ)=  208  0.25 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE*  0.25 


JJJ,V(JJJ)*  209 

0.25 

JJJ,V(JJJ)*  210 

0.24 

JJJ,V(JJJ)*  211 

0.22 

JJJ,V(JJJ)=  212 

0.20 

JJJ,V(JJJ)*  213 

0.18 

JJJ,V(JJJ)=  214 

0.17 

JJJ,V(JJJ)=  215 

0.17 

JJJ,V(JJJ)*  216 

0.16 

JJJ,V(JJJ)*  217 

0.15 

JJJ,V(JJJ)=  218 

0.13 

JJJ,V(JJJ)=  219 

0.10 

JJJ,V(JJJ)*  220 

0.08 

JJJ,V(JJJ)*  221 

0.06 

JJJ,V(JJJ)*  222 

0.04 

JJJ,V(JJJ)*  223 

0.04 

JJJ,V(JJJ)*  224 

0.05 

C 
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JJJ,V(JJJ)=  234  0.12 

JJJ,V(JJJ)=  235  0.12 

JJJ,V(IJJ)=  236  0.11 
JJJ.VaJJ)-  237  0.11 

JJJ,V(nj)=  238  0.12 
JJJ,VaJJ)=  239  0.13 

JJJ,VaJJ)=  240  0.14 

JJJ.V(JJJ)=  241  0.14 

JJJ,VaJJ)=  242  0.13 

JJJ,V(JJJ)=  243  0.13 

JJJ,V(JJJ)=  244  0.12 

JJJ,V(JJJ)=  245  0.11 

JJJ,V(IIJ)=  246  0.11 
JJJ,V(JJJ)=  247  0.11 

JJJ.V(JJJ)=  248  0.11 

JJJ,V(JJJ)=  249  0.12 

JJJ,V(JJJ)=  250  0.13 

JJI. V(JJJ)=  251  0.14 

JJJ, V(JJJ)=  252  0.15 

JJJ,V(JJJ)=  253  0.16 

JJJ,V(JJJ)=  254  0.17 

JJJ,V(JJJ)=  255  0.17 

JJJ,V(JJJ)=  256  0.15 

JJJ,V(JJJ)=  257  0.14 

JJJ,VaJJ)=  258  0.12 
JJJ,V(JJJ)=  259  0.12 

JJJ,V(JJJ)=  260  0.12 
JJJ,V(JJJ)=  261  0.12 

JJJ,V(JJJ)=  262  0.13 

JJJ,V(JJJ)=  263  0.14 

JJJ,V(JJJ)=  264  0.15 

JJJ,V(JJJ)=  265  0.16 

JJJ,V(JJJ)=  266  0.15 

JJJ,V(JJJ)=  267  0.13 

JJJ,V(JJJ)=  268  0.10 

JJJ,V(JJJ)=  269  0.06 

JJJ,V(JJJ)=  270  0.03 

JJJ,V(JJJ)=  271  0.01 

JJJ,VaJJ)=  272  0.00 

JJJ,V(JJJ)=  273  0.00 

JJJ,V(JJJ)=  274  0.00 

JJJ,V(JJJ)=  275  0.00 

JJJ,V(JJJ)=  276  0.00 

JJJ,V(JJJ)=  277  0.01 

JJJ,V(JJJ)=  278  0.01 

JJJ,V(JJJ)=  279  0.01 

JJJ,V(JJJ)=  280  0.00 

JJJ,V(JJJ)=  281  -0.01 
JJJ,V(JJJ)=  282  -0.02 
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292 

0.09 

JJJ,V(JJJ)= 

293 

0.10 

294 

0.10 

JJJ,V(JJJ)« 

295 

0.10 

296 

0.12 

JJJ.V(JJJ)« 

297 

0.14 

JJJ.VdJJ)* 

298 

0.16 

299 

0.17 

JJJ,V(IJJ)= 

300 

0.18 

JJJ,V(IJJ)= 

301 

0.20 

JJJ,V(JJJ)= 

302 

0.21 

jjj.vau)* 

303 

0.22 

304 

0.22 

jjj.vajj)= 

305 

0.21 

306 

0.20 

JJJ.V(JJJ)= 

307 

0.19 

JJJ,V(JJJ)= 

308 

0.18 

JJJ,V(JJJ)= 

309 

0.17 

JJJ,V(JJJ)= 

310 

0.16 

JJJ,V(JJJ)= 

311 

0.15 

JJJ,V(JJJ)= 

312 

0.14 

JJJ,V(JJJ)= 

313 

0.13 

JJJ.V(JJJ)= 

314 

0.11 

JJJ,V(JJJ)= 

315 

0.08 

JJJ,V(WJ)= 

316 

0.06 

317 

0.03 

JJJ,V(JJJ)= 

318 

0.01 

JJJ,V(JJJ)= 

319 

-0.01 

JJJ,V(JJJ)= 

320 

-0.01 

JJJ,V(JJJ)= 

321 

0.00 

JJJ,V(JJJ)= 

322 

0.00 

JJJ,V(JJJ)= 

323 

0.01 

JJJ,V(JJJ)= 

324 

0.01 

JJJ,V(JJJ)= 

325 

0.01 

JJJ,V(JJJ)= 

326 

0.01 

JJJ,V(JJJ)= 

327 

0.00 

JJJ,V(JJJ)= 

328 

-0.02 

JJJ,V(JJJ)= 

329 

-0.03 

JJJ,V(JJJ)= 

330 

-0.04 

JJJ,V(JJJ)= 

331 

-0.06 

JJJ,V(JJJ)= 

332 

-0.08 

JJJ,V(JJJ)= 

333 

-0.10 

JJJ,V(JJJ)= 

334 

-0.11 

JJJ,V(JJJ)= 

335 

-0.13 

JJJ,V(JJJ)= 

336 

-0.14 

JJJ,V(JJJ)= 

337 

-0.15 

JJJ,V(JJJ)= 

338 

-0.16 

JJJ,V(JJJ)= 

339 

-0.17 

340 

-0.13 
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JJJ,V(IJJ)=  350 

-0.03 

JJJ,VaJJ)=  351 

-0.03 

JJJ.VdJJ)*  352 

•0.03 

JJJ.V(JJJ)=  353 

-0.03 

JJJ,V(JJJ)=  354 

■0.03 

JJJ,V(JJJ)=  355 

-0.02 

JJJ,V(JJJ)=  356 

•0.02 

JJJ,V(JJJ)=  357 

-0.02 

JJJ,V(JJJ)=  358 

-0.01 

JJJ,V(JJJ)=  359 

0.00 

JJJ,V(JJJ)=  360 

0.01 

JJJ.V(JJJ)=  361 

0.02 

JJJ.V(JJJ)=  362 

0.03 

J1J,V(JJJ)=  363 

0.05 

JJJ.V(JJJ)=  364 

0.06 

JJJ.V(JJJ)=  365 

0.06 

JJJ.V(JJJ)=  366 

0.06 

JJJ,V(JJJ)=  367 

0.06 

JJJ.V(JJJ)=  368 

0.06 

JJJ,V(JJJ)=  369 

0.07 

JJJ.V(JJJ)=  370 

0.08 

JJJ,VaJJ)=  371 

0.09 

JJJ,V(JJJ)=  372 

0.11 

JJJ,V(JJJ)=  373 

0.12 

JJJ,V(JJJ)=  374 

0.11 

375 

0.10 

JJJ,V(JJJ)=  376 

0.09 

377 

0.07 

JJJ,V(JJJ)=  378 

0.06 

JJJ,V(JJJ)=  379 

0.06 

JJJ,V(JJJ)=  380 

0.05 

JJJ,V(JJJ)=  381 

0.04 

JJJ,V(JJJ)=  382 

0.03 

JJJ,V(JJJ)=  383 

0.02 

JJJ,V(JJJ)=  384 

0.01 

JJJ,V(JJJ)=  385 

-0.01 

JJJ,V(JJJ)=  386 

-0.03 

JJJ,V(JJJ)=  387 

-0.05 

JJJ,V(JJJ)=  388 

-0.07 

JJJ,V(JJJ)=  389 

-0.08 

JJJ,V(JJJ)=  390 

-0.10 

JJJ,V(JJJ)=  391 

-0.10 

JJJ,V(JJJ)=  392 

-0.10 

JJJ.V(JJJ)=  393 

-0.09 

JJJ.V(JJJ)=  394 

-0.08 

JJJ,V(JJJ)=  395 

-0.07 

JJJ,V(JJJ)=  396 

-0.05 

JJJ,V(JJJ)=  397 

•0.04 

JJJ,V(JJJ)=  398 

-0.03 
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JJJ,VaJJ)=  408  -0.13 
JJJ,VaJJ)*=  409  -0.16 
JJJ,VaJJ)=  410  -0.18 
JJJ,V(JJJ)=  411  -0.19 
JJJ,V(IJJ)=  412  -0.20 
JJJ,V(JJJ)=  413  -0.22 
JJJ,V(JJJ)=  414  -0.22 
JJJ,V(JJJ)=  415  -0.22 
416  -0.21 
JJJ,V(JJJ)=  417  -0.19 
JJJ,V(JJJ)=  418  -0.17 
JJJ,VaJJ)=  419  -0.14 
JJJ,V(JJJ)=  420  -0.11 
JJJ,VaJJ)=  421  -0.08 
JJJ,V(IJJ)=  422  -0.05 
JJJ,V(JJJ)=  423  -0.03 
JJJ.V(JJJ)=  424  -0.03 
JJJ,V(JJJ)=  425  -0.02 
JJJ.V(JJJ)=  426  -0.02 
JJJ,V(JJJ)=  427  -0.01 
JJJ.V(JJJ)=  428  0.00 

JJJ,V(JJJ)=  429  0.01 

JJJ,V(JJJ)=  430  0.01 

JJJ,V(JJJ)=  431  0.01 

JJJ,V(JJJ)=  432  0.00 

JJJ,V(JJJ)=  433  -0.02 
JJJ,V(JJJ)=  434  -0.04 
JJJ,V(JJJ)=  435  -0.06 
JJJ.V(JJJ)=  436  -0.08 
JJJ,V(JJJ)=  437  -0.09 
JJJ,V(JJJ)=  438  -0.10 
JJJ,V(JJJ)=  439  -0.10 
JJJ,V(JJJ)=  440  -0.09 
JJJ,V(JJJ)=  441  -0.08 
JJJ,V(JJJ)=  442  -0.07 
JJJ,V(JJJ)=  443  -0.07 
JJJ,V(JJJ)=  444  -0.06 
JJJ,V(JJJ)=  445  -0.06 
JJJ,V(JJJ)=  446  -0.05 
JJJ,V(JJJ)=  447  -0.05 
JJJ,V(JJJ)=  448  -0.06 
JJJ,V(JJJ)=  449  -0.06 
JJJ,V(JJJ)=  450  -0.06 
JJJ,V(JJJ)=  451  -0.05 
JJJ,V(JJJ)=  452  -0.05 
JJJ,V(JJJ)=  453  -0.04 
JJJ,V(JJJ)=  454  -0.04 
JJJ,V(JJJ)=  455  -0.04 
JJJ,V(JJJ)=  456  -0.04 
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JJJ.VaiJ)*  466 

<0.10 

JJJ,V(JJJ)=  467 

-0.07 

JJJ,V(JJJ)=  468 

-0.05 

JJJ.VaJJ)*  469 

-0.03 

470 

■0.01 

JJJ,V(JJJ)=  471 

0.01 

472 

0.02 

JJJ,VaJJ)=  473 

0.02 

JJJ,V(JJJ)=  474 

0.02 

JJ3,V(JJJ)=  475 

0.02 

JIJ,V(JJJ)=  476 

0.02 

JJJ,V(IJJ)=  477 

0.04 

JJJ,V(JJJ)=  478 

0.07 

JJJ,V(JJJ)=  479 

0.10 

JJJ,V(JJJ)=  480 

0.11 

JJJ,V(JJJ)=  481 

0.11 

JJJ,V(JJJ)=  482 

0.09 

JJJ,V(JJJ)=  483 

0.06 

JJJ,V(JJJ)=  484 

0.02 

JJJ,V(JJJ)=  485 

0.00 

JJJ,V(JJJ)=  486 

-0.03 

IJJ,V(IJJ)=  487 

-0.04 

JJJ,V(JJJ)=  488 

-0.04 

JJJ,V(JJJ)=  489 

-0.02 

JJJ,V(JJJ)=  490 

0.00 

JJJ,V(JJJ)*  491 

0.03 

JJJ,V(JJJ)=  492 

0.06 

JJJ,V(JJJ)=  493 

0.08 

JJJ.V(JJJ)=  494 

0.10 

JJJ.V(JJJ)=  495 

0.11 

JJJ,VaJJ)=  496 

0.11 

JJJ,V(JJJ)=  497 

0.11 

JJJ,V(JJJ)=  498 

0.10 

JJJ,V(JJJ)=  499 

0.09 

JJJ,V(JJJ)=  500 

0.08 

JJJ,V(JJJ)=  501 

0.07 

JJJ,V(JJJ)=  502 

0.09 

JJJ,V(JJJ)=  503 

0.11 

JJJ,V(JJJ)=  504 

0.14 

JJJ,V(JJJ)=  505 

0.15 

JJJ,V(JJJ)=  506 

0.15 

JJJ,V(JJJ)=  507 

0.14 

JJJ,V(JJJ)=  508 

0.11 

JJJ,V(JJJ)=  509 

0.07 

JJJ,V(JJJ)=  510 

0.03 

JJJ,V(JJJ)=  511 

0.00 

JJJ,V(JJJ)=  512 

-0.01 

JJJ,V(JJJ)=  513 

-0.01 

JJJ,V(JJJ)=  514 

0.01 
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JJJ,VaJJ)=  582  0.07 

JJJ.V(JJJ)=  583  0.07 

JJJ.V(JJJ)=  584  0.08 

JJJ.VaiJ)*  585  0.11 

JJJ,V(IJJ)=  586  0.13 

587  0.15 

JJJ.VdJJ)*  588  0.16 

JJJ.VaJJ)*  589  0.17 

JJJ,V(JJJ)=  590  0,17 

JJJ.V(JJJ)=  591  0.15 

JJJ,V(JJJ)=  592  0.14 

JJJ.V(JJJ)=  593  0.13 

JJJ,V(JJJ)=  594  0.13 

JJJ,V(JJJ)=  595  0.13 

JJJ,V(JJJ)=  596  0.14 

JJJ,V(JJJ)=  597  0.15 

JJJ,V(JJJ)=  598  0.15 

JJJ,V(JJJ)=  599  0.15 

JJJ.V(JJJ)=  600  0.15 

JJJ.V(JJJ)=  601  0.13 

JJJ,V(JJJ)=  602  0.11 

JJJ.V(JJJ)=  603  0,08 

JJJ,V(JJJ)=  604  0.06 

605  0.04 

606  0.02 

JJJ.V(JJJ)=  607  0.01 

JJJ,V(JJJ)=  608  0.02 

JJJ,V(JJJ)=  609  0.03 

610  0.04 

JJJ,V(JJJ)=  611  0.05 

JJJ,V(JJJ)=  612  0.05 

JJJ,V(JJJ)=  613  0.06 

JJJ,V(JJJ)=  614  0.06 

JJJ,V(JJJ)=  615  0.05 

JJJ,V(JJJ)=  616  0.04 

JJJ,V(JJJ)=  617  0.03 

JJJ,V(JJJ)=  618  0.04 

JJJ,V(JJJ)=  619  0.05 

JJJ,V(JJJ)=  620  0.07 

JJJ,V(JJJ)=  621  0.08 

JJJ,V(JJJ)=  622  0.09 

JJJ,V(JJJ)=  623  0.09 

JJJ,V(JJJ)=  624  0.09 

JJJ,V(JJJ)=  625  0.07 

JJJ,V(JJJ)=  626  0.05 

JJJ,V(JJJ)=  627  0.03 

JJJ,V(JJJ)=  628  O.OI 
JJJ,V(JJJ)=  629  0.01 

JJJ,V(JJJ)=  630  0.01 
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JJJ,V(JJJ)=  640 

0.15 

JJJ,VaJJ)=  641 

0.15 

JJJ,V(JJJ)=  642 

0.15 

JJJ,VaJJ)=  643 

0.15 

JJJ.V(JJJ)=  644 

0.16 

JJJ,VaJJ)=  645 

0.16 

646 

0.15 

JJJ,VaJJ)=  647 

0.14 

JJJ,V(JJJ)=  648 

0.12 

JJJ.VaJJ)=  649 

0.09 

JJJ.V(JJJ)=  650 

0.07 

JJJ,V(JJJ)=  651 

0.04 

JJJ.V(JJJ)=  652 

0.02 

JJJ,V(JJJ)=  653 

0.00 

JJJ,V(JJJ)=  654 

-0.01 

JJJ,V(JJJ)*  655 

-0.01 

JJJ,V(JJJ)=  656 

-0.01 

JJJ,V(JJJ)=  657 

0.01 

JJJ,V(JJJ)=  658 

0.02 

JJJ,V(JJJ)=  659 

0.03 

JJJ,V(JJJ)=  660 

0.04 

JJJ.V(JJJ)=  661 

0.05 

!  JJJ,V(JJJ)=  662 

0.05 

JJJ,V(JJJ)=  663 

0.04 

JJJ,V(JJJ)=  664 

0.03 

JJJ,V(JJJ)=  665 

0.02 

JJJ,V(JJJ)=  666 

0.01 

JJJ,V(JJJ)=  667 

0.01 

JJJ,V(JJJ)=  668 

0.02 

JJJ,V(JJJ)=  669 

0.02 

JJJ,V(JJJ)=  670 

0.03 

JJJ,V(JJJ)=  671 

0.04 

JJJ,V(JJJ)=  672 

0.03 

JJJ,V(JJJ)=  673 

0.02 

JJJ,V(JJJ)=  674 

0.00 

JJJ,V(JJJ)=  675 

-0.01 

JJJ,V(JJJ)=  676 

-0.03 

JJJ,V(JJJ)=  677 

-0.05 

JJJ,V(JJJ)=  678 

■0.06 

JJJ,V(JJJ)=  679 

-0.07 

JJJ,V(JJJ)=  680 

-0.06 

JJJ,V(JJJ)=  681 

-0.05 

JJJ,V(JJJ)=  682 

-0.03 

JJJ,V(JJJ)=  683 

-0.02 

JJJ,V(JJJ)=  684 

-0.01 

JJJ,V(JJJ)=  685 

0.00 

JJJ,V(JJJ)=  686 

0.00 

JJJ,V(JJJ)=  687 

0.00 

JJJ,V(JJJ)=  688 

-0.01 
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JJJ,V(JJJ)=  698 

0.02 

699 

-0.01 

JJJ,V(JJJ)=  700 

-0.04 

JJJ,V(JJJ)=  701 

-0.06 

JJJ.V(JJJ)=  702 

-0.08 

JJJ,V(JJJ)=  703 

-0.09 

JJJ.V(JJJ)=  704 

-0.09 

JJJ,V(JJJ)=  705 

-0.07 

JJJ,V(JJJ)=  706 

-0.05 

JJJ,V(JJJ)=  707 

-0.03 

JJJ,V(JJJ)=  708 

-0.01 

JJJ,V(JJJ)=  709 

0.01 

JJJ,V(JJJ)=  710 

0.03 

JJJ,V(JJJ)=  711 

0.04 

JJJ,V(JJJ)=  712 

0.04 

JJJ.V(JJJ)=  713 

0.04 

JJJ.V(JJJ)=  714 

0.04 

JJJ,V(JJJ)=  715 

0.05 

JJJ,V(JJJ)=  716 

0.06 

JJJ,V(JJJ)=  717 

0.08 

JJJ,V(JJI)=  718 

0.09 

JJI,V(JJJ)=  719 

0.10 

JJJ,V(JJJ)=  720 

0.10 

JJJ,V(JJJ)=  721 

0.08 

JJJ,V(JJJ)=  722 

0.06 

JJJ,V(JJJ)=  723 

0.04 

JJJ,V(JJJ)=  724 

0.01 

JJJ,V(JJJ)=  725 

-0.02 

JJJ,V(JJJ)=  726 

-0.04 

JJJ,V(JJJ)=  727 

-0.06 

JJJ,V(JJJ)=  728 

-0.07 

JJJ,V(JJJ)=  729 

-0.07 

JJJ,V(JJJ)=  730 

-0.06 

731 

-0.05 

JJJ,V(JJJ)=  732 

-0.04 

JJJ,V(JJJ)=  733 

-0.03 

JJJ.V(JJJ)=  734 

-0.03 

JJJ,V(JJJ)=  735 

-0.03 

JJJ,V(JJJ)=  736 

-0.03 

JJJ,V(JJJ)=  737 

-0.03 

JJJ,V(JJJ)=  738 

-0.03 

JJJ.V(JJJ)=  739 

-0.03 

JJJ,V(JJJ)=  740 

-0.02 

JJJ,V(JJJ)=  741 

-0.01 

JJJ,V(JJJ)=  742 

0.00 

JJJ,V(JJJ)=  743 

O.o; 

JJJ,V(JJJ)=  744 

0.02 

JJJ,V(JJJ)=  745 

0.01 

JJJ,V(JJJ)=  746 

-0.01 
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JJJ,V(JJJ)=  756 

-0.16 

JJJ,V(JJJ)=  757 

-0.17 

JJJ,V(JJJ)=  758 

-0.17 

JJJ,V(JJJ)=  759 

-0.17 

JJJ,V(JJJ)=  760 

-0.17 

JJJ,V(JJJ)=  761 

-0.17 

JJJ,V(JJJ)=  762 

-0.17 

JJJ,V(JJJ)=  763 

-0.16 

JJJ,V(JJJ)=  764 

-0.15 

JJJ.V(JJJ)=  765 

-0.14 

JJJ,V(JJJ)=  766 

-0.12 

767 

-0.10 

JJJ,VaJJ)=  768 

-0.08 

JJJ.V(JJJ)=  769 

-0.07 

JJJ.V(JJJ)=  770 

-0.07 

JJJ.V(JJJ)=  771 

-0.07 

JJJ.V(J;J)=  772 

-0.07 

JJJ.V(JJJ)=  773 

-0.08 

JJJ,V(JJJ)=  774 

-0.08 

775 

-0.08 

JJJ,V(JJJ)=  776 

-0.07 

JJJ.V(JJJ)=  777 

-0.06 

JJJ,V(JJJ)=  778 

-0.04 

JJJ,V(JJJ)=  779 

-0.04 

JJJ,V(JJJ)=  780 

-0.03 

JJJ,V(JJJ)=  781 

-0.03 

782 

-0.04 

JJJ,V(JJJ)=  783 

-0.05 

JJJ,V(JJJ)=  784 

-0.07 

JJJ,V(JJJ)=  785 

-0.08 

JJJ,V(JJJ)=  786 

-0.10 

JJJ,V(JJJ)=  787 

-0.11 

JJJ,V(JJJ)=  788 

-0.11 

JJJ,V(JJJ)=  789 

-0.11 

JJJ,V(JJJ)=  790 

-0.11 

JJJ,V(JJJ)=  791 

-0.10 

JJJ,V(JJJ)=  792 

-0.10 

JJJ,V(JJJ)=  793 

-0.10 

JJJ,V(JJJ)=  794 

-0.11 

JJJ,V(JJJ)=  795 

-0.12 

JJJ,V(JJJ)=  796 

-0.12 

JJJ,V(JJJ)=  797 

-0.13 

JJJ,V(JJJ)=  798 

-0.14 

JJJ,V(JJJ)=  799 

-0.14 

JJJ,V(JJJ)=  800 

-0.12 

JJJ,V(JJJ)=  801 

-0.11 

JJJ,V(JJJ)=  802 

-0.09 

JJJ,V(JJJ)=  803 

-0.07 

JJJ,V(JJJ)=  804 

-0.06 

JJJ,V(JJJ)=  814 

0.03 

JJJ,V(JJJ)=  815 

0.03 

JJJ,VaJJ)=  816 

0.04 

JJJ,V(JJJ)=  817 

0.03 

JJJ,V(JJJ)=  818 

0.03 

JJJ,V(JJJ)=  819 

0.02 

JJJ,V(JJJ)=  820 

0.01 

JJJ,V(JJJ)=  821 

-0.01 

JJJ,V(JJJ)=  822 

-0.01 

JJJ,V(JJJ)=  823 

-0.01 

JJJ,V(JJJ)=  824 

0.00 

JJJ.V(JJJ)=  825 

0.02 

JJJ,V(JJJ)=  826 

0.05 

JJJ,V(JJJ)=  827 

0.07 

JJJ,V(JJJ)=  828 

0.08 

JJJ,V(JJJ)=  829 

0.09 

JJJ,V(JJJ)=  830 

0.09 

JJJ.V(JJJ)=  831 

0.08 

JJJ,V(JJJ)=  832 

0.07 

JJJ,V(JJJ)=  833 

0.06 

JJJ,V(JJJ)=  834 

0.06 

JJJ,V(JJJ)=  835 

0.05 

JJJ,V(JJJ)=  836 

0.05 

JJJ,V(JJJ)=  837 

0.05 

JJJ,V(JJJ)=  838 

0.05 

;JJ,V(JJJ)=  839 

0.04 

JJJ.V(JJJ)=  840 

0.03 

841 

0.01 

JJJ,V(JJJ)=  842 

-0.01 

JJJ,V(JJJ)=  843 

-0.03 

JJJ,V(JJJ)=  844 

-0.04 

JJJ,V(JJJ)=  845 

-0.05 

JJJ,V(JJJ)=  846 

-0.05 

JJJ,V(JJJ)=  847 

-0.04 

JJJ,V(JJJ)=  848 

-0.03 

JJJ,V(JJJ)=  849 

0.00 

JJJ,V(JJJ)=  850 

0.02 

JJJ,V(JJJ)=  851 

0.04 

JJJ,V(JJJ)=  852 

0.05 

JJJ,V(JJJ)=  853 

0.04 

JJJ,V(JJJ)=  854 

0.03 

JJJ,V(JJJ)=  855 

0.02 

JJJ,V(JJJ)=  856 

0.00 

JJJ,V(JJJ)=  857 

■0.01 

JJJ,V(JJJ)=  858 

-0.01 

JJJ,V(JJJ)=  859 

-0.01 

JJJ,V(JJJ)=  860 

-0.01 

JJJ.V(JJJ)=  861 

0.00 

00 

II 

> 

0.01 
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JJJ.V(JJJ)= 

872 

0.03 

JJJ.V(JJJ)= 

873 

0.06 

JJJ,V(JJJ)= 

874 

0.09 

JJJ,V(JJJ)= 

875 

0.12 

JJJ,V(JJJ)= 

876 

0.13 

JJJ,V(JJJ)= 

877 

0.14 

JJJ,V(JJJ)= 

878 

0.13 

JJJ,V(JJJ)= 

879 

0.12 

880 

0.10 

jjj,vajj)= 

881 

0.10 

882 

0.10 

jjj,vajj)= 

883 

0.11 

JJJ,V(JJJ)= 

884 

0.11 

JJJ.V(JJJ)= 

885 

0.12 

JJJ.V(JJJ)= 

886 

0.13 

JJJ.V(JJJ)= 

887 

0.12 

JJJ,V(JJJ)= 

888 

0.11 

JJJ,V(JJJ)= 

889 

0.09 

JJJ,V(JJJ)= 

890 

0.06 

JIJ,V(JJJ)= 

891 

0.03 

JJJ.V(JJJ)= 

892 

0.01 

JJJ.V(JJJ)= 

893 

-0.01 

JJJ.V(JJJ)= 

894 

-0.02 

JJJ,V(JJJ)= 

895 

-0.01 

896 

0.00 

JJJ,V(JJJ)= 

897 

0.01 

898 

0.04 

JJJ.V(JJJ)= 

899 

0.06 

JJJ,V(JJJ)= 

900 

0.09 

JJJ,V(JJJ)= 

901 

0.10 

j;j,v(jjj)= 

902 

0.11 

JJJ,V(JJJ)= 

903 

0.11 

JJJ,V(JJJ)= 

904 

0.11 

JJJ,V(JJJ)= 

905 

0.11 

JJJ,V(JJJ)= 

906 

0.12 

JJJ,V(JJJ)= 

907 

0.13 

JJJ,V(JJJ)= 

908 

0.13 

JJJ,V(JJJ)= 

909 

0.14 

JJJ,V(JJJ)= 

910 

0.13 

JJJ,V(JJJ)= 

911 

0.13 

JJJ,V(JJJ)= 

912 

0.12 

JJJ,V(JJJ)= 

913 

0.10 

JJJ,V(JJJ)= 

914 

0.07 

JJJ,V(JJJ)= 

915 

0.04 

JJJ,V(JJJ)= 

916 

0.02 

JJJ,V(JJJ)= 

917 

0.00 

JJJ,V(JJJ)= 

918 

-0.01 

JJJ,V(JJJ)= 

919 

-0.02 

JJJ,V(JJJ)= 

920 

-0.02 
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JJJ.V(JJJ)=  930 

0.08 

JJJ,VaJJ)=  931 

0.09 

JJJ,V(JJJ)=  932 

0.09 

JJJ,V(JJJ)=  933 

0.09 

JJJ,V(JJJ)=  934 

0.09 

JJJ,V(JJJ)=  935 

0.09 

JJJ,V(JJJ)=  936 

0.08 

JJJ,V(JJJ)=  937 

0.06 

JJJ,V(JJJ)=  938 

0.04 

JJJ,V(JJJ)=  939 

0.01 

JJJ,V(JJJ)=  940 

-0.01 

JJJ.V(JJJ)=  941 

-0.02 

JJJ,V(JJJ)=  942 

-0.03 

JJJ,VaJJ)=  943 

-0.03 

JJJ,V(JJJ)=  944 

-0.02 

JJJ,VaJJ)=  945 

-0.01 

JJJ,V(JJJ)=  946 

-0.02 

JJJ,V(JJJ)=  947 

-0.02 

JJJ,V(JJJ)=  948 

-0.02 

JJJ,V(JJJ)=  949 

-0.02 

JJJ,V(JJJ)=  950 

-0.02 

JJJ,V(JJJ)=  951 

-0.03 

JJJ,V(JJJ)=  952 

-0.04 

JJJ,V(JJJ)=  953 

-0.04 

JJJ,V(JJJ)=  954 

-0.04 

JJJ,V(JJJ)=  955 

-0.03 

JJJ,V(JJJ)=  956 

-0.01 

957 

0.01 

958 

0.03 

JJJ,V(JJJ)=  959 

0.04 

JJJ,V(JJJ)=  960 

0.04 

JJJ,V(JJJ)=  961 

0.04 

JJJ,V(JJJ)=  962 

0.03 

JJJ,V(JJJ)=  963 

0.02 

JJJ.V(JJJ)=  964 

0.01 

JJJ,V(JJJ)=  965 

-0.01 

JJJ,V(JJJ)=  966 

-0.01 

JJJ,V(JJJ)=  967 

-0.01 

JJJ,V(JJJ)=  968 

0.00 

JJJ,V(JJJ)=  969 

0.01 

JJJ,V(JJJ)=  970 

0.01 

JJJ,V(JJJ)=  971 

0.00 

JJJ,V(JJJ)=  972 

0.00 

JJJ,V(JJJ)=  973 

0.00 

JJJ,V(JJJ)=  974 

-0.01 

JJJ,V(JJJ)=  975 

-0.02 

JJJ,V(JJJ)=  976 

-0.04 

JJJ,V(JJJ)=  977 

-0.05 

JJJ,V(JJJ)=  978 

-0.06 
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jjj.vajj)* 

jjj,vajj)= 

JJJ.V(JJJ)= 

JJJ,V(JJJ)* 

JJJ,V(JJJ)= 

JJJ,V(JJJ)* 

jjj.vajj)* 

JJJ,V(JJJ)=* 

JJJ,V(JJJ)= 

jjj.vajj)* 

JJJ.V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ.V(JJJ)= 

JJJ,V(JJJ)= 

JJJ.V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)=: 

JJJ,V(JJJ)* 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 


-0.07 

-0.08 

■0.09 

-0.08 

-0.07 

-0.05 

-0.04 

-0.03 

-0.02 

•0.01 

0.00 

0.00 

-0.01 

-0.01 

-0.02 

-0.03 

-0.03 

-0.03 

-0.02 

-0.02 

-0.03 

-0.04 

-0.05 

-0.06 

-0.07 

-0.07 

-0.07 

-0.06 

-0.04 

-0.01 

0.01 

0.02 

0.03 

0.04 

0.03 

0.02 

0.01 

0.01 

0.00 

-0.01 

-0.03 

-0.04 

-0.04 

-0.05 

-0.06 

-0.08 

-0.10 

-0.11 

-0.12 
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JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ.V(JJJ)= 

jjj.vajj)= 

JJJ,V(JJJ)= 

JJJ,V(JIJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

.  JJJ,V(JJJ)= 
JJJ,V(JJJ)= 
JJJ,V(JJJ)= 
jjj,vajj)= 
JJJ,V(JJJ)= 
JJJ,V(JJJ)= 
JJJ,V(JJJ)= 
JJJ,V(JJJ)= 
JJJ,V(JJJ)= 
JJJ,V(JJJ)= 
JJJ,V(JJJ)= 
JJJ,V(JJJ)= 
JJJ,V(JJJ)= 
JJJ,V(JJJ)= 
JJJ,V(JJJ)= 
JJJ,V(JJJ)= 
JJJ,V(JJJ)= 
JJJ,V(JJJ)= 
JJJ,V(JJJ)= 
JJJ,V(JJJ)= 
JJJ,V(JJJ)= 
JJJ,V(JJJ)= 
JJJ,V(JJJ)= 
JJJ,V(JJJ)= 
JJJ.V(JJJ)= 
JJJ,V(JJJ)= 
JJJ.V(JJJ)= 


1046 

1047 

1048 

1049 

1050 

1051 

1052 

1053 

1054 

1055 

1056 

1057 

1058 

1059 

1060 
1061 
1062 

1063 

1064 

1065 

1066 

1067 

1068 

1069 

1070 

1071 

1072 

1073 

1074 

1075 

1076 

1077 

1078 

1079 

1080 
1081 
1082 

1083 

1084 

1085 

1086 

1087 

1088 

1089 

1090 

1091 

1092 

1093 

1094 


-0.11 

-0.12 

-0.12 

-0.12 

-0.12 

-0.12 

-0.12 

-0.13 

-0.13 

-0.14 

-0.15 

-0.15 

-0.16 

-0.17 

-0.19 

-0.19 

-0.20 

-0.19 

-0.17 

-0.14 

-0.12 

-0.09 

-0.07 

-0.06 

-0.06 

-0.06 

-0.06 

-0.05 

-0.03 

-0.02 

-0.02 

-0.02 

-0.03 

-0.05 

-0.07 

-0.09 

-0.11 

-0.13 

-0.15 

-0.17 

-0.18 

-0.18 

-0.17 

-0.15 

-0.14 

-0.12 

O.ll 

-0.09 

-0.09 
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JJJ,V(JJJ)=  1104 

-0.07 

1105 

-0.09 

JJJ,V(JJJ)=  1106 

-0.11 

JJJ,V(JJJ)=  1107 

-0.13 

JJJ,V(JJJ)=  1108 

-0.15 

JJJ,V(JJJ)=  1109 

-0.15 

JJJ.V(JJJ)=  1110 

-0.15 

jjj,v(jjj)=  nil 

-0.13 

JJJ,V(JJJ)=  1112 

-0.10 

JJJ,VajJ)=  1113 

-0.06 

JJJ,V(JJJ)=  1114 

-0.03 

JJJ,VaJJ)=  1115 

0.00 

JJJ,V(JJJ)=  1116 

0.02 

JJJ,V(JJJ)=  1117 

0.04 

JJJ,V(JJJ)=  1118 

0.04 

JJJ.V(JJJ)=  1119 

0.04 

JJJ,V(JJJ)=  1120 

0.03 

JJJ,V(JJJ)=  1121 

0.02 

JJJ,V(JJJ)=  1122 

0.03 

JJJ,V(JJJ)=  1123 

0.04 

JJJ,V(JJJ)=  1124 

0.07 

JJJ,V(JJJ)=  1125 

0.10 

JJJ.V(JJJ)=  1126 

0.12 

JJJ,V(JJJ)=  1127 

0.13 

JJJ.V(JJJ)=  1128 

0.12 

JJJ,V(JJJ)=  1129 

0.10 

1130 

0.07 

JJJ,V(JJJ)=  1131 

0.04 

JJJ.V(JJJ)=  1132 

0.01 

JJJ,V(JJJ)=  1133 

0.00 

JJJ,V(JJJ)=  1134 

-0.01 

JJJ,V(JJJ)=  1135 

0.00 

JJJ,V(JJJ)=  1136 

0.02 

JJJ,V(JJJ)=  1137 

0.05 

JJJ,V(JJJ)=  1138 

0.07 

JJJ,V(JJJ)=  1139 

0.09 

JJJ,V(JJJ)=  1140 

0.10 

JJJ,V(JJJ)=  1141 

0.11 

JJJ,V(JJJ)=  1142 

0.11 

JJJ,V(JJJ)=  1143 

0.10 

JJJ,V(JJJ)=  1144 

0.09 

JJJ,V(JJJ)=  1145 

0.08 

JJJ,V(JJJ)=  1146 

0.09 

JJJ,V(JJJ)=  1147 

0.10 

JJJ,V(JJJ)=  1148 

0.12 

JJJ,V(JJJ)=  1149 

0.14 

JJJ,V(JJJ)=  1150 

0.16 

JJJ,V(JJJ)=  1151 

0.17 

JJJ,V(JJJ)=  1152 

0.16 
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JJJ.V(JJJ)=  1162 

0.07 

1163 

0.08 

1164 

0.09 

JJJ.V(JJJ)*  1165 

0.09 

1166 

0.09 

1167 

0.10 

JJJ,V(JJJ)*  1168 

0.11 

1169 

0.12 

JJJ.VdJJ)*  1170 

0.14 

JJJ,V(JJJ)=  1171 

0.18 

1172 

0.21 

1173 

0.24 

1174 

0.26 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 

0.26 

AR  SHUTDOWN  STEP  AT  1 174 

1175 

0.28 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE* 

0.28 

AR  SHUTDOWN  STEP  AT  1 175 

1176 

0.27 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE* 

0.27 

AR  SHUTDOWN 

STEP  AT  1176 

JJJ,V(JJJj*  1177 

0.24 

JJJ,V(JJJ)*  1178 

0.20 

JJJ,V(JJJ)*  1179 

0.16 

JJJ,V(JJJ)*  1180 

0.14 

JJJ,V(JJJ)*  1181 

0.13 

JJJ,V(JJJ)«  1182 

0.13 

JJJ,V(JJJ)*  1183 

0.15 

JJJ,V(JJJ)=  1184 

0.17 

JJJ,V(JJJ)*  1185 

0.20 

JJJ,V(JJJ)=  1186 

0.23 

JJJ,V(JJJ)*  1187 

0.24 

JJJ,V(JJJ)*  1188 

0.24 

JJJ,V(JJJ)=  1189 

0.24 

JJJ,V(JJJ)=  1190 

0.23 

JJJ,V(JJJ)*  1191 

0.22 

JJJ,V(JJJ)=  1192 

0.21 

JJJ,V(JJJ)=  1193 

0.21 

JJJ,V(JJJ)=  1194 

0.23 

JJJ,V(JJJ)=  1195 

0.25 

JJJ,V(JJJ)*  1196 

0.27 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE* 

0.27 

JJJ,V(JJJ)=  1197 

0.30 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE* 

0.30 

JJJ,V(JJJ)*  1198 

0.32 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE* 

0.32 

KRIGING  SHUTDOWN  STEP  AT  1 198 

JJJ,V(JJJ)*  1199 

0.34 

1  EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE* 

0.34 

Appandix  C  Output  Raxults 


KRIGING  SHUTDOWN  STEP  AT  1202 

AR  SHUTDOWN  STEP  AT  1202 

JJJ.V(JJJ)=  1203 

0.23 

JJJ.VdJJ)*  1204 

0.21 

1205 

0.19 

JJJ,V(JJJ)=  1206 

0.18 

JJJ,V(JJJ)»  1207 

0.19 

JJJ.VaJJ)*  1208 

0.21 

JJJ.VaJJ)*  1209 

0.22 

1210 

0.23 

JJJ,V(JJJ)*  1211 

0.24 

JJJ,V(IJJ)=  1212 

0.23 

JJJ,VaJJ)=  1213 

0.21 

JJJ.V(JJJ)=  1214 

0.18 

JJJ,V(JJJ)=  1215 

0.14 

JJJ,V(JJJ)=  1216 

0.12 

JJJ,V(JJJ)=  1217 

0.11 

JJJ,V(JJJ)=  1218 

0.11 

JJJ,V(JJJ)=  1219 

0.11 

JJJ.V(JJJ)=  1220 

0.12 

JJJ,V(JJJ)=  1221 

0.14 

JJJ.V(JJJ)=  1222 

0.15 

JJJ,V(JJJ)=  1223 

0.17 

JJJ,V(JJJ)=  1224 

0.17 

JJJ,V(JJJ)=  1225 

0.17 

1226 

0.14 

JJJ,V(JJJ)=  1227 

0.13 

JJJ,V(JJJ)=  1228 

0.11 

JJJ,V(JJJ)=  1229 

0.10 

JJJ,V(JJJ)=  1230 

0.10 

JJJ,V(JJJ)=  1231 

0.12 

JJJ,V(JJJ)=  1232 

0.13 

JJJ,V(JJJ)=  1233 

0.14 

JJJ,V(JJJ)=  1234 

0.16 

JJJ,V(JJJ)=  1235 

0.16 

JJJ,V(JJJ)=  1236 

0.16 

JJJ,V(JJJ)=  1237 

0.15 

JJJ,V(JJJ)=  1238 

0.13 

JJJ,V(JJJ)=  1239 

0.12 

JJJ,V(JJJ)=  1240 

0.11 

JJJ,V(JJJ)=  1241 

0.11 

JJJ,V(JJJ)=  1242 

0.11 

JJJ,V(JJJ)=  1243 

0.11 

JJJ,V(JJJ)=  1244 

0.12 

JJJ,V(JJJ)=  1245 

0.12 

JJJ,V(JJJ)=  1246 

0.12 

JJJ,V(JJJ)=  1247 

0.12 

JJJ,V(JJJ)=  1248 

0.11 

JJJ,V(JJJ)=  1249 

0.09 
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JJJ,V(JJJ)=  1317  0.07 

JJJ,V(JJJ)=  1318  0.09 

JJJ,V(JJJ)=  1319  0.08 

JJJ,V(JJJ)=  1320  0.06 

JJJ.V(JJJ)=  1321  0.02 

JJJ,V(JJJ)=  1322  -0.04 
JJJ,V(JJJ)=  1323  -0.09 
JJJ.V(JJJ)=  1324  -0.14 

1325  -0.17 
JJJ.V(JJJ)=  1326  -0.18 
JJJ,V(JJJ)=  1327  -0.17 
JJJ,V(JJJ)=  1328  -0.14 
JJJ,V(JJJ)=  1329  -0.11 
JJJ,V(JJJ)=  1330  -0.08 
JJJ.V(JJJ)=  1331  -0.06 

JJI, V(JJJ)=  1332  -0.04 

JJJ. V(JJJ)=  1333  -0.03 

JJJ,V(JJJ)=  1334  -0.03 
JJJ,V(JJJ)=  1335  -0.02 
JJJ,V(JJJ)=  1336  0.00 

1337  0.04 

JJJ,V(JJJ)=  1338  0.08 

JJJ,V(JJJ)=  1339  0.13 

JJJ,V(JJJ)=  1340  0.17 

JJJ,V(JJJ)=  1341  0.19 

JJJ,V(JJJ)=  1342  0.20 

JJJ,V(JJJ)=  1343  0.18 

JJJ,V(JJJ)=  1344  0.15 

JJJ.V(JJJ)=  1345  0.10 

JJJ,V(JJJ)=  1346  0.04 

1347  -0.01 
JJJ,V(JJJ)=  1348  -0.05 
JJJ,V(JJJ)=  1349  -0.07 
JJJ,V(JJJ)=  1350  -0.07 
JJJ,V(JJJ)=  1351  -0.06 
JJJ,V(JJJ)=  1352  -0.05 
JJJ,V(JJJ)=  1353  -0.04 
JJJ,V(JJJ)=  1354  -0.04 
JJJ,V(JJJ)=  1355  -0.05 
JJJ,V(JJJ)=  1356  -0.06 
JJJ,V(JJJ)=  1357  -0.08 
JJJ,V(JJJ)=  1358  -0.11 
JJJ,V(JJJ)=  1359  -0.14 
JJJ,V(JJJ)=  1360  -0.16 
JJJ,V(JJJ)=  1361  -0.16 
JJJ,V(JJJ)=  1362  -0.13 
JJJ,V(JJJ)=  1363  -0.08 
JJJ,V(JJJ)=  1364  -0.02 
JJJ,V(JJJ)=  1365  0.02 
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1375  -0.14 
JJJ,V(JJJ)=  1376  -0.10 
JJJ,V(JJJ)=  1377  -0.05 
JJJ,V(JJJ)=  1378  -0.01 
JJJ,V(JJJ)=  1379  0.01 

JJJ.V(JJJ)=  1380  0.02 
1381  0.01 

JJJ,V(JJJ)=  1382  -0.01 
JJJ,V(JJJ)=  1383  -0.03 
JJJ,V(JJJ)=  1384  -0.03 
JJJ,V(JJJ)=  1385  -0.02 
JJJ,V(JJJ)=  1386  0.02 

JJJ,V(JJJ)=  1387  0.07 

JJJ,V(JJJ)=  1388  0.12 

JJJ,V(JJJ)=  1389  0.18 

JJJ,V(JJJ)=  1390  0.21 

JJJ,V(JJJ)=  1391  0.22 

JJJ,V(JJJ)=  1392  0.20 

JJJ,V(JJJ)=  1393  0.17 

JJJ.V(JJJ)=  1394  0.12 

JJJ,V(JJJ)=  1395  0.06 

JJJ,V(JJJ)=  1396  0.02 

JJJ,V(JJJ)=  1397  -0.01 
JJJ,V(JJJ)=  1398  -0.03 
1399  -0.03 
JJJ,V(JJJ)=  1400  -0.01 
JJJ,V(JJJ)=  1401  0.01 

JJJ,V(JJJ)=  1402  0.04 

JJJ,V(JJJ)=  1403  0.05 

JJJ,V(JJJ)=  1404  0.04 

JJJ,V(JJJ)=  1405  0.02 

JJJ,V(JJJ)=  1406  -0.01 
JJJ,V(JJJ)=  1407  -0.04 
JJJ,V(JJJ)=  1408  -0.06 
JJJ,V(JJJ)=  1409  -0.07 
JJJ,V(JJJ)=  1410  -0.06 
JJJ,V(JJJ)=  1411  -0.03 
JJJ,V(JJJ)=  1412  0.02 

JJJ,V(JJJ)=  1413  0.07 

JJJ,V(JJJ)=  1414  0.1 1 

JJJ,V(JJJ)=  1415  0.13 

JJJ,V(JJJ)=  1416  0.13 

JJJ,V(JJJ)=  1417  0.11 

JJJ,VajJ)=  1418  0.06 

JJJ,VaJJ)=  1419  0.01 

JJJ,V(JJJ)=  1420  -0.04 
JJJ,V(JJJ)=  1421  -0.08 
JJJ,V(JJJ)=  1422  -0.10 
JJJ,V(JJJ)=  1423  -0.11 
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JJJ,V(JJJ)=  1433 

-0.06 

JJJ,V(JJJ)=  1434 

-0.04 

JJJ,V(JJJ)=  1435 

-0.01 

JJJ,V(JJJ)=  1436 

0.04 

JJJ,V(JJJ)=  1437 

0.09 

JJJ.V(JJJ)=  1438 

0.15 

JJJ,V(JJJ)=  1439 

0.19 

JJJ,V(JJJ)=  1440 

0.21 

JJJ.V(JJJ)=  1441 

0.21 

JJJ,V(JJJ)=  1442 

0.19 

JJJ,V(JJJ)=  1443 

0.15 

JJJ,V(JJJ)=  1444 

0.11 

JJJ,V(JJJ)=  1445 

0.08 

JJJ,V(JJJ)=  1446 

0.06 

JJJ,V(JJJ)=  1447 

0.05 

JJJ,V(JJJ)=  1448 

0.06 

JJJ,V(JJJ)=  1449 

0.07 

JJJ,V(JJJ)=  1450 

0.10 

JJJ,V(JJJ)=  1451 

0.13 

JJJ,V(JJJ)=  1452 

0.15 

JJJ,V(JJJ)=  1453 

0.16 

JJJ,V(JJJ)=  1454 

0.15 

JJJ,V(JJJ)=  1455 

0.12 

JJJ,V(JJJ)=  1456 

0.09 

JJJ,V(JJJ)=  1457 

0.06 

JJJ,V(JJJ)=  1458 

0.04 

JJJ,V(JJJ)=  1459 

0.04 

JJJ,V(JJJ)=  1460 

0.06 

JJJ,V(JJJ)=  1461 

0.09 

JJJ,V(JJJ)=  1462 

0.13 

JJJ,V(JJJ)=  1463 

0.17 

JJJ,V(JJJ)=  1464 

0.21 

JJJ,V(JJJ)=  1465 

0.22 

JJJ,V(JJJ)=  1466 

0.21 

JJJ,V(JJJ)=  1467 

0.19 

JJJ,V(JJJ)=  1468 

0.16 

JJJ,V(JJJ)=  1469 

0.13 

JJJ,V(JJJ)=  1470 

0.11 

JJJ,V(JJJ)=  1471 

0.09 

JJJ,V(JJJ)=  1472 

0.10 

JJJ,V(JJJ)=  1473 

0.11 

JJJ,V(JJJ)=  1474 

0.12 

JJJ,V(JJJ)=  1475 

0.14 

JJJ,V(JJJ)=  1476 

0.15 

JJJ,V(JJJ)=  1477 

0.15 

JJJ,V(JJJ)=  1478 

0.14 

JJJ,V(JJJ)=  1479 

0.12 

JJJ,V(JJJ)=  1480 

0.09 

JJJ,V(JJJ)=  1481 

0.05 
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JJJ,V(JJJ)=  1491 

0.16 

JJJ,V(JJJ)=  1492 

0.14 

JJJ,V(JJJ)=  1493 

0.12 

JJJ,V(«J)=  1494 

0.10 

JJJ,V(JJJ)=  1495 

0.09 

JJJ,V(JJJ)=  1496 

0.09 

JJJ,V(JJJ)=  1497 

0.10 

JJJ.VajJ)=  1498 

0.12 

JJJ.V(JJJ)=  1499 

0.13 

JJJ.V(JJJ)=  1500 

0.15 

JJJ,V(JJJ)=  1501 

0.16 

JJJ,V(JJJ)=  1502 

0.16 

JJJ,V(JJJ)=  1503 

0.15 

JJJ.V(JJJ)=  1504 

0.13 

JJJ,V(JJJ)=  1505 

0.09 

JJJ,V(JJJ)=  1506 

0.06 

JJJ,V(JJJ)=  1507 

0.03 

JJJ,V(JJJ)=  1508 

0.01 

JJJ,V(JJJ)=  1509 

0.01 

JJJ,V(JJJ)=  1510 

0.03 

JJJ.V(JJJ)=  1511 

0.06 

JJJ,V(JJJ)=  1512 

0.10 

JJJ,V(JJJ)=  1513 

0.13 

JJJ,V(JJJ)=  1514 

0.15 

JJJ,V(JJJ)=  1515 

0.16 

1516 

0.16 

JJJ,V(JJJ)=  1517 

0.16 

JJJ,V(JJJ)=  1518 

0.15 

JJJ,V(JJJ)=  1519 

0.15 

JJJ,V(JJJ)=  1520 

0.15 

JJJ,V(JJJ)=  1521 

0.17 

JJJ,V(IJJ)=  1522 

0.19 

JJJ,V(JJJ)=  1523 

0.21 

JJJ,V(JJJ)=  1524 

0.23 

JJJ,V(JJJ)=  1525 

0.25 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 

0.25 

JJJ,V(JJJ)=  1526 

0.27 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 

0.27 

JJJ,V(JJJ)=  1527 

0.27 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 

0.27 

JJJ,V(JJJ)=  1528 

0.25 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 

0.25 

JJJ,V(JJJ)=  1529 

0.23 

JJJ,V(JJJ)=  1530 

0.20 

JJJ,V(JJJ)=  1531 

0.16 

JJJ,V(JJJ)=  1532 

0.12 

JJJ,V(JJI)=  1533 

0.09 

JJJ,V(JJJ)=  1534 

0.06 

JJJ,V(JJJ)=  1535 

0.06 
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JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ.V(JJJ)= 

jjj,vajj)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

jjj.vajj)= 

JJJ.V(JJJ)= 

JJJ,V(JjJ)= 

JJJ.V(JJJ)= 

JJJ,V(IJJ)= 

jjj.va«)= 

jjj,vajj)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

IJJ.V(JJJ)= 

JJJ.V(JJJ)= 

JJJ,V(JJJ)= 

JJJ.V(JJJ)= 

JJJ.V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ.V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ.V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 

JJJ,V(JJJ)= 


1545 

1546 

1547 

1548 

1549 

1550 

1551 

1552 

1553 

1554 

1555 

1556 

1557 

1558 

1559 

1560 

1561 

1562 

1563 

1564 

1565 

1566 

1567 

1568 

1569 

1570 

1571 

1572 

1573 

1574 

1575 

1576 

1577 

1578 

1579 

1580 

1581 

1582 

1583 

1584 

1585 

1586 

1587 

1588 

1589 

1590 

1591 

1592 

1593 


0.11 

0.12 

0.13 

0.14 

0.15 

0.16 

0.17 

0.16 

0.15 

0.14 

0.11 

0.07 

0.04 

0.00 

-0.02 

-0.03 

-0.03 

-0.03 

-0.02 

0.00 

0.00 

O.OQ 

0.00 

0.00 

0.02 

0.03 

0.06 

0.08 

0.11 

0.13 

0.17 

0.20 

0.22 

0.23 

0.22 

0.20 

0.17 

0.14 

0.11 

0.10 

0.09 

0.10 

0.11 

0.13 

0.14 

0.14 

0.14 

0.14 

0.14 


KRIGING  SHUTDOWN  STEP  AT  1600 

AR  SHUTDOWN  STEP  AT  1600 

JJJ,VaJJ)=  1601 

0.30 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE* 

0.30 

KRIGING  SHUTDOWN  STEP  AT  1601 

AR  SHUTDOWN  STEP  AT  1601 

JJJ,V(JJJ)=  1602 

0.30 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE* 

0.30 

KRIGING  SHUTDOWN  STEP  AT  1602 

AR  SHUTDOWN  STEP  AT  1602 

JJJ,V(JJJ)=  1603 

0.29 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE* 

0.29 

KRIGING  SHUTDOWN  STEP  AT  1603 

AR  SHUTDOWN  STEP  AT  1603 

JJJ,V(JJJ)=  1604 

0.27 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE* 

0.27 

KRIGING  SHUTDOWN  STEP  AT  1604 

AR  SHUTDOWN 

STEP  AT  1604 

JJJ,V(JJJ)=  1605 

0.24 

JJJ,V(JJJ)=  1606 

0.21 

JJJ,V(JJJ)=  1607 

0.18 

JJJ,V(IJJ)=  1608 

0.15 

JJJ,V(JJJ)=  1609 

0.13 

JJJ,V(JJJ)=  1610 

0.12 

JJJ,V(JJJ)=  1611 

0.12 

JJJ,V(JJJ)=  1612 

0.13 

JJJ,V(JJJ)*  1613 

0.15 

JJJ,V(JJJ)*  1614 

0.16 

JJJ,V(JJJ)=  1615 

0.16 

JJJ,y(JJJ)=  1616 

0.16 

JJJ,V(JJJ)=  1617 

0.15 

JJJ,V(JJJ)=  1618 

0.14 

JJJ,V(JJJ)=  1619 

0.13 

JJJ,V(JJJ)=  1620 

0.13 

JJJ,V(JJJ)=  1621 

0.14 

JJJ,V(JJJ)=  1622 

0.17 

JJJ,V(JJJ)=  1623 

0.20 

JJJ,V(JJJ)=  1624 

0.23 

JJJ,V(JJJ)=  1625 

0.27 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE* 

0.27 

JJJ,V(JJJ)=  1626 

0.29 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE* 

0.29 

KRIGING  SHUTDOWN  STEP  AT  1626 

AR  SHUTDOWN  STEP  AT  1626 

JJJ,V(JJJ)*  1627 

0.31 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE* 

0.31 

KRIGING  SHUTDOWN  STEP  AT  1627 

AR  SHUTDOWN  STEP  AT  1627 

JJJ,V(JJJ)=  1628 

0.31 
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KRIGING  SHUTDOWN  STEP  AT 

1630 

JJJ,V(JJJ)=  1631 

0.27 

EXCEEDANCE  OF  THRESHOLD 

FOUND, 

KRIGING  SHUTDOWN  STEP  AT 

1631 

JJJ,V(JJJ)=  1632 

0.25 

EXCEEDANCE  OF  THRESHOLD 

FOUND, 

JJJ,V(JJJ)=  1633 

0.24 

JJJ,V(JJJ)=  1634 

0.23 

JJJ.V(JJJ)=  1635 

0.23 

JJJ,V(J1I)=  1636 

0.25 

JJJ.V(JJJ)=  1637 

0.27 

EXCEEDANCE  OF  THRESHOLD  FOUND, 
KRIGING  SHUTDOWN  STEP  AT  1637 
AR  SHUTDOWN  STEP  AT  1637 
JJJ.V(JJJ)=  1638  0.29 

EXCEEDANCE  OF  THRESHOLD  FOUND. 
KRIGING  SHUTDOWN  STEP  AT  1638 
AR  SHUTDOWN  STEP  AT  1638 
JJJ,V(JJJ)=  1639  0.31 

EXCEEDANCE  OF  THRESHOLD  FOUND, 
KRIGING  SHUTDOWN  STEP  AT  1639 
AR  SHUTDOWN  STEP  AT  1639 
JJJ,V(JJJ)=  1640  0.32 

EXCEEDANCE  OF  THRESHOLD  FOUND, 
KRIGING  SHUTDOWN  STEP  AT  1640 
AR  SHUTDOWN  STEP  AT  1640 
JJJ,V(JJJ)=  1641  0.31 

EXCEEDANCE  OF  THRESHOLD  FOUND, 
KRIGING  SHUTDOWN  STEP  AT  1641 
AR  SHUTDOWN  STEP  AT  1641 
JJJ,V(JJJ)=  1642  0.30 

EXCEEDANCE  OF  THRESHOLD  FOUND, 
KRIGING  SHUTDOWN  STEP  AT  1642 
AR  SHUTDOWN  STEP  AT  1642 
JJJ,V(JJJ)=  1643  0.29 

EXCEEDANCE  OF  THRESHOLD  FOUND, 
KRIGING  SHUTDOWN  STEP  AT  1643 
AR  SHUTDOWN  STEP  AT  1643 
JJJ,V(JJJ)=  1644  0.29 

EXCEEDANCE  OF  THRESHOLD  FOUND, 
KRIGING  SHUTDOWN  STEP  AT  1644 
AR  SHUTDOWN  STEP  AT  1644 
JJJ,V(JJJ)=  1645  0.30 

EXCEEDANCE  OF  THRESHOLD  FOUND, 
KRIGING  SHUTDOWN  STEP  AT  1645 
AR  SHUTDOWN  STEP  AT  1645 
JJJ,V(JJJ)=  1646  0.30 

EXCEEDANCE  OF  THRESHOLD  FOUND, 
KRIGING  SHUTDOWN  STEP  AT  1646 


VALUE= 


VALUE= 


VALUE= 


VALUE= 


VALUE= 


VALUE= 


VALUE= 


VALUE= 


VALUE= 


VALUE= 


VALUE = 


VALUE = 


0.27 

0.25 

0.27 

0.29 

0.31 

0.32 

0.31 

0.30 

0.29 

0.29 

0.30 

0.30 
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1649  0.36 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 

0.36 

KRIGING  SHUTDOWN  STEP  AT  1649 

AR  SHUTDOWN  STEP  AT  1649 

JJJ.VdJJ)*  1650  0.40 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 

0.40 

KRIGING  SHUTDOWN  STEP  AT  1650 

AR  SHUTDOWN  STEP  AT  1650 

JJJ,V(JJJ)=  1651  0.43 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 

0.43 

KRIGING  SHUTDOWN  STEP  AT  1651 

AR  SHUTDOWN  STEP  AT  1651 

JJJ,V(IJJ)=  1652  0.45 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 

0.45 

KRIGING  SHUTDOWN  STEP  AT  1652 

AR  SHUTDOWN  STEP  AT  1652 

JJJ,V(JJJ)=  1653  0.46 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 

0.46 

KRIGING  SHUTDOWN  STEP  AT  1653 

AR  SHUTDOWN  STEP  AT  1653 

JJJ,V(JJJ)=  1654  0.45 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 

0.45 

KRIGING  SHUTDOWN  STEP  AT  1654 

AR  SHUTDOWN  STEP  AT  1654 

JJJ,V(JJJ)=  1655  0.43 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 

0.43 

KRIGING  SHUTDOWN  STEP  AT  1655 

AR  SHUTDOWN  STEP  AT  1655 

JJJ,VaJJ)=  1656  0.40 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 

0.40 

KRIGING  SHUTDOWN  STEP  AT  1656 

AR  SHUTDOWN  STEP  AT  1656 

JJJ,V{JJJ)=  1657  0.37 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 

0.37 

KRIGING  SHUTDOWN  STEP  AT  1657 

AR  SHUTDOWN  STEP  AT  1657 

JJJ,V(JJJ)=  1658  0.33 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 

0.33 

KRIGING  SHUTDOWN  STEP  AT  1658 

AR  SHUTDOWN  STEP  AT  1658 

JJJ,V(JJJ)=  1659  0.30 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 

0.30 

KRIGING  SHUTDOWN  STEP  AT  1659 

AR  SHUTDOWN  STEP  AT  1659 

JJJ,V(JJJ)=  1660  0.29 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 

0.29 

KRIGING  SHUTDOWN  STEP  AT  1660 

AR  SHUTDOWN  STEP  AT  1660 

JJJ,V(JJJ)=  1661  0.30 
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KRIGING  SHUTDOWN  STEP  AT  1663 
AR  SHUTDOWN  STEP  AT  1663 
JJJ,V(JJJ)=  1664  0.35 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE=  0.35 
KRIGING  SHUTDOWN  STEP  AT  1664 
AR  SHUTDOWN  STEP  AT  1664 
JJJ,V(JJJ)=  1665  0.36 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE  =  0.36 

KRIGING  SHUTDOWN  STEP  AT  1665 
AR  SHUTDOWN  STEP  AT  1665 
JJJ.V(IJJ)=  1666  0.36 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE=  0.36 
KRIGING  SHUTDOWN  STEP  AT  1666 
AR  SHUTDOWN  STEP  AT  1666 
JJJ,V(JJI)=  1667  0.35 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE=  0.35 
KRIGING  SHUTDOWN  STEP  AT  1667 
AR  SHUTDOWN  STEP  AT  1667 
JJJ,V(JJJ)=  1668  0.33 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE=  0.33 
KRIGING  SHUTDOWN  STEP  AT  1668 
AR  SHUTDOWN  STEP  AT  1668 
JJJ,V(JJJ)=  1669  0.30 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE  =  0.30 

KRIGING  SHUTDOWN  STEP  AT  1669 
AR  SHUTDOWN  STEP  AT  1669 
JJJ,V(JJJ)=  1670  0.28 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE=  0.28 
KRIGING  SHUTDOWN  STEP  AT  1670 
AR  SHUTDOWN  STEP  AT  1670 
JJJ,V(JJJ)=  1671  0.27 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE=  0.27 
KRIGING  SHUTDOWN  STEP  AT  1671 
AR  SHUTDOWN  STEP  AT  1671 
JJJ,V(JJJ)=  1672  0.26 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE=  0.26 
KRIGING  SHUTDOWN  STEP  AT  1672 
AR  SHUTDOWN  STEP  AT  1672 
JJJ,V(JJJ)=  1673  0.28 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE=  0.28 
KRIGING  SHUTDOWN  STEP  AT  1673 
AR  SHUTDOWN  STEP  AT  1673 
JJJ,V(JJJ)=  1674  0.31 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE=  0.31 
KRIGING  SHUTDOWN  STEP  AT  1674 
AR  SHUTDOWN  STEP  AT  1674 
JJJ,V(JJJ)=  1675  0.35 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE=  0.35 
KRIGING  SHUTDOWN  STEP  AT  1675 
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JJJ,V(IJJ)=  1678  0.44 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE = 

0.44 

KRIGING  SHUTDOWN  STEP  AT  1678 

AR  SHUTDOWN  STEP  AT  1678 

JJJ,V(JJJ)=  1679  0.43 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE  = 

0.43 

KRiGING  SHUTDOWN  STEP  AT  1679 

AR  SHUTDOWN  STEP  AT  1679 

JJJ,VaJJ)=  1680  0.39 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE = 

0.39 

AR  SHUTDOWN  STEP  AT  1680 

JJJ,V(JJJ)=  1681  0.35 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 

0.35 

KRIGING  SHUTDOWN  STEP  AT  1681 

AR  SHUTDOWN  STEP  AT  1681 

JJJ,V(JJJ)=  1682  0.29 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 

0.29 

KRIGING  SHUTDOWN  STEP  AT  1682 

AR  SHUTDOWN  STEP  AT  1682 

JJJ,V(JJJ)=  1683  0.25 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 

0.25 

KRIGING  SHUTDOWN  STEP  AT  1683 

AR  SHUTDOWN  STEP  AT  1683 

JJJ,V(JJJ)=  1684  0.23 

JJJ,V(JJJ)*  1685  0.23 

JJJ,V(JJJ)=  1686  0,24 

JJJ,V(JJJ)=  1687  0.27 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 

0.27 

AR  SHUTDOWN  STEP  AT  1687 

JJJ,V(JJJ)=  1688  0.30 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 

0.30 

KRIGING  SHUTDOWN  STEP  AT  1688 

AR  SHUTDOWN  STEP  AT  1688 

JJJ,V(JJJ)=  1689  0.34 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 

0.34 

KRIGING  SHUTDOWN  STEP  AT  1689 

AR  SHUTDOWN  STEP  AT  1689 

JJJ,V(JJJ)=  1690  0.35 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 

0.35 

KRIGING  SHUTDOWN  STEP  AT  1690 

AR  SHUTDOWN  STEP  AT  1690 

JJJ,V(JJJ)=  1691  0.35 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 

0.35 

KRIGING  SHUTDOWN  STEP  AT  1691 

AR  SHUTDOWN  STEP  AT  1691 

JJJ,V(JJJ)=  1692  0.34 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 

0.34 

KRIGING  SHUTDOWN  STEP  AT  1692 

AR  SHUTDOWN  STEP  AT  1692 
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EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE = 
1699  0.29 

0.25 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 
JJJ.V(JJJ)=  1700  0.33 

0.29 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 
JJJ,V(JJJ)*  1701  0.37 

0.33 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 
JJJ,V(JJJ)=  1702  0.39 

0.37 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 
JJJ,V(JJJ)=  1703  0.40 

0.39 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 
JJJ,V(JJJ)=  1704  0.38 

0.40 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 
JJJ,V(JJJ)=  1705  0.34 

0.38 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 
JJJ,V(JJJ)=  1706  0.29 

0.34 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 
JJJ,V(JJJ)=  1707  0.25 

0.29 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 
JJJ.V(JJJ)=  1708  0.23 

JJJ,VaJJ)=  1709  0.22 

JJJ,V(JJJ)=  1710  0.23 

JJJ,V(JJJ)=  1711  0.27 

0.25 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 
JJJ,V(JJJ)=  1712  0.32 

0.27 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 
JJJ,V(JJJ)=  1713  0.37 

0.32 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE* 
JJJ,V(JJJ)=  1714  0.40 

0.37 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE* 
JJJ,V(JJJ)=  1715  0.41 

0.40 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE* 
JJJ,V(JJJ)*  1716  0.40 

0.41 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE* 
JJJ,V(JJJ)*  1717  0.37 

0.40 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE* 
JJJ,V(JJJ)=  1718  0.35 

0.37 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE* 
JJJ,V(JJJ)*  1719  0.33 

0.35 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE* 
JJJ,V(JJJ)=  1720  0.31 

0.33 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE* 
JJJ,V(JJJ)=  1721  0.31 

0.31 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE* 
KRIGING  SHUTDOWN  STEP  AT  1721 

JJJ,V(JJJ)*  1722  0.32 

0.31 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE* 
KRIGING  SHUTDOWN  STEP  AT  1722 

JJJ,V(JJJ)=  1723  0.33 

0.32 
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KRIGING  SHUTDOWN  STEP  AT  1725 
AR  SHUTDOWN  STEP  AT  1725 
JJJ,V(JJJ)=  1726  0.40 

EXCEEDANCE  OF  THRESHOLD  FOUND.  VALUE=  0.40 
KRIGING  SHUTDOWN  STEP  AT  1726 
AR  SHUTDOWN  STEP  AT  1726 
JJJ,V(JJJ)=  1727  0.41 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE=  0  41 
KRIGING  SHUTDOWN  STEP  AT  1727 
AR  SHUTDOWN  STEP  AT  1727 
JJJ.V(JJJ)=  1728  0.40 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE=  0.40 
KRIGING  SHUTDOWN  STEP  AT  1728 
AR  SHUTDOWN  STEP  AT  1728 
JJJ,V(JJJ)=  1729  0.36 

EXCEEDANCE  OF  THRESHOLD  FOUND.  VALUE=  0.36 
KRIGING  SHUTDOWN  STEP  AT  1729 
AR  SHUTDOWN  STEP  AT  1729 
JJJ,V(JJJ)=  1730  0.31 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE=  0.3 1 
AR  SHUTDOWN  STEP  AT  1730 
JJJ,V(JJJ)=  1731  0.25 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE=  0.25 

JJJ,V(JJJ)=  1732  0.21 

nj,V(JJJ)=  1733  0.18 

JJJ,V(JJJ)=  1734  0.17 

JJJ,V(JJJ)=  1735  0.19 

JJJ,V(JJJ)=  1736  0.23 

JJJ,V(JJJ)=  1737  0.27 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE=  0.27 
KRIGING  SHUTDOWN  STEP  AT  1737 
AR  SHUTDOWN  STEP  AT  1737 
JJJ,V(JJJ)=  1738  0.30 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE=  0.30 
KRIGING  SHUTDOWN  STEP  AT  1738 
AR  SHUTDOWN  STEP  AT  1738 
JJJ,V(JJJ)=  1739  0.30 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE  =  0.30 

KRIGING  SHUTDOWN  STEP  AT  1739 
AR  SHUTDOWN  STEP  AT  1739 
JJJ,V(JJJ)=  1740  0.29 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE  =  0.29 

AR  SHUTDOWN  STEP  AT  1740 
JJJ,V(JJJ)=  1741  0.26 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE=  0.26 

AR  SHUTDOWN  STEP  AT  1741 

JJJ,V(JJJ)=  1742  0.22 

JJJ,V(JJJ)=  1743  0.18 

JJJ.V(JJJ)=  1744  0.15 
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EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 

JJJ,V(JJJ)=  1753 

0.23 

JJJ,V(JJJ)=  1754 

0  19 

JJJ,VaJJ)=  1755 

0.14 

HJ,V(JJJ)=  1756 

0.09 

JJJ.V(JJJ)=  1757 

0.06 

JJJ,V(JJJ)=  1758 

0.04 

JJJ,V(JJJ)=  1759 

0.04 

JJJ,V(JJJ)=  1760 

0.06 

JJJ,V(JJJ)=  1761 

0.09 

JJJ,V(JJJ)=  1762 

0.12 

JJJ,VaJJ)=  1763 

0.13 

JJJ,V(JJJ)=  1764 

0.13 

JJJ,V(JJJ)=  1765 

0.11 

JJJ,V(JJJ)=  1766 

0.09 

JJJ,V(JJJ)=  1767 

0.06 

JJJ,V(JJJ)=  1768 

0.04 

JJJ,V(JJJ)=  1769 

0.04 

JJJ,V(JJJ)=  1770 

0.05 

JJJ,V(JJJ)=  1771 

0.07 

JJJ,V(JJJ)=  1772 

0.09 

JJJ,V(JJJ)=  1773 

0.12 

JJJ,V(JJJ)=  1774 

0.15 

JJJ,V(JJJ)=  1775 

0.18 

JJJ,V(JJJ)=  1776 

0.18 

JJJ,V(JJJ)=  1777 

0.18 

JJJ,V(JJJ)=  1778 

0.16 

JJJ,V(JJJ)=  1779 

0.13 

1780 

0.10 

JJJ,V(JJJ)=  1781 

0.06 

JJJ,V(JJJ)=  1782 

0.04 

JJJ,V(JJJ)=  1783 

0.03 

JJJ,V(JJJ)=  1784 

0.04 

JJJ,V(JJJ)=  1785 

0.06 

JJJ,V(JJJ)=  1786 

0.09 

JJJ,V(JJJ)=  1787 

0.11 

JJJ,V(JJJ)=  1788 

0.13 

JJJ,V(JJJ)=  1789 

0.13 

JJJ,V(JJJ)=  1790 

0.12 

JJJ,V(JJJ)=  1791 

0.09 

JJJ,V(JJJ)=  1792 

0.06 

JJJ,V(JJJ)=  1793 

0.04 

JJJ,V(JJJ)=  1794 

0.03 

JJJ,V(JJJ)=  1795 

0.04 

JJJ,V(JJJ)=  1796 

0.07 

JJJ,V(JJJ)=  1797 

0.09 

JJJ,V(JJJ)=  1798 

0.12 

JJJ,V(JJJ)=  1799 

0.15 

JJJ,V(JJJ)=  1800 

0.17 

0.26 
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JJJ.V(JJJ)=  1810 

0.16 

JJJ,V(JJJ)=  1811 

0.18 

JJJ,V(JJJ)=  1812 

0.20 

JJJ,V(JJJ)=  1813 

0.20 

JJJ,V(JJJ)=  1814 

0.19 

JJJ,V(JJJ)=  1815 

0.16 

JJJ.V(JJJ)=  1816 

0.12 

JJJ.V(JJJ)=  1817 

0.08 

JJJ,V(JJJ)=  1818 

0.06 

1819 

0.04 

JJJ.V(JJJ)=  1820 

0.03 

JJJ,V(JJJ)=  1821 

0.04 

JJJ,V(JJJ)=  1822 

0.04 

JJJ,V(JJJ)=  1823 

0.05 

JJJ,V(JJJ)=  1824 

0.06 

JJJ,V(JJJ)=  1825 

0.06 

JJJ,V(JJJ)=  1826 

0.06 

JJJ,V(JJJ)=  1827 

0.04 

JJJ,V(JJJ)=  1828 

0.02 

JJJ.V(JJJ)=  1829 

0.00 

JJJ,V(JJJ)=  1830 

-0.02 

JJJ,V(JJJ)=  1831 

-0.04 

JJJ,V(JJJ)=  1832 

-0.04 

JJJ,V(JJJ)=  1833 

-0.04 

JJJ,V(JJJ)=  1834 

-0.03 

JJJ,V(JJJ)=  1835 

-0.01 

JJJ.V(JJJ)=  1836 

0.00 

JJJ,V(JJJ)=  1837 

0.02 

JJJ,V(JJJ)=  1838 

0.03 

JJJ,V(JJJ)=  1839 

0.03 

JJJ,V(JJJ)=  1840 

0.01 

JJJ,V(JJJ)=  1841 

-0.01 

JJJ,V(JJJ)=  1842 

-0.03 

JJJ,V(JJJ)=  1843 

-0.04 

JJT,V(JJJ)=  1844 

-0.04 

JJJ,V(JJJ)=  1845 

-0.04 

JJJ,V(JJJ)=  1846 

-0.03 

JJJ,V(JJJ)=  1847 

-0.01 

JJJ,V(JJJ)=  1848 

0.00 

JJJ,V(JJJ)=  1849 

0.01 

JJJ,V(JJJ)=  1850 

0.02 

JJJ,V(JJJ)=  1851 

0.02 

JJJ,V(JJJ)=  1852 

0.01 

JJJ,V(JJJ)=  1853 

0.00 

JJJ,V(JJJ)=  1854 

-0.01 

JJJ,V(JJJ)=  1855 

-0.02 

JJJ,V(JJJ)=  1856 

-0.02 

JJJ,V(JJJ)=  1857 

-0.02 

JJJ,V(JJJ)=  1858 

0.00 
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JJJ,V(JJJ)=  1868 

0.16 

JJJ,V(JJJ)=  1869 

0.17 

JJJ.V(JJJ)=  1870 

0.19 

JJJ.V(JJJ)=  1871 

0.20 

JJJ,V0JJ)=  1872 

0.21 

JJJ,V(JJJ)=  1873 

0.22 

JJJ,V(JJJ)=  1874 

0.23 

JJJ.V(JJJ)=  1875 

0.23 

JIJ,V(JJJ)=  1876 

0.23 

JJJ.V(JJJ)=  1877 

0.21 

JJJ,V(JJJ)=  1878 

0.20 

JJJ.V(JJJ)=  1879 

0.18 

JJJ,V(JJJ)=  1880 

0.17 

JJJ.V(JJJ)=  1881 

0.16 

JJJ,V(JJJ)=  1882 

0.15 

JJJ,V(JJJ)=  1883 

0.15 

JJJ,V(JJJ)=  1884 

0.15 

JJJ,V(JJJ)=  1885 

0.16 

JJJ,V(JJJ)=  1886 

0.17 

JJJ,V(JJJ)=  1887 

0.18 

JJJ,V(JJJ)=  1888 

0.17 

JJI.V(JJJ)=  1889 

0.16 

JJJ,V(JJJ)=  1890 

0.14 

JJJ,V(jJJ)=  1891 

0.l2 

JJJ,V(JJJ)=  1892 

0.10 

JJJ,V(JJJ)=  1893 

0.09 

JJJ,V(JJJ)=  1894 

0.08 

JJJ,V(JJJ)=  1895 

0.07 

JJJ,V(JJJ)=  1896 

0.07 

JJJ,V(JJJ)=  1897 

0.07 

JJJ,V(JJJ)=  1898 

0.07 

JJJ,V(JJJ)=  1899 

0.08 

JJJ,V(JJJ)=  1900 

0.09 

JJJ,V(JJJ)=  1901 

0.09 

JJJ,V(JJJ)=  1902 

0.09 

JJJ,V(JJJ)=  1903 

0.08 

JJJ,V(JJJ)=  1904 

0.06 

JJJ,V(JJJ)=  1905 

0.04 

JJJ,V{JJJ)=  1906 

0.02 

JJJ,V(JJJ)=  1907 

0.01 

JJJ,V(JJJ)=  1908 

0.00 

JJJ,V(JJJ)=  1909 

0.00 

JJJ,V(JJJ)=  1910 

0.00 

JJJ,V(JJJ)=  1911 

0.01 

JJJ,V(JJJ)=  1912 

0.01 

JJJ,V(JJJ)=  1913 

0.01 

JJJ,V(JJJ)=  1914 

0.01 

JJJ,V(JJJ)=  1915 

000 

JJJ,V(JJJ)=  1916 

-0.02 
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JJJ.V(JJJ)=  1926  0.11 

1927  0.13 

1928  0.14 

JJJ,VaJJ)*  1929  0.15 

JJJ,V(JJJ)«  1930  0.16 

1931  0.18 

JJJ.VdJJ)*  1932  0.19 

1933  0.21 

JJJ,V(JJJ)=  1934  0.24 

1935  0.27 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 

0.27 

1936  0.30 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 

0.30 

JJJ,V(JJJ)=  1937  0.31 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 

0.31 

JJJ,V(JJJ)=  1938  0.32 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 

0.32 

JJJ,V(JJJ)=  1939  0.32 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 

0.32 

JJJ,V(JJJ)=  1940  0.31 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE = 

0.31 

1941  0.29 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 

0.29 

JJJ,V(JJJ)*  1942  0.27 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE* 

0.27 

JJJ,V(JJJ)«  1943  0.26 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE* 

0.26 

JJJ,V(JJJ)*  1944  0.25 

JJJ,V(JJJ)=  1945  0.24 

JJJ,V(JJJ)=  1946  0.24 

JJJ,V(JJJ)=  1947  0.25 

JJJ,V(JJJ)=  1948  0.26 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE* 

0.26 

JJJ,V(JJJ)=  1949  0.27 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE* 

0.27 

JJJ,V(JJJ)=  1950  0.29 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE* 

0.29 

JJJ,V(JJJ)=  1951  0.29 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE* 

0.29 

JJJ,V(JJJ)=  1952  0.29 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE* 

0.29 

KRIGING  SHUTDOWN  STEP  AT  1952 

JJJ,V(JJJ)*  1953  0.30 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE* 

0.30 

KRIGING  SHUTDOWN  STEP  AT  1953 

AR  SHUTDOWN  STEP  AT  1953 

JJJ,V(JJJ)*  1954  0.30 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE* 

0.30 

KRIGING  SHUTDOWN  STEP  AT  1954 
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JJJ,V(JJJ)=  1957  0.30 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE  =  0.30 

KRIGING  SHUTDOWN  STEP  AT  1957 
AR  SHUTDOWN  STEP  AT  1957 
JJJ.V(JJJ)=  1958  0.30 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE=  0.30 
KRIGING  SHUTDOWN  STEP  AT  1958 
AR  SHUTDOWN  STEP  AT  1958 
IJJ,V(JJJ)=  1959  0.31 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE=  0.31 
KRIGING  SHUTDOWN  STEP  AT  1959 
AR  SHUTDOWN  STEP  AT  1959 
JJJ,V(JJJ)=  1960  0.32 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE  =  0.32 

KRIGING  SHUTDOWN  STEP  AT  1960 
AR  SHUTDOWN  STEP  AT  1960 
JJJ,V(JJJ)=  1961  0.33 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE=  0.33 
KRIGING  SHUTDOWN  STEP  AT  1961 
AR  SHUTDOWN  STEP  AT  1961 
JJJ,V(JJJ)=  1962  0.34 

EXCEEDANCE  OF  THRESHOLD  FOUND.  VALUE=  0.34 
KRIGING  SHUTDOWN  STEP  AT  1962 
AR  SHUTDOWN  STEP  AT  1962 
JJJ,V(JJJ)=  1963  0.33 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE=  0.33 
KRIGING  SHUTDOWN  STEP  AT  1963 
AR  SHUTDOWN  STEP  AT  1963 
JJJ,V(JJJ)=  1964  0.31 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE=  0.31 
KRIGING  SHUTDOWN  STEP  AT  1964 
AR  SHUTDOWN  STEP  AT  1964 
JJJ,V(JJJ)=  1965  0.28 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE=  0.28 
KRIGING  SHUTDOWN  STEP  AT  1965 
AR  SHUTDOWN  STEP  AT  1965 


JJJ,V(JJJ)= 

1966 

0.23 

JJJ,V(JJJ)= 

1967 

0.19 

JJJ,V(JJJ)= 

1968 

0.15 

JJJ,V(JJJ)* 

1969 

0.12 

JJJ,V(JJJ)= 

1970 

0.09 

JJJ,V(JJJ)= 

1971 

0.07 

JJJ,V(JJJ)= 

1972 

0.07 

JJJ,V(JJJ)= 

1973 

0.07 

JJJ,V(JJJ)= 

1974 

0.08 

JJJ,V(JJJ)* 

1975 

0.09 

JJJ,V(JJJ)= 

1976 

0.09 

JJJ,V(JJJ)= 

1977 

0.08 

JJJ,V(JJJ)= 

1978 

0.06 
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JJJ,V(JJJ)=  1988 

0.19 

JJJ,V(JJJ)=  1989 

0.18 

JJJ,V(JJJ)=  1990 

0.15 

1991 

0.13 

JJJ,V(JJJ)=  1992 

0.10 

1993 

0.07 

JJJ,V(JJJ)=  1994 

0.05 

JJJ.V(JJJ)=  1995 

0.04 

JJJ,V(JJJ)=  1996 

0.05 

JJJ,V(JJJ)=  1997 

0.07 

JJJ.V(JJJ)=  1998 

0.10 

JJJ,V(JJJ)=  1999 

0.13 

JJJ,V(JJJ)=  2000 

0.16 

JJJ,V(JJJ)=  2001 

0.18 

JJJ,V(JJJ)=  2002 

0.19 

JJJ,V(JJJ)=  2003 

0.19 

JJJ,V(JJJ)=  2004 

0.18 

JJJ,V(JJJ)=  2005 

0.18 

JJJ,V(JJJ)=  2006 

0.18 

JJJ,V(JJJ)=  2007 

0.19 

JJJ,V(JJJ)=  2008 

0.21 

JJJ,V(JJJ)=  2009 

0.24 

JJJ.V(JJJ)=  2010 

0.28 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE = 

0.28 

2011 

0.32 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE = 

0.32 

JJJ,V(JJJ)=  2012 

0.35 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 

0.35 

JJJ,V(JJJ)=  2013 

0.35 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 

0.35 

JJJ,V(JJJ)=  2014 

0.34 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE = 

0.34 

JJJ,V(JJJ)=  2015 

0.32 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE = 

0.32 

JJJ,V(JJJ)=  2016 

0.28 

EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE= 

0.28 

JJJ,V(JJJ)=  2017 

0.24 

JJJ,V(JJJ)=  2018 

0.19 

JJJ,V(JJJ)=  2019 

0.16 

JJJ,V(JJJ)=  2020 

0.14 

JJJ,V(JJJ)=  2021 

0.14 

JJJ,V(JJJ)=  2022 

0.15 

JJJ,V(JJJ)=  2023 

0.17 

JJJ,V(JJJ)=  2024 

0.20 

JJJ,V(JJJ)=  2025 

0.21 

JJJ,V(JJJ)=  2026 

0.21 

JJJ,V(JJJ)=  2027 

0.20 

JJJ,V(JJJ)=  2028 

0.17 

JJJ,V(JJJ)=  2029 

0.14 
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EXCEEDANCE  OF  THRESHOLD  FOUND,  VALUE*  0.26 

JJJ,V(JJJ)=  2038 

0.24 

JJJ,V(JJJ)=  2039 

0.20 

JJJ.VaJJ)=  2040 

0.15 

JJJ.VaJJ)*  2041 

0.09 

JJJ.V(JJJ)=  2042 

0.04 

JJJ,V(JJJ)=  2043 

-0.01 

JJJ.V(JJJ)=  2044 

-0.03 

JJJ,V(JJJ)=  2045 

-0.04 

JJJ.V(JJJ)=  2046 

-0.02 

JJJ,V(JJJ)=  2047 

0.00 

JJJ.V(JJJ)=  2048 

0.02 

JJJ,V(JJJ)=  2049 

0.04 

JJJ.V(JJJ)=  2050 

0.04 

JJJ,V(IJJ)=  2051 

0.03 

JJJ,V(JJJ)=  2052 

0.00 

JJJ,V(JJJ)=  2053 

-0.02 

JJJ,V(JJJ)=  2054 

-0.04 

JJJ.V(JJJ)=  2055 

-0.06 

JJJ,V(JJJ)=  2056 

-0.05 

JJJ,V(JJJ)=  2057 

-0.02 

JJJ.V(JJJ)=  2058 

0.03 

JJJ,V(JJJ)=  2059 

0.08 

JJJ,VaJJ)=  2060 

0.13 

JJJ,V(IJJ)=  2061 

0.16 

2062 

0.17 

JJJ,V(JJJ)=  2063 

0.15 

JJJ,V(JJJ)=  2064 

0.12 

JJJ.V(JJJ)=  2065 

0.07 

JJJ,V(JJJ)=  2066 

0.03 

JJJ,V(JJJ)*  2067 

-0.01 

JJJ,V(JJJ)=  2068 

-0.04 

JJJ,V(JJJ)=  2069 

-0.05 

JJJ,V(JJJ)=  2070 

-0.05 

JJJ,V(JJJ)=  2071 

-0.03 

JJJ,V(JJJ)=  2072 

0.00 

2073 

0.02 

JJJ,V(JJJ)=  2074 

0.03 

JJJ,V(JJJ)=  2075 

0.02 

JJJ,V(JJJ)=  2076 

-0.01 

JJJ.VaJJ)*  2077 

-0,04 

JJJ,V(JJJ)=  2078 

-0.07 

JJJ,V(JJJ)*  2079 

-0.09 

JJJ,V(JJJ)*  2080 

-0.09 

JJJ,V(JJJ)=  2081 

-0.08 

2082 

-0.04 

JJJ,V(JJJ)*  2083 

0.01 

JJJ,V(JJJ)=  2084 

0.07 

JJJ,V(JJJ)=  2085 

0.12 
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JJJ,V(JJJ)=  2095 

-0.03 

JJJ,V(JJJ)=  2096 

-0.01 

JJJ,V(JJJ)=  2097 

0.00 

JJJ,V(J)J)=  2098 

0.01 

JJJ,V(JJJ)=  2099 

0.00 

JJJ.V(JJJ)=  2100 

-0.03 

JJJ,V(JJJ)=  2101 

-0.07 

JJJ,V(JJJ)=  2102 

-0.10 

JJJ,V(JJJ)=  2103 

-0.13 

JJJ,V(JJJ)=  2104 

-0.15 

JJJ.VaiJ)=  2105 

-0.14 

JJJ.VaJJ)=  2106 

-0.12 

JJJ,V(JJJ)=  2107 

-0.07 

JJJ,V(JJJ)=  2108 

-0.01 

JJJ,V(JJJ)=  2109 

0.05 

JJj,V(JJJ)=  2110 

0.09 

JJJ,V(JJJ)=  2111 

0.12 

JJJ,V(JJJ)=  2112 

0.12 

JJJ,V(JJJ)=  2113 

0.11 

JJJ.V(JJJ)=  2114 

0.08 

JJJ.V(JJJ)=  2115 

0.04 

JJJ,V(JJJ)=  2116 

0.01 

JJJ,V(JJJ)=  2117 

-0.02 

JJJ.V(JJJ)=  2118 

-0.04 

JJJ,V(JJJ)=  2119 

-0.04 

JJJ,V(JJJ)=  2120 

-0.03 

JJJ.V(JJJ)=  2121 

-0.01 

JJJ.VaiJ)*  2122 

0.00 

2123 

0.01 

JJJ,V(JJJ)*  2124 

-0.01 

JJJ,V(JJJ)=  2125 

-0.04 

JJJ,V(JJJ)=  2126 

-0.07 

JJJ,V(JJJ)=  2127 

-0.11 

JJJ,V(JJJ)=  2128 

-0.14 

JJJ,V(JJJ)=  2129 

-0.17 

JJJ,V(JJJ)=  2130 

-0.17 

JJJ,V(JJJ)=  2131 

-0.14 

JJJ,V(JJJ)=  2132 

-0.09 

JJJ,V(JJJ)=  2133 

-0.04 

JJJ,V(JJJ)=  2134 

0.01 

2135 

0.04 

JJJ,V(JJJ)=  2136 

0.06 
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